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A COMPUTER PROGRAM FOR CALCULATING INVISCID, ADIABATIC FLOW 
ABOUT BLUNT BODIES TRAVELING AT SUPERSONIC AND 
HYPERSONIC SPEEDS AT ANGLE OF ATTACK 

By Richard W. Barnwell and Ruby M. Davis 
Langley Research Center 

SUMMARY 

This paper is a user's manual for a computer program which calculates inviscid 
plane, axisymmetric, and three-dimensional flow about blunt bodies traveling at super- 
sonic and hypersonic speeds in a uniform free stream. An exact time-dependent finite- 
difference method of second-order accuracy which is described in NASA TN D-6283 is 
used. The bodies which can be treated include plane and axisymmetric bodies with sharp 
shoulders and smooth nonaxisymmetric bodies. Equilibrium-air and perfect-gas thermo- 
dynamic models can be used. A procedure for approximating equilibrium gases with the 
perfect-gas model is described. The results of the program include the shock-wave loca- 
tion and the flow properties at a number of grid points on the body surface, on the shock 
wave, and in the region between the body and shock. The paper contains descriptions of 
the input and output, a Fortran IV program listing, and an input and output for a sample 
case. 


INTRODUCTION 

This paper is a user's manual for a computer program which calculates inviscid 
plane, axisymmetric, and three-dimensional flow about blunt bodies traveling at super- 
sonic and hypersonic speeds in a uniform free stream. The bodies which can be treated 
include plane and axisymmetric bodies with sharp shoulders and smooth nonaxisymmetric 
bodies. Equilibrium-air and perfect-gas thermodynamic models can be used. A curve 
fit is used in connection with the equilibrium-air model. A procedure for approximating 
equilibrium gases with the perfect-gas model is described. The results of the program 
include shock-wave location and the flow properties at a number of grid points on the body 
surface, on the shock wave, and in the region between the shock and body. 

The numerical method which is employed is an exact time- dependent finite- 
difference method of second-order accuracy and is described in detail in references 1 
and 2. Results for steady flow are obtained with time- dependent methods by applying 
them to the solution of initial-value problems. An approximate method is employed to 



determine a starting solution at the initial time. The time-dependent method is used to 
calculate the solution at subsequent time steps. Results for steady flow are obtained 
after many time steps when the difference in the solutions at consecutive time steps is 
sufficiently small. 

A Fortran IV program listing is given in appendix A, and a sample input and output 
are given in appendixes B and C, respectively. The minimum number of computer stor- 
age locations required to execute the program on the CDC 6600 computer is 70 OOOg. 

SYMBOLS 

speed of sound 

bluntness of initial shock shape 
time-step coefficient 

length characteristic of subsonic region of blunt-body flow field 
internal energy (heat of formation chosen so that e = 0 at T = 0° K) 
function plotted in figure 9 

function which accounts for nonorthogonal nature of computational coordinate 
system 

static enthalpy (heat of formation chosen so that h = 0 at T = 0° K) 

index which indicates method of solution at points of curvature discontinuity 
on body 

number: j = 0 for plane and axisymmetric flow fields; j = 1 for three- 
dimensional flow fields 

curvature of ith segment of reference-surface generator 
kj,k 2 ,k 3 ,k 4 values of l where curvature of reference surface changes 

Z,m,n grid indices for X-, Y-, and <f> -coordinates 


a 

B s 

C T 

d 

e 

f( e / e max) 

g 

h 

I 

j 

Ki 
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M 

Ni 

P 

R 

R b,i 

R s 

r 

Si 

T 

u,v,w 

V 

X 


X 

x o 

x l’ x 2’ x 3 

Y 


yb 

ybo 


Mach number 

number of mesh spacings AX on ith reference-surface segment 
static pressure 

perpendicular distance from axis of initial shock shape 
radius of curvature of ith body segment 

♦ 

nose radius of initial shock shape 

perpendicular distance from axis of reference surface 

arc length of ith body segment 

absolute temperature 

velocity components 

magnitude of total velocity 

distance along reference-surface generator from axis 

distance along axis of reference surface from intersection with body surface 
distance shown in figure 7 
Cartesian coordinates 

quantity proportional to distance from body along lines perpendicular to ref- 
erence surface: Y = 0 at body; Y = 1 at shock 

distance from reference surface to body along normals to reference surface 

constant value of y^ 

angle of attack 
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r 

y 

AX,AY,A4> 
At 
5 * 
e . 

I 

e 

\ 

p 

T 

$ 

X 

Xo 

<jO 

Subscripts: 


ratio of specific heats 

ratio of static enthalpy h to internal energy e (heat of formation chosen 
so that h = e = 0 at T = 0° K) 

mesh spacings for X-, Y-, and $ -coordinates 

mesh spacing for time given by equation (8) 

distances between shock and body along normal to reference surface 
damping coefficient 

reference length used to normalize distance 

angle between normal to reference surface and axis of reference surface 

scale factor for X-coordinate 

density 

time 

azimuthal angle 

distance along axis of initial shock from intersection with shock 
distance shown in figure 7 

distance from origin of xj,X 2 ,Xg coordinate system 


f largest value of grid indices Z,m,n 

max maximum value 

min minimum value 

s • quantity evaluated immediately behind normal shock 
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t 


quantity evaluated at stagnation point 


00 quantity evaluated in free stream 

Barred quantities are dimensional, and unbarred quantities are nondimensionaL 

DISCUSSION 
Type of Bodies Treated 

Supersonic flow fields about two classes of two- and three-dimensional bodies can 
be calculated with this program. Class I bodies are symmetric about an axis and can 
have discontinuities in both slope and curvature. Class II bodies must have continuous 
slope and curvature. Two-dimensional class II bodies are not required to be symmetric, 
but three-dimensional class II bodies must be symmetric with respect to a plane. 

The generator of a class I body surface can have a maximum of four points of dis- 
continuity in curvature or three points of discontinuity in curvature and one point of dis- 
continuity in slope. The portion of the generator of a class I body immediately upstream 
of a discontinuity in slope where the flow becomes sonic (point A in fig. 1(a)) must be a 
segment of constant curvature (circular arc or straight-line segment). Also, the portions 
of the generator on either side of a discontinuity in curvature where the flow-property dis- 
tributions have discontinuities in slope (point B in fig. 1(a)) must be segments of constant 
curvature. These requirements occur because of the way in which the program is for- 
mulated. The curvature of other portions of the surface generator of a class I body can 
vary continuously. 

It is required that the flow fields about three-dimensional class II bodies have one 
plane of symmetry. Therefore, the free-stream velocity vector must be parallel to the 




(b) Class II. 


Figure 1. - Types of bodies treated. 
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plane of symmetry of the body. It should be noted that three-dimensional class I bodies 
are axisymmetric so that the flow fields about them always have one plane of symmetry. 

Nondimensional Quantities 

In this computer program, all lengths are normalized with some reference length £ 
which is chosen by the user. The velocity components are made nondimensional with the 
magnitude of the free-stream velocity V^; the density, with the free-stream density p^; 

_ — o 

and the pressure, with the product . The speed of sound is made nondimensional 

— — 2 
with V^, the enthalpy and internal energy are normalized with , and the time is non- 

dimensionalized with the quantity ^V m . It should be noted that the nondimensional free- 

stream pressure is 

p 1 

p.— 9 (1) 

PooVoo Voo M J 

where y ^ is the free-stream ratio of specific heats, and M w is the free-stream Mach 
number. 


Basic Geometry 

Coordinate system .- A nonorthogonal computational coordinate system is used. 

This coordinate system is oriented with respect to a reference surface which is sym- 
metric about an axis and which has a generator composed of segments of constant curva- 
ture (circular arcs and straight-line segments). The computer program is set up to use 
a maximum of five segments. An example of the computational grid for a three- 
dimensional flow field is shown in figure 2. 

The region on one side of the symmetry plane is divided with a number of evenly 
spaced half-planes which pass through the axis. The angle $, which is used to identify 
the various half-planes, is measured from the leeward side of the symmetry plane. Each 
of the half-planes is subdivided with lines of constant X and Y. The coordinate X is 
the distance from the axis along the reference-surface generator. The grid lines of con- 
stant X are straight lines which are perpendicular to the reference surface and which 
are spaced so that the distance along the reference-surface generator between adjacent 
lines is the constant AX. The line X = X max must be located so that it is downstream 
of the sonic line in each half-plane of constant <£. The quantity Y is a normalized 
variable which is measured along the lines normal to the reference surface and has the 
value 0 at the body surface and 1 at the shock wave. 

As shown in figure 3, the distance from the body to the shock wave along the lines of 
constant X in each plane of constant $ is designated as 8, and the distance from the 
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reference surface to the body surface along these lines is called y^. The shock-layer 
thickness 5 is a function of time r as well as of X and $ because the position of 
the shock wave adjusts with time during the course of the computation. In general, the 
distance y^ is a function of X and <&. The velocity components which are used are 
oriented with respect to the reference surface rather than the body surface. The com- 
ponent u lies in the plane of constant $ and is perpendicular to the normal to the ref- 
erence surface; the component v is in the direction of the normal to the reference sur- 
face; and the component w is perpendicular to the plane of constant <3>. 



Figure J. - Velocity components, auxiliary coordinates, and geometric parameters. 


An auxiliary coordinate system x,r is used to locate points in the planes of con- 
stant <f>. In this system, which is shown in figure 3, the coordinate r is the perpen- 
dicular distance from the axis of the reference surface, and the coordinate x is the dis- 
tance from the body along the axis in the upstream direction. The angle between the 
normal to the reference surface and the direction of the axis is denoted by 6 . 

The points of discontinuity in curvature and slope on the body-surface generator of 
class I bodies must lie on normals to the reference surface which pass through the points 
of discontinuity on the reference-surface generator. Because of the way in which the 
computer program is set up, the distance yj 3 (X,<i>) from the reference surface to the 
body surface along normals to the reference surface must have a constant value y^ Q for 
the segment of the body surface located immediately upstream of a sharp sonic corner 
(point A in fig. 1(a)) and for the segments on either side of a point of discontinuity in cur- 
vature on the body-surface generator where the flow-property distributions have discon- 
tinuities in slope (point B in fig. 1(a)). It should be noted that when y^(X,$) = y^ Q , 
straight-line segments on the reference-surface generator are paired with parallel 
straight-line segments on the body-surface generator, circular arcs on the reference- 
surface generator are paired with concentric circular arcs on the body-surface generator, 
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and a circular arc on the reference-surface generator with a radius of y^ 0 subtends a 
sharp corner on the body-surface generator. Also, when y^(X,<t>) = y^ Q , the velocity 
components are oriented with respect to the body surface as well as the reference surface. 

Grid indices .- The grid indices which are used for the X-, Y-, and 4? -coordinates 
are Z, m, and n, respectively. The following grid indices are used for the boundary 


locations : 

Axis (X = 0) Z = 1 

Line X = Xm ^y Z ~ Zf 

Body m =. 1 

Shock m = m-f 

Leeward side of symmetry plane (<i> = 0) n = 1 

Windward side of symmetry plane (4? = 180°) n = nf 


The curvature of the reference-surface generator can be discontinuous at four 
points. The indices of these points are Z = kj, k 2 , k 3 , and k 4 where kf < k 2 < k 3 < k,j. 
The grid indices for the body shown in figure 2 are given in figure 4. 


Half-plane of constant $ 



Windward 


Figure 4. - Grid indices for sharp- shouldered axisymmetric body. 

The values of the indices Zf, mj, nf, kj, k 2 , k 3 , and k^ must be supplied by 
the user. Because of the way in which the program is set up, the value of mj cannot be 
less than 3. For axisymmetric flow or symmetric two-dimensional flow, only one plane 
of constant is needed; therefore, nf can be set equal to 1. For nonsymmetric two- 
dimensional flow nf must have a value of 2. The program is written so that each 
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segment of the reference-surface generator except the first must have a length of not 
less than 3 AX. The first segment must have a length of at least 2 AX. Therefore, 
the indices of the points of discontinuity must satisfy the relations 

kj g 3 

ki = ki_i + 3 

The unused indices k^ should be set equal to some number greater than if. 

Geometry Input for Class I Bodies 

Consider class I bodies for which y b (X,4>) = y^. Let Sf and f (i = 1,. . .,5) 
be the arc lengths and radii of curvature of the segments of the body-surface generator. 
The curvature of each of the corresponding segments of the reference surface is deter- 
mined with the equation 


(i = 2,3,4) 


(2) 


Ki = 


R 


b,i + ^bo 


(i = 1,. • -,5) 


(3) 


The body shown in figure 4 is a class I body with a generator which, for the purposes of 
this program, is considered to be composed of three segments. The second segment is 
associated with the corner and has an arc length S 2 and a radius of curvature R b 2 
of zero. However, the ratio 2 is a ^ nite number equal to the angle between the 

lines normal to the upstream and downstream surfaces at the corner. In this program, 
the curvatures of convex and concave arcs are positive and negative, respectively. The 
radius of curvature of an arc has the same sign as the curvature of that arc. For exam- 
ple, the arcs on the generators of the reference surface and body surface shown in fig- 
ure 4 are convex and hence have positive curvatures and radii of curvature. All arc 
lengths are considered to be positive. 

All arcs of the reference-surface generator must be integral multiples of AX. If 
the body surface is located a constant distance y bo from the reference surface, the mesh 
spacing AX is related to the quantities Si, Rfc,^, and y bo by the equation 


AX = 


Sj K,i + ybo\ 

^i\ ^b,i / 


(4) 


where , is the number of mesh spacings associated with the ith segment of the body 
surface. If the body-surface generator has only one segment with a radius of curvature 
Rjj j, the values for the quantities Sj, Nj, and y bo can be chosen independently. If 
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the body-surface generator has two segments with radii of curvature R^i and R b 2 , 
the values of the quantities Sj, S 2 , Nj, and N 2 can be chosen independently, but the 
value of y^Q must satisfy the equation 


s ll 

(Vi + O 

1 S2 ( 

fe b ,2 + ybo) 

Ni 

\ R b,l / 

1 n 2 \ 

v R b ,2 j 


If the body-surface generator has more than two segments of constant curvature, all of 
the quantities R b Sj, and Nj (i = 1 , 2 , etc.) cannot be chosen independently if 
y]g(X,$) = y ^ 0 for all the segments. The index lq is related to the number Nj of 
mesh spacings AX on the ith arc by the equations 

k! = 1 + Nj_ 

k i = k i-l + N i (i = 2,3,4) 


Spherical cap- cylinder .- Consider the spherical cap-cylinder shown in figure 4. 

Let the scaling length be the radius of the cylinder r^ so that the nondimensional radius 

of the cylinder is 1, and assume that y b (X,4>) = y bo as shown. Let the nose radius have 
the value \/2. It follows that the arc length Sj has the value (n/4)\[2. The radius Rjj 2 
and arc length S 2 of the corner are both zero, but the ratio S2/Rb ; 2 is tf/4. The 
radius of the segment of the body-surface generator for the cylinder R^ 3 is °°, and 
the arc length S 3 is unspecified. From figure 4 it is seen that kj = 7 and k 2 = 10. 

It follows from equations (5) and (4), respectively, that y bo = ]/2 and AX = \/ 27 r/l 2 . 

From equation (3) it is found that the curvature for the segments of the reference-surface 
generator are Kj = |/2/4, K 2 = 1 / 2 / 2 , and K 3 = 0. 

Spherically blunted cone-cylinder .- Consider a spherically blunted cone-cylinder 
with a sharp corner at the junction of the cone and cylinder, a cone semiapex angle of 60°, 
and a ratio of nose radius to base radius equal to 0.5 as shown in figure 5. The scaling 
length is chosen to be the cylinder radius r^. The indices kj, k 2 , and k 3 are chosen 
to have values of 4, 10, and 13, respectively. Also, Sj = 77 / 12 ; Rj-,2 = °°; S 2 = 

R b 3 = 0; S 3 = 0; and R b 4 = °°. From equations (5) and (4), the values of the quantities 
y^Q and AX are found to be 


ybo = 


1/V3 

2 \ 77 


-0.32699334 
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k 3 = 13 


k,= 10 ✓ / 

2 A A0 / 
/ v 



Figure 5- - Reference surface for spherically blunted 
cone-cylinder with cone semiapex angle of 60 ° and 
ratio of nose radius to base radius equal to 0 . 5 . 


and 


AX = jL(l + 2y b J) = \/3/l2 * 0.14433757 


From equation (3) it is found that 

Kj = 2lL ~ 1.2091996 
3\/3 

K 2 = 0 

Ko =— |S — * 3.0581662 
3\j3 - 7 t 


K 4 = 0 


It should be noted that the product must not differ from 1 by more than 1 x 10" 6 

if the sharp-corner option is to be activated. The angle subtended by the third segment is 


^bo 


3][3 - TT 


tt ~ 1.3242248 rad * 75.872493° 
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The function y b (X,<f>) is equal to the constant y bo for values of X not greater than 
12 AX (in other words, l = k 3 = 13). For values of X greater than 12 AX, the equa- 
tion for y b (X,<3>) is 

y b (X,$) = y bQ - (X - 12 AX)tan(A0 - = y b (X) 

Geometry Input for Class II Bodies 

For class II bodies it is recommended that a one-segment reference-surface gen- 
erator be used so that curvature discontinuities do not occur on the reference surface. 
However, it should be noted that this is a recommendation and not a requirement. The 
expression for y b (X,<h) must be derived for each body and reference surface. It should 
be noted that y b is defined as the distance from the reference surface to the body sur- 
face along normals to the reference surface and is positive when the reference surface 
lies outside the body as it does in figures 2, 3, 4, and 5. 

A prolate spheroid is an example of a class II body. This shape is generated by 
rotating an ellipse about its major axis. The body which is considered here has semi- 
major and semiminor axes of 1 and 2/3, respectively, and is alined so that the major axis 
is normal to the free-stream direction. This body is shown in figure 6. 


x 3 



Figure 6. - Prolate spheroid with axis ratio of 2/3. 
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The coordinates xj, xg, and Xg are related to the angles 0 and $ and the distance 
from the origin o> by the equations 

Xj = o> cos 0 

xg = sin 0 sin <f> 

X 3 = oj sin 0 cos $ 

The equation for the prolate spheroid is 



It follows that the distance from the origin to a point on the surface is - — 

\j 9-5 sin^ 0 sin^i> 

The reference surface is chosen to be a sphere with the curvature K-^ = 1. Therefore, 
the distance from the reference surface to the body along lines normal to the refer- 
ence surface is 


y b (x,*) = 1 - j== . 

y 9 - 5 sin ^0 sin^<f> 

Initial Values 

The shape of the initial shock wave which is used in this computer program is a 
conic section of revolution with the axis of symmetry alined with the free-stream direc- 
tion. The plane of symmetry of a typical flow field is shown in figure 7. 



Figure 7. - Coordinate system for initial shock shape. 
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The quantity a is the angle of attack. The x,r coordinate system has already been 
introduced. The coordinate system is oriented with respect to the initial shock as 
shown in the figure. The initial shock shape is specified by the equation 

R 2 = -2R sX - B s y 2 (6) 

where R s is the nose radius of curvature and B s is the bluntness. The bluntness val- 
ues are -1 for a hyperbola, 0 for a parabola, 1/2 for a prolate ellipse, 1 for a circle, and 
2 for an oblate ellipse. It should be noted that the major axis of a prolate ellipse and the 
minor axis of an oblate ellipse are alined with the y-axis. The initial shock shape for a 
given angle of attack a is specified by the parameters x 0 , y Q , R s> and B s . These 
quantities must be supplied by the user. An approximate location for the shock wave can 
be determined from experimental or theoretical results for similar bodies and flight con- 
ditions or from empirical formulas such as those of Kaattari (refs. 3 and 4). 

It is assumed that the shock is not moving initially. Therefore, the flow properties 
at the initial shock wave can be determined from the assumed shape and the Rankine- 
Hugoniot equations. The initial surface-pressure distribution is assumed to be Newtonian. 
There are several exceptions. At sharp shoulders where the flow would otherwise be sub- 
sonic the pressure is adjusted to a sonic value. If stagnation pressure is not achieved 
anywhere on the surface, the pressures are scaled up so that the largest pressure is the 
stagnation pressure. This situation arises for a flat-face cylinder at angle of attack. In 
shielded or wind-shadow regions the surface pressure is assumed to be equal to that at 
the last unshielded point in the same plane of constant 4?. There is a provision for spec- 
ifying a minimum surface pressure p m ^ n for the initial solution. The density at the 
surface is determined with the normal-shock entropy, and the surface velocity components 
are determined from the total enthalpy and the Newtonian streamlines. The flow proper- 
ties at points between the shock and surface are determined by interpolation. 

It is the opinion of the authors that the most important part of establishing the ini- 
tial solution is positioning the shock wave as accurately as possible. In many cases, a 
considerable amount of computer time is required to move the shock wave a long distance. 
It is not of particular importance to determine the flow properties at the grid points with 
great accuracy. If the shock wave is accurately positioned, it will require about the same 
amount of computer time to obtain a converged solution from very good initial data for the 
flow properties as from fair data. As an illustration of this feature, consider a converged 
solution. Now, let this solution be altered at one grid point in the subsonic region. If the 
time-dependent calculation is continued, this disturbance will travel outward from the 
grid point, in the form of acoustic waves as shown in figure 8. If d is a characteristic 
length for the subsonic portion of the flow field and is the speed of sound at the 
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Origin of — 

acoustic 

disturbance 



Figure 8. - Acoustic waves propagating 
from point of disturbance. 

stagnation point, the time required to reestablish steady flow is of the order of 1 to 
10 times the characteristic time d/af. 

Thermodynamic Models 

The flow of perfect gases and equilibrium air can be calculated with this computer 
program. The task of the thermodynamic routine, which is contained in a subroutine of 
the program, is to determine the speed of sound a and y, the ratio of static enthalpy 
per unit mass h to internal energy per unit mass e, for given values of e and the 
density p. It should be noted that the heat of formation of the gas must be chosen so 
that h = e = 0 for T = 0° K. 

Perfect gas .- For a perfect gas, y is a constant which is supplied by the user, 
and the speed of sound is simply 

a = \j9iy - l)e 

The program is formulated so that y can be given one value in the free stream and 
another value in the shock layer between the shock wave and the body surface. This pro- 
vision permits the approximation of real-gas phenomena with a perfect-gas model. It 
should be noted that if the free stream is cold, the free-stream value of y is the same 
as the free-stream ratio of specific heats y m . A value of y which is suggested for use 
in the shock layer is the normal-shock value y s - This quantity is related to the normal- 
shock density ratio P s /Poo, the free-stream Mach number M m , and by the equation 
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This equation is obtained by solving equations (23) in reference 5. 

Equilibrium air .- The ratio y = h/e and the speed of sound for equilibrium air are 
determined with a curve fit which is discussed in references 2 and 5. This curve fit is 
valid for the following ranges of the density p and internal energy per unit mass e: 

10" 4 5 p/p o ^ 10 and e/ltT 0 = 1500 

where the reference quantities have the values 

p Q = 1.292 kg/m3 (2.507 x 10" 3 slugs/ft 2 ) 

RT 0 = 78.40 kJ/kg (8.439 x 10 5 ft 2 /sec2) 


The user must supply values of the free-stream density in kilograms per cubic meter and 
the free-stream velocity in meters per second. The user must also supply the value 
(1.405) for the free-stream ratio of specific heats y^ and an approximate value for y 
in the shock layer. This value of y, which is used to calculate the initial solution, can 
be obtained from equation (7). The density ratio Pg/p^, in equation (7) can be obtained 
from tables of the equilibrium normal-shock properties of air for various velocities and 
altitudes, such as those in references 6 and 7. 

It should be noted that the gas -model subroutine is written so that other equilibrium - 
gas models can be added to the program. 
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Figure 9.- Function which gives influence of damping coefficient on time-step coefficient. 


Computational Parameters 

The stability of the numerical calculation is governed by the time-step coefficient 
C T and the damping coefficient e, which are specified by the user. The equation for 
the time step is 


At = 


Crg 

V + a 


min(X AX, 6 AY,r A$) 


(8) 


where V is the magnitude of the velocity, a is the speed of sound, g is a function of 
the independent variables and the shock-layer thickness which accounts for the nonorthog- 
onal nature of the coordinate system, and A is the scale factor for the X-coordinate. 

The function g, which is determined within the program, is obtained in appendix B of ref- 
erence 2. The time step is calculated in subroutine DELTAT for each cycle. 

Let the parameter j have values of 0 for plane or axisymmetric flow and 1 for 
three-dimensional flow. The damping coefficient must satisfy the inequalities 
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( 9 ) 


0 < e < e 


max " 


1 

8(j + 2) 


The time-step coefficient C T must satisfy the inequality 



The function f(e/e max ), which is derived in reference 2, is plotted in figure 9. As a gen- 
eral rule, the coefficient C T should be as large as possible. It has been found that a 
practical value of the damping coefficient is e/e max = 0.75. 

There are two options for calculating the flow at points on the body surface where 
the curvature changes discontinuously and where the flow is subsonic. These options 
involve the manner in which the scale factor in the governing equations is treated. In 
general, the first option (I = 1) should be used. However, if the flow-property distribu- 
tions have severe discontinuities in slope at surface points where the flow is subsonic and 
where the curvature changes discontinuously, the second option (I = 0) must be chosen. 

PROGRAM INPUT 


The computer program has 16 fixed-point and 24 floating-point input parameters. 
The symbolic program name, the symbol used in the discussion section, and a brief 
description for each parameter are as follows: 

Program name Symbol Description 


IYREF1 

LF, MF, NF Zj,mj,nf 

DX AX 

KURVE1, KURVE2, k^k^kg,!^ 

KURVE3, KURVE4 


CURVE 1 , CURVE 2, 
CURVE 3, CURVE4, 
CURVE 5 


Geometry index: 0, two dimensions; 
1, three dimensions 

Largest values of grid indices for 
X-, Y-, and ^-coordinates (see 
fig. 4) 

Mesh spacing for X-coordinate 

Values of X index where reference- 
surface generator changes curva- 
ture discontinuously (see inequali- 
ties (2)) 

Curvature of segments of reference- 
surface generator (see eq. (3)) 
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Program name 


Symbol 


Description 


YREF0 


YREF1, YREF2, 
YREF3, YREF4 

IYREF2, IYREF3, 
IYREF4 

RS 

BS 

X0 


CX0 


IGAS 

DENS0, VEL0 


GAMMA 

GAMM 


PINF 

ALPHAD 

KP0 


^bo 


R s 

B s 

x 0 


*o 




Poo 


a 


Perpendicular distance from reference 
surface to body surface (see figs. 2 
and 4) 

Unused floating-point parameters 

Unused fixed-point parameters 

Nose radius of initial shock-wave shape 

Bluntness of initial shock-wave shape 

Distance along axis of reference sur- 
face from body surface to intersec- 
tion with axis of initial shock-wave 
shape (see fig. 7) 

Distance along axis of initial shock- 
wave shape from shock surface to 
intersection with axis of reference 
surface (see fig. 7) 

Gas-model index: 0, perfect gas; 1, 
equilibrium air 

Free-stream density, and magnitude of 
velocity in kg/m3 and m/sec, respec- 
tively, which are used in equilibrium- 
air thermodynamic model 

Free-stream ratio of specific heats 

Ratio of static enthalpy h to internal 
energy e in region between shock 
and body (see section entitled 
"Thermodynamic Models") 

Nondimensional free-stream pressure 
(see eq. (1)) 

Angle of attack in deg 

First time cycle at which output is 
printed 
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Program name 

Symbol 

Description 

KPP 


Number of time steps between successive 
printing cycles 

KF 


Index for final time step 

IYREF0 

I 

Index governing computation of body points 
where curvature is discontinuous (see 
section entitled "Computational 
Parameters”) 

CT 

c T 

Time-step coefficient (see inequality (10)) 

EPSI 

e 

Damping coefficient (see inequality (9)) 

PMIN 

^min 

Minimum nondimensional surface pressure 
for initial solution 


The program is set up to use reference-surface generators with a maximum of five 
segments of constant curvature and four points of curvature discontinuity. If it is desired 
to use fewer than the maximum number of segments, the unused indices KURVE4, 
KURVE3, and so forth should be set equal to some number greater than LF. Some value 
should be given to all unused parameters such as the curvature CURVES, CURVE4, and 
so forth of the unused segments of the reference-surface generator and the quantities 
VEL0 and DENS0 when a perfect gas is being treated. 

The distance y b (X,4>) from the reference surface to the body surface along lines 
normal to the reference surface is specified in subroutine START. The equation for y b 
which is used in this program is 

y b (X,$)=y b0 (11) 

Thus, class I bodies which are located a constant distance from the reference surface 
along lines normal to the reference surface can be treated without any changes to the 
computer program. The symbolic program statement for equation (11) is 
YREF(L,N) = YREF0. This statement immediately precedes statement number 4 in 
subroutine START. 

In order to treat class n bodies and class I bodies which are not located a constant 
distance from the reference surface, it is necessary to replace the statement for 
YREF(L,N) in the program. It should be noted that the angles 6(X) and $ are known 
at this point in the program and have the program names THETA(L) and PHI. Several 
unused fixed-point and floating-point parameters have been provided for use as the need 
arises. It is suggested that these parameters can be used in the new program statement 
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for YREF(L,N), It is emphasized that the parameter YREF0 should not be used in any 
statement other than the one discussed above. 

The program treats class I bodies with corners. The curvature of the arc of the 
reference-surface generator subtending the corner is l/y^ 0 (1./YREF0). In order to 
actuate the sharp-corner option of the program, the product of the value of the curvature 
of the arc subtending the corner and the value of YREF0 can differ from 1 by no more 
than 1 x 10" 6 . 


PROGRAM OUTPUT 

The first quantities to be output are the geometry parameters for the reference and 


body surfaces. These are 



Program name 

Notation 

Description 

THETA(L) 

0(X) 

Angle between normal to reference surface 
and free-stream direction (see fig. 3) 

CURV2(L) 

K 

Curvature of reference surface 

XREF(L), RREF(L) 


Coordinates x and r of points on reference 
surface 

YREF(L,N) 

y b (x,$) 

Distance from reference surface to body 
surface along normals to reference 
surface 

PYBPX(L,N) 

3y b /0x 

Partial derivative 

PYBPPH(L,N) 

9 y b / 8 $ 

Partial derivative 

The complete solution at each grid point is output for the initial cycle (K = 0) and the 
regular printing cycles (K = KP0, KP0 + KPP, KPO + 2KPP, . . .). The output for each 
grid point is composed of four fixed-point numbers and 12 floating-point numbers. Four 
additional floating-point numbers are output for points on the shock wave. For each 
printing cycle the nondimensional time which has elapsed since the calculation started is 
printed. The unit time is lj V^. The program name for this quantity is TIME. 

Program name 

Notation 

Description 

K 


Time index 

L 

l 

Index for X- coordinate (l = 1 at X = 0 


and L = LF at X = X max ) 
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Program name 

Notation 

M 

m 

N 

n 

UN, VN, WN 

u,v,w 

RH0 

P 

P(L,M,N) 

P 

H(L,M,N) 


PHI 

$ 

RAD, XAD 

r,x 


GAMM 

V 

ENTR0P 

P 

AMACH 

M 

DELTA(1,L,N) 

5(X,$,t) 

PDPT(1,L,N) 

86/0T 

PDPX(1,L,N) 

06/0X 

PDPPH(1,L,N) 

96/8$ 


Description 

Index for Y-coordinate (M = 1 at body 
and M = MF at shock) 

Index for ^-coordinate (N = 1 at $ = 0° 
and N = NF at $ = 180°) 

X-, Y-, and ^-components of velocity 
made nondimensional with 

Density made nondimensional with 

Pressure made nondimensional with 



Total enthalpy made nondimensional with 



Azimuthal angle 

Perpendicular distance from axis of ref- 
erence surface and distance along axis 
from body surface in upstream direc- 
tion (see fig. 7) 

Ratio of static enthalpy h to internal 
energy e 

Function of entropy for perfect gas 

Mach number 

Distance between shock and body along 
normal to reference surface 

Component of shock velocity in direction 
normal to reference surface 

Partial derivative 

Partial derivative 


PROGRAM DETAILS 
Arrays Used in Program 

A number of arrays are used in the program. Many of these arrays are functions 
of L,M, and N, the indices for the X-, Y-, and $- coordinates. The arrays which depend 
on L,M, and N are 
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A1(L,M,N), A2(L,M,N), A3(L,M,N), A4(L,M,N), A5(L,M,N), H(L,M,N), P(L,M,N), 
R(L,M,N), U(L,M,N), V(L,M,N), W(L,M,N), C1A(L,I1,N), C2A(L,I1,N), C3A(L,I1,N), 
C4A(L,I1,N), C5A(L,I1,N), C1A0(I2,I1,N), C2A0(I2,I1,N), C3A0(I2,I1,N), C4A0(I2,I1,N), 
C5A0(I2,I1,N), DELTA(I3,L,N), PDPT(I3,L,N), PDPX(I3,L,N), PDPPH(I3,L,N), D1A(L,I1), 
D2A(L,I1), D3A(L,I1), D4A(L,I1), D5A(L,I1), YREF(L,N), PYBPX(L,N), PYBPPH(L,N), 
THETA(L) , CURV2(L), RREF(L), XREF(L) 

The arrays are dimensioned for maximum values of L,M, and N of 17, 6, and 9, 
respectively. These maximum values can be changed by the user if necessary. The 
maximum values of the indices II, 12, and 13 are 3, 4, and 2, respectively, and should not 
be changed. 


Error Statements 

The program contains five error statements which are located in three of the sub- 
routines. The triggering of any of these statements results in the termination of the 
calculation. 

Subroutine PELT AT .- The error statement in this subroutine is triggered when the 
nondimensional quantity V + a in equation (8) is less than 1 x 10"^. When this occurs, 
the calculation is clearly unstable. The quantities which are printed in the error state- 
ment are the time index, the position indices of the point where the instability occurs, and 

Y6 - y b , X, 8, r, AX, AY, A<f>, X AX, r A«&, min(X AX, 8 AY,r A4>) 

The error statement is followed by the complete output at the previous time step. The 
suggested remedy for avoiding the instability is to decrease the coefficient C T . If the 
instability occurs repeatedly for smaller and smaller values of Cr, a more fundamental 
problem which is peculiar to the particular case may be present. 

Subroutine SHOCK .- There are two error statements in this subroutine, both of 
which indicate a breakdown in the solution of a cubic equation at the shock wave. Let 
CSSHK be the cosine of the angle between the outward normal to the shock and the free- 
stream direction, and let CSBX, CSBY, and CSBP be the direction cosines of the shock 
normal. The breakdown of the solution generally occurs because the shock wave develops 
a kink, the magnitude of CSSHK becomes so small that the discontinuity can no longer be 
interpreted as a shock wave, or the flow at the downstream boundary has become sub- 
sonic. The suggested remedy is to change the location of the initial shock wave. 

BLK(l): The quantities which are printed are the time index K, the position indices 
L,M, and N of the point of instability, and the quantities F4, F5, CSSHK, CSBX, CSBY, 
CSBP, F9, and F10, where F4, F5, F9, and F10 are functions of the flow properties. The 
error statement is followed by the full output for the time cycle at which the error occurs. 


24 



The quantity F9 should always be positive. The error statement is triggered when the 
value of F9 becomes less than 1 x 10"6. 

BLK(2): The quantities which are printed are those listed for the previous error 
statement and the quantity Fll, which should never be greater than 1 in magnitude. The 
error statement is triggered when the magnitude of Fll exceeds 1. 

Subroutine START .- There are two error statements in this subroutine, both of 
which indicate that the shape which was input for the initial shock wave cannot, in fact, be 
interpreted as a shock wave. If one of these errors occurs, the user should change the 
shape of the initial shock wave. 

ER(1): The quantities which are printed are the time index K = 0, the position 
indices L,M, and N of the point where the error occurs, and CSSHK, the cosine of the 
angle between the normal to the initial shock shape and the free-stream direction. It 
should be noted that CSSHK should always be negative. The error statement is triggered 
when CSSHK becomes greater than -1.0 x 10"6. 

ER(2): The quantities which are printed are the time index K = 0, the position 
indices L,M, and N of the point where the error occurs, and the shock density ratio. The 
error statement is triggered when the density ratio becomes less than 1. 

Conservation of Total Enthalpy 

One application of the present method is the generation of an initial data line for a 
method-of-characteristics calculation of the supersonic flow on the afterbody. For such 
a calculation, it is necessary that the total enthalpy be conserved, and even small excur- 
sions can jeopardize the characteristics calculation. In the present program the total 
enthalpy is a computed rather than a constrained quantity. The ability of the method to 
conserve this quantity is considered to be an indication of the accuracy and proximity to 
convergence of the solution. In general, the total enthalpy is conserved to within a per- 
cent or so for a converged solution. However, at points on the body surface downstream 
of a sharp shoulder, the error in the total enthalpy is much larger. This is probably due 
to the fact that, the flow there is highly expanded, and the entropy layer lies between the 
surface (m = 1) and the first line of grid points in the shock layer (m = 2) and hence is not 
resolved. This discrepancy in the results for the total enthalpy is of concern from the 
point of view of establishing an initial data line for a characteristics calculation, but it is 
not so important from a purely physical point of view since the flow at the surface down- 
stream of a sharp corner is viscous dominated and probably separated so that the inviscid 
solution in this region does not represent reality. 

The total enthalpy at grid points on the surface and in the region between the shock 
and surface can be constrained to a constant value with relatively simple changes to sub- 
routine PDE. These changes are given in terms of Fortran IV statements: 
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(1) At the beginning of the subroutine define the total enthalpy as 

HT0T=GAMMA*PINF/(GAMMA- 1.)+. 5 

(2) If the equilibrium -air gas model is being used, place the following statements 
immediately after statement 26 in order to estimate the value of the ratio y = h/e at 
the grid location 

HXX=(HT0T*R(L,M,N)-,5*(U(L,M,N)**2+V(L,M,N)**2+W(L,M,N)**2)/R(L,M,N))/P(L,M,N) 
GAMM=HXX/ (HXX- 1 . ) 

Omit these statements if the perfect-gas model is being used. 

(3) Replace the statement following statement 30 with the statement 

ENERS=(HT0T-.5*(U(L,M,N)**2+V(L,M,N)**2+W(L,M,N)**2)/R(L,M,N)**2)/GAMM 

(4) Replace the statement following statement 31 with the statement 

H(L,M,N)=HT0T 

It is not recommended that the total enthalpy be constrained in this manner unless 
it is absolutely necessary since the effect is to ignore the corrective influence of the 
differential-energy equation and hence retard convergence. No attempt should be made 
to constrain the total enthalpy at the shock since this could interfere with the solution of 
a cubic equation for the shock speed. 

It has been found that in highly expanded regions of flow the pressure can assume 
small negative values. In particular, this has been found to be true at grid points on the 
surface downstream of a sharp corner, starting several mesh spacings from the corner. 
Requiring that the total enthalpy be constant does not prevent the negative pressures. 
Constraining the surface entropy to be constant does not help since the effect is to 
decrease the value of the density and hence cause instabilities rather than increase the 
value of the pressure. It is recommended that characteristic initial data lines for bodies 
with sharp corners be located just downstream of the corner in order to avoid the nega- 
tive pressure values. 


Alternate Expressions for Shock Slopes 

Under certain conditions the converged results of many time-dependent, blunt-body 
methods for the shock-layer thickness (distance between the shock and body) do not vary 
smoothly with distance along the body but contain oscillations. The wavelength of these 
oscillations is usually two mesh spacings. Although no oscillations of this type are 
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exhibited by the results computed by the present method and reported in references 1 
and 2, it is possible that they may occur in other cases. 

In some instances the oscillations are due to the type of finite-difference expres- 
sions which are used to represent the derivatives of the shock-layer thickness with 
respect to the coordinates. Under such circumstances, improved results can sometimes 
be obtained by substituting other finite-difference expressions for these derivatives. 
Unfortunately, the expressions which are less likely to produce oscillations are also gen- 
erally of a lower order of formal accuracy. 

In the present program three-point central-difference formulas are used to evaluate 
the derivatives 98/9X and 96/9$ if the velocity components tangent to the shock and 
in the X- and $ -directions, respectively, are less than the local speed of sound. These 
formulas are of second-order accuracy. Two-point upwind-difference formulas are used 
if the velocity components are greater than or equal to the local sonic speed. Although 
these formulas are only of first-order accuracy, they correctly prevent the influence of 
downstream points on the shock-layer -thickness derivatives in supersonic regions. 

In the event that central-difference formulas permit the development of oscillations 
in the shock-layer-thickness distribution with respect to the X- or $ -coordinate, it is 
suggested that upwind-difference procedures be employed for that coordinate in the sub- 
sonic as well as the supersonic regions to calculate the shock-layer-thickness derivatives. 
These changes in the differencing procedure can be effected by modifying subroutine 
SL0PE. In order to obtain windward difference expressions for the derivative 96/9X 
at all points, remove the second statement after statement 6 in subroutine SL0PE and 
replace the next statement with the statement 

IF(PDPX(KKK,L,N).LT.O.)G0 T0 7 

Statement 8 should also be removed. In order to obtain upwind-difference expres- 
sions for the derivative 96/9$ at all points, remove the third ’IF" statement following 
statement 9 and replace the next statement with the statement 

IF(PDPPH(KKK,L,N).LT.O.)G0 T0 10 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 2, 1971. 
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APPENDIX A 


PROGRAM LISTING 


The computational program is listed in this appendix in the Fortran IV computer 
language. This program consists of the main executive program and 11 subroutines. 


PROGRAM EXECt INPUT=20 1» OUTPUT » TAPE 5= INPUT , TAPE6=0UTPUT ) 

C**** CAVIS-DLD-A2636-APPLIED MECHANICS 
C**** CltK eapnwell 

C***« TIME DEPENDENT METHOD FOR COMPUTING FLOW ABOUT BLUNT BODIES 
0**** TRAVELING AT SUPERSONIC SPEEDS AT ANGLE OF ATTACK 

COMMON IGAS,ISND,J,K,XK,KURVE1,KURVE2,KURVE3.,KURVE4,LF,K,MF»MM,N,N 
IF, LLO, DENSO » VELO , KP , KPP ,KF ,LHI ,LCORN 

COMMON AL PH A , B S »CSALP ?CUR VE1 »CUR VE2 »CURVE3 , CURVE A, CURVE 5 » CXO,OPHI , 
* lDT,DX,DY,ElNF,ENERSfEPSI,GAMM,GAMMA,PINF, PM IN, RHO, R S, SNALP,SND,XO , 
2YREFO »YREFl,YREF2*YREF3,YRt:F4 

COMMON AI(17»6,9) ,A2(17,6,9),A3( 17 ,6, 9 ) , A4 ( 17, 6 , 9 ) , A5 { IT , 6,9) ,H( 17 
1,6,9) ,P(17»6,9J,R(17,6,9) ,Utl 7 ,6 ,9 ) , VI17 ,6 ,9 ) , W( 17, 6, 9 ) , CIA( 17, 3, 9 
2),C2A(17,3,9) ,C3A(17,3,9) ,C4A( 17,3,9),C5A( 17,3,9) »CIAQ(4,3,9I ,C2AO 
' 3 ( A , 3 , 9 ) »C3A0(4,3,9),C4A0(4,3,9),C5AO<4,3,9) 

COMMON D1A(17,3),02A(17,3),D3A(17,3),D4A(17,3) ,D5A( 17,3) ,CURV2(17) 
1,THETA(17),RREF(17),XREF( 17) ,YREFl 17,9 ), PYBPXI 1 7, 9) , PYBPPHI 17,9) ,D 
ZELTAt 2,17,9) ,PCPT( 2,17,9) ,PDPX 12 ,17,9) ,PCPPH 2,17,9 ) 

COMMON AAA (5, 5)»AAB(5)»B1A{3) ? B2A ( 3) ,B3A( 3) »B4A( 3) , B5A{ 3 ) ,RRXF(5I , 
1XRXF ( 5 ) ,CT,TIME,IYREF0,IYREF1,IYBEF2,IYREF3,IYREF4 
COMMON /ERROR/ER (12) ,BLK(12) 

CCMMCN/eLCCKl/KPO 

DATA ( ER (I ) ,1=1,12 )/3HERl, 3FER2 ,3HER3 *3HEP4 »3HER5»3HER6» 3HER7, 3HER8 

I, 3FER9, 4HER10»4HER11»4HER 12/ 

CATA ( BLKt l ) ,I=l»12)/4HBLKl,4HBLK2,4KBLK3,4HBLK4,4HBLK5,4HBLK6,4HB 
1LK7,4HBLK8,4HBLK9,5HBLK10»5H8LK11,5HBLK12/ 
NAMELIST/NAME/IYREF1,LF,MF,NF,DX,KLRVE1,KURVE2,KURVE3,KURVE4, CURVE 

II, CURVE 2 ,CURVE3,CURVE4,CURVE5 ,IYREF2,IYREF3,IYREF4, YREFO, YREF1 , YRE 
2F2,YREF3,YREF4,RS,BS f >D,CX0, 1 GAS, DENSO, VELC, GAMMA, GAMM,P INF, ALPHAD 
3,KPO,KPP,KF,IYREFO,CT,EPSI ,PMIN 

80 C FORMAT { //* DA VI5-DLD-A2636-APPLI ED MECHAN IC S-9-22-69* ) 

801 FORMAT ( * DICK BARNWELL *) 

802 FORMAT ( * TIME DEPENDENT METHOD FOR COMPUTING FLOW ABOUT BLUNT BODI 
1ES TRAVELING AT SUPERSONIC SPEEDS AT ANGLE OF ATTACK*//) 

803 FORMAT! 1H1//) 

804 FORMAT t 2A10) 

805 FORMAT ( * A2636 COMPLETED AT CENTER AT *,2A10) 

DIMENSION I DATE { 2 ) 

CALL CAYT1MC ICATE) 

PRINT 8C4, IDATEt PRINT 803$ PRINT 804, IDATE 
PRINT 000$ PRINT 801$ PRINT 802 
1 FEAC(5,NAME) 

IFIEOF, 5)9999, 9990 


9990 CONTINUE 
KP=KPC 

AL PH A= ALPHAD/ 57,25 577951 
WRITE (6, NAME) 

CALL START 
CALL RESULT 

2 KK= 1 
K=K+ 1 

IF(K« GT »KF)GQ TC 1 
CALL DELTAT 
CALL A 

3 M=MF 

CO 4 J= 1,3 
MM=M+J-3 
CALL C 

4 CONTINUE 
CALL SHOCK 
MM7HI=MF-1 

CC 7 MM7=1,MMZHI 
M=MF“MM Z 

1F(M*EC*1.0R.MMZ.EQ.1)G0 TO 6 
DO 5 112*1*2 
J=4-I IZ 
CALL TRANSC 

5 CONTINUE 
J = 1 

MM=M-1 
CALL C 

6 CALL PDE 

7 CONTINUE 
IF(KK.EQ.2)G0 TG 8 
KK = 2 

. GO TO 3 

8 IF(K*NE.KF)G0 TO 2 
CALL RESULT 
KP=KP+KPP 

GO TO 2 

9999 CALL, CAYT I M < ICATE )$ PRINT 605,IDATE 
STCP C1C1 
END 
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SUEROLTINE A 

CCMMON IGAS,ISN0»J,K # KK f KURVEl,KURVE2,KURVE3»KURVE4,LF,M,MF,MM»N,N 
lF»LLO»DENSO?VELO»KP»KPP,KF»LHI»LCORN 
CCMMON ALPHA, BS»CSALP, CUR VE1, CUR VE2,CURVE3,CURVE4,CURVE5,CXO,DPHI» 
1CT,DX,CY,EINF,ENERS,EPSI,GAMM,GAMMA,PINF,FMN,RHC,RS,SNALP,SND,X0, 
2YREFO,YREFl,YREF2,YREF3,YREF4 

COMMON *1(17,6,9} ,A2( 17,6,9) »A3117,6,9J» A4( 17,6,91, A5( 17,6*9) ,H(17 
1,6,9) ,P(17,6, 9), R( 17,6,9) ,U( 17, 6,9), V( 17, 6, 9) ,H( 17,6,9) ,C1A( 17,3,9 
2) ,C2A( 17,3,9) ,C3A( 17,3,91 ,C4A ( 17 , 3 ,9 > , C5A ( 17 , 3 , 9 ) , C 1AO ( 4, 3, 9) , C2AO 
3 (4,3, 9) ,C3AO( 4,3,9) ,C4A0( 4,3,9) ,C5AO( 4,3, 9 ) 

CCMMCN C1A(17,3),C2A(17,3 ) ,D3A( 17 , 3) ,04A ( 17 ,3 ) ,D5A< 17,3) ,CURV2(171 
l , THETA ( 17) , RREF( 17) ,XREF( 17) ,YREF( 17,9), PYBPXt 17,9) ,PYBPPH( 17,9),D 
2 ELT A ( 2,17,9), PDPT( 2,17,9) ,POPX ( 2 , 17 , 9 ) ,PDPPH( 2 ,17,9 ) 

CCMMON AAA (5, 5),AAE(5) ,B1 A(3) ,B2A(3) , 83A ( 3 ) ,B4A ( 3 ) , 85A(3) ,RRXF(5) , 
1XRXF (5) ,CT,TIME,IYREFQ» IYREFI , JYREF2, 1YREF3, I YREF4 
CCMMCN/ ERPOR/EP (12 ) ,BLK( 12) 

MFM1=MF-1 
ISND= 1 

CC 11 M= 1 »MFM 1 
l LC»1 

CO 10 N=1 ,NF 
I K Z-0 

CO 9 l«LLC,LF 

IFtM.EG.l.AND.L.GT . LCCRN, AND.LoLT .LHI )G0 TC 9 
IF(M.EQ.l)ClA(LfM,NI = W 
CO 7 11 = 1 ,5 
LL*L+I 1-3 

IF(LL .LE.LF )G0 TO 2 
DO 1 JJ=1*5 
IIM1=II-1 
1 1 M2= 1 1-2 

AAA( JJ,II)*2* *AAA ( J J , 1 1 Ml )-A AA ( J J , 1 1 M2 ) 

1 CONTINUE 
CC TC 7 

2 IF(L.EQ.LL0.CR.II.E0.5)G0 TO 4 
IF(M®EQ«l®ANC®L®EQ®LHI )G0 TO 4 
CO 3 J J = 1 ,5 

I IP 1= 1 1 + 1 

AAA(JJ,II)=AAA(JJ» IIP1 ) 

3 CONTINUE 
CC TO 7 

4 IF(LLpLT„l)GC TO 5 


LLL=LL 

nn*n 

Fl=l. 

C-C TO 6 

5 LLL=2“LL 
NN=NF+1-N 
Fl*-1. 

6 AAA(1, II)=R(LLL,M,NN) 

AAA(2, II)=F1*U(LLL,M,NN> 

AAA<3,U)*V(LLL,M,NN) 

MA(4,I I)=W(LLL,M,NN) 

AAA<5,II)=R(LLL,M,NN)*H(LLL,M,NN>-P(LLL,M,NN) 

7 CONTINUE 

IF( M® NE® 1 ) GO TC 70 
LPl=L+l 

IF(KUBVE1.NE.LP1.AND.KURVE2.NE.LP1.AND.KURVE3.NE.LP1.AND.KURVE4.NE 
l.LPl.ANO®IKZ.EQ.O)GO TO 70 
IF ( I KZ . EO» 1 ) GO TO 74 
IF C IKZ® EQ® 2 ) GO TO 75 
IF( IKZ«EQ®3)G0 TO 72 
L=L + 1 

PHC=R(L,M,N> 

ENERS = H(L,M,N)-(P(L,M,N)+.5*(U(L,M,N)**2+V(L,M,N)**2+W(L ,M,N)**2> / 
1R(L,M,N))/R(L ,M,N) 

CALL GAS 

CFECK9=U(L,M,N)/(R(L,M, N) *SN01 
L =L~ 1 

IF(LP1.NE.LC0RN)G0 TO 69 
I KZ=1 
GO TO 73 

69 IF ( CHECK9*LT • 1 •*AND®CHECK9®GT • -1 ® * AND® IYPEFC® NE® C )GC TO 70 
CHECK 8-CUR V2 ( L )-CUR V2 ( LP1 ) 

IF(CHECK9*GE.1..AND*CHECK8.GT.C. I IKZ= 1 

IF (CHECK9.LE.-1.® AND.CFECK8.LT ®0. ) IKZ=2 
IF I IKZ.EO.O) IKZ=1 
GC TO 73 

70 CO 80 J J= 1 , 5 

AAB(JJ)=AAA< JJ,3)-EPSI*(6.*AAA(JJ ,3)-4.*<AAA( JJ,2)+AAA( JJ,4) )+AAA( 
1JJ,1)+AAA( JJ,5) ) 

33 CONTINUE 
C-C TO 8 

72 IKZ-O 

73 CC 83 J J= 1 » 5 
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AAB< J J> =AAA< JJ,3) ♦EPS I* ( A A A C J J ,2 ) -2. * A A A ( J J # 3 ) + A AA ( JJ ,4 ) ) 

83 CCNTINUE 
CG 70 e 

74 I J 1 = 2 
IJ2 = 1 

IF (|_ * NE*LCCP N )GG TO 76 
IF(CHECK9,LT«1« )C1A(L ,N 1=0. 

GO TO 76 

75 I J 1=4 
I J2=5 

76 CC 86 JJ=1,5 

IFU_«EC.LCDRNoAND.CHECK9.LT.»l.)G0 TO 860 

AAB( JJ)=AAA(JJ,3)-EPSI*(AAA< J J ,3 ) -2.* AAAI J J, I Jl) *AAA( JJ, IJ2M 
GO TO 86 

860 AA6I JJ)=AAA( JJ,3> 

86 CCNTINUE 
IK Z=3 

8 Al(l*M,N)=AAB(l) 

A2(L,M»N)=AAB(2) 

A3 ( L * R , M = AAB (3) 

A4U_,R»N)=AAB(4) 

A5IL,M, N)=AABI5I 

9 CONTINUE 
LL0=2 

10 CCNTINUE 
2 1 CCNTINUE 

IF(NF.LE.2)GC TO 220 
CO 22 R=1,MFW1 
NUP= 1 

CC 21 L = 1 »LF 

IF (M* EC. 1. AND.L.GT .LCQRN. AND. !_• LT« LHI ) GO TC 21 
CO 20 N=1,NUP 
CO 18 11=1,5 
NN=N+ T 1—3 

IF(NN.GT.NF.AN0.II.EC.5)GC TC 12 
IFIN.EQ.l.OR® II.EQ.51G0 TO 15 
! IP1= I I + l 
F2=i. 

GO TO 13 

12 F2=-l a 
TF(N.LT.NF) I I P 1=3 
IFIN.EC.NF) IIP 1=1 

13 CO- 14 J J=1 , 5 


IF (JJ«NE*4)F3=1« 

IF(JJ«,EQ.4)F3=F2 
AAA(JJ,1I ) = F3*AAA{ JJ, II PI } 

14 CCNTINUE 
GO TO 18 

15 IFCNN.LT.1IGC TO 16 

F4= la 

NNN=NN 
GO TO 17 

16 F4=- 1 » 

NNN=2-NN 

17 A A A ( 1 , 1 I ) =R ( l V M,NNN) 

AAA(2, I I»=U(L,M,NNN) 

AAA(3,II)*V(L »H»NNN) 

AAA(4, m=F4*WtL,M,NNN) 

A A A ( 5 , 1 1 1 =R( L ,NNNJ *H(Lt R ,NNN l-P (L , M, NNN ) 

18 CCNTINUE 

CC 19 J J=1 » 5 

AAR( JJ)=-EPSI*(6.*AAA( JJ, 3)-4,*< AAAI JJ,2)-*AAA< JJ,4> l+A AA ( J J , 1 ) +A A A 
1 ( JJ , 5 I ) 

19 CCNTINUE 

F5= 1 . 

IF(M,EQ.1)F5=C1A(L,M,N) 

A1(L,N,N)=A1(L,N,N) +F5*AA6 ( l ) 

A2IUM, K)=A2(L,M«N)+F5*AAe<2l 
A3(L,M,N»=A3(L,P, N)+AAB(3 ) 

A4(L,N,N) =A4 (L,M,N)+AAB(4 I 
A5(L,M,M = A5(L,M,N)4F5*AAE(5I 

20 CCNTINUE 
NUP=NF 

21 CCNTINUE 

22 CCNTINUE 

220 IF I MF .LE.3 I GO TO 39 
NUP= 1 

CO 38 L = 1 ,LF 
RLC= 1 

IFIL.GT.LCORN.ANO.L.LT.LH I )MLC = 2 
L N 1= L - 1 

IFlt.EQ.DLMl =1 
CO 37 N= 1 »NUP 
CO 36 M=ML0,RFR1 

1 1 L G— 1 
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IF(M.EC.iniL0*3 
IF(M.EC.2)I ILQ=2 
IIUP=5 

IF(M.EQ«MFM1)T IUP=4 

00 25 1 1= I ILC * I IUP 

IF(M.EG.ML0.CR.II.EQ*5)G0 TO 24 

1 1 PI* 1 1 +1 

00 23 JJ=1,5 

AAA( J J* II)=AAA(JJt I I P 1 I 

23 CONTINUE 
CO TO 25 

24 NM=M+II-3 
APM=MM 

AL AHB=1 *-«-CURV2( LM 1 ) * ( (AMM-l. )*DY*DELTA ( I ,L *N1 -YRFF ( L , N ) ) 

AAAtl ,II)=R(L»f'M,N)*ALAMB 
AAA (2*11) =U(L *MMtN)*ALAMB 
AAAI 3, II)*V(L .MM, N)*ALAMB 
AAA (4* I I)=W(L ,MM,N)*ALAMB 

AAA(5* I I>=(R<L,MM,N)*Hf l,MM*N)-P(L ,MM,N) )* ALAMB 

25 CONTINUE 

ALAHB=1.+CURV2(LM1I*( (AM-1. )*DY*C£LTA(l f L,N)-YREF(L ,N) ) 
IFCM.EQ.DGO TO 27 
IF ( M. EQ .2 ) GO TC 31 
IF (M. EQ.MFM1 ) GC TO 32 
00 26 J J=1 * 5 

AAB(JJ)=-EP$I%(6.*AAA(JJ, 3)-4.*( AA A ( J J *2 ) 4 A A A ( J J , 4 ) I +A AA ( JJ , 1 ) 4AAA 
1 ( J J * 5 ) ) /ALAMB 

26 CONTINUE 
GO TO 35 

27 12=4 
13=5 

29 CO 30 J J= 1 * 5 

AAB( JJ)=-EPSI*(AAA< J J,3 )-2.*AAA ( JJ,I2)+AAA (JJ,I3) > 

30 CONTINUE 
GC TC 35 

31 11=2 
12=3 
13=4 
14=5 

GO TO 33 

32 11=4 
12=3 
13=2 


14=1 

33 DO ’ 34 J J=1 *5 

AAB( JJ)=-EPSI*(AAA(JJ, I4)-AAA( JJ , 1 1 )-3.* < A AA( J J , 1 3 ) -AAA ( J J, 1 2 ) > >/A 
1LA*B 

34 CONTINUE 

35 F 5= 1 . 

IF ( M« EG* 1 ) F5=C 1A{ L*P *N) 

Al(L*M*M=Al(LtM,N)4F5*AAB(l) 

A2(L,M,N)=A2(L,M»N14F5*AAE(2> 

A3(L,M,N)=A3( L*M,N)4AA8(3 ) 

A4 ( L * M * N ) =A4 (L*M*N)+AAB(4) 

A5(L,M,M=A5( L,M,N)+F5*AAB(5) 

36 CONTINUE 

37 CCNTINUE 
NUP=NF 

38 CCNTINUE 

39 RETURN 
END 



SUBROUTINE C 

COMMON IGAS, lSND,J»K,KK # KURVEl,KURVE2,KUPVE3fKURVE4,LF,M,MF,MM,N,N 
1F,110,OEN$0,VELO,KP,KPP,KF,IHI ,LCCRN 
COMMON ALPHA, BS,CSALP,CUR VE1 ,CL'R VE2»CURVE3 ,CUR VE4,CURVE5 ,CXO ,DPHI , 
lCT,DX,DY f EINF,ENERS,EPSI,GAMM,GAMMA,PINF, PMIN, RHO, RS, SNALP,SNO, XO , 
2YREFO , YREF1 »YREF2 »YREF3»YREF4 

COMMON AI(17,6»9)*A2(17,6,9),A3< 17,6,9 ) , A4 ( 17,6,9) ,A5( 17,6,9) ,H( 17 
1,6,9) ,P( 17,6, 91, R( 17,6,9) ,U( 17,6 ,9) ,V( 17,6,9) ,H(17,6,9) ,CIA(17 ,3,9 
2 ),C2A( 17,3,9) ,C3A( 17,3,9) ,C4A{ 17, 3 ,9) ,C5A C 17, 3,91 ,C 1AO ( 4 , 3 , 9) , C2AO 
3(4,3, 9) ,C3A0 (4,3,9), C4AQ ( 4*3»9),C5AQ(4*3,9) 

COMMON 01A(17,3),C2A<17,3»,D3A(17,3> ,D4Atl7,3) ,D5A(17,3> ,CURV2(17) 
1,THETA( 17),RREF(17) ,XREF( 17) , YREF ( 17, 9 ) , PYBPX( 17,9) ,PYBPPH( 17,9) ,0 
2FLTA(2,17,9) , PCPTC 2 , 17,9) , PDPX (2 , 17 ,9) ,PCPPH(2, 17,9 ) 

COMMON AAA l 5 , 5) , AAB ( 5 ) » B1 A ( 3 ) , B2A ( 3) ,B3A ( 3) ,B4A( 3) , B5A ( 3 ) ,RRXF<5) , 
1XRXF ( 5 ) f CT,TIME,IYREFO,IYREFl, IYREF2, I YR EF3 , I YREF4 
CCMMCN/ERROR/ER(12),BLK(12) 

A MM = MM 
LL0=1 
F 3=0 « 

DO 4 N= 1 , NF 
CO 3 L=LLO,LF 
Y=(AHM-lo )*DY 

F AD*RREF(L) + ( Y+OEITA (KK,L , N) -YPEF (L ,N ) )*SIN (THETA (L ) ) 
F4=PDPX(KK,L,N) 

L J=L 

1 ALAMB=1.+CURV2(LJ ) *( Y*OELTA( KK ,L , N ) -YREF ( L ,N ) ) 

IFU.EC.DGO TC 10 

F3=(PYBPPH(L,N)+Y*PDPPH(KK,l,N))*ALAMB/RAC 
10 F 1=Y*P0PT (KK, L ,N ) *ALAMB 
F2=PYBPXU.,N>*Y*F4 

C10=ALAMB*V(L,MM,N)-F1*R<L,MM,N)-F2*U(L,MM,N)-F3*W( L, MM, N ) 
C2C=C1Q*U(L,MM,N)/R(L,MM, N)-F2*P(L,MM,N) 
C3G=C10*V(L,MM,N)/R(L,MM,M+ALAMB*P(L,MM,N) 
C4C=C1C*H(L,MM,N)/R(L,MM,N)—F3^P(L,MM,N) 
C5C=C1C*H(L,MM,N)+F1*P(L,MM,N) 

IF (L • GT «L J ) GO TO 2 
CIAU,J,NI*C1C 
C 2 A ( L >4 »N)=C20 
C3A(L,J,N )=C30 
C4A(L,J,N)=C4G 
C5A(L,J,N) =C50 

IF ( L. NE» KURVE 1 «AND« L»NE»KURVE2« AND*L •NE,KURVE 3® AND* L* NEo KURVE4 )G0 


1 TO 3 
LJ=L-1 

F4=F4*< 1*+CURV2(LJ)*(DELTA<KK,L,N)-YREF(L,N)) )/(l.+CURV2 CL )*( DELTA 
1(KK,L,N)-YREF(L,N) ) ) 

GC TO 1 

2 IF (L • EQ#KURVE 1 I JJ J*1 
IF(L.EC.KGRVE2UJJ=2 
IFIL.EC.KURVE3) JJJ=3 
IF{L.EC„KURVE4)JJJ=4 
C1A0( JJJ,J,N) =C10 
C2AG(JJJ,J»N)=C20 
C3AC(JJJ,J,N)-C30 
C4A0( JJJ,J,N)=C40 
C5A0( JJJ, J,N)*C50 

IF { MM« EQ« 1, ANO.IYREFO.EQ.O)C3A(L,2,N)=V(l,2,N)/RU,2,N) 

3 CONTINUE 
LLC=2 

4 CONTINUE 
RETURN 
END 
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SUBROUTINE 0 

CCMMON IGAS,ISNC,J»K,KK,KURVE1,KURVE2,KURVE3,KURVE4,LF,M,MF,MM,N, N 
IF f ttO, DENSO * VELO, KP , KPP,KF,LHT ,LCCPN 
COMMON ALPHA, BS ,CSAL P , CUR VE1 ,CUR VE2 ,CURVE 3 ,CURV E4 ,C UR VE 5 , CXO ,DPHI , 
ICT,DX,DY,EINF , ENERS ,EPS I , G AMN , GAMP A, PINF, PM IN , RHO, RS , SN ALP, SND , XO , 
2YREF0,YREF1,YREF2,YREF3,YREF4 

COMMON Al< 17,6,9), A2( 17, 6,9), A3 < 17,6,9), A4 <17,6,9 ),A5( 17,6,9) ,H( 17 
1,6,9) ,P(17, 6,9) ,R( 17,6, 9) ,U( 17, 6,9), V< 17, 6,9) ,W( 17, 6,9), C1A< 17,3,9 
2) ,C2A ( 11,3,9) ,C3A( 17,3,9) ,C4A< 17,3,9) ,C5 A ( 17, 3,9) ,C 1AG ( 4, 3, 9) ,C2AO 
3 (4, 3, 9) ,C3AC(4,3,9),C4A0(4,3,9),C5A0(4,3,9) 

CCMMON D1A( 17 ,3),D2A( 17,3) ,D3A( 17,3) ,04A 117,3) ,05A( 17,3) ,CURV2<17) 
1 , THETA { 17) ,RREF(17) ,XREF( 17) , YREF ( 17, 9 ) , PYBPXI 17,9) ,PYBPPH( 17,9) ,D 
2 FLTA( 2,17,9) , PCPT ( 2 , 17 , 9) ,PDPX (2 ,17,9 ) ,PDPFFf 2, 17,9 ) 

COMMON AAA{5,5),AAB(5),B1A{3),B2A( 2) ,B3A(3) , B 4A ( 3) , B5A { 3 ) ,RRXF ( 5 ) , 
1XRXF ( 5 ) ,CT,TIME,IYREFC,IYREF1,IYREF2,IYREF3,IYREF4 
CGMMCN / ERROR/ EP( 12) ,BLK(12) 

AM=M 

CO 10 I J=l,3 

IF{N*EQ*NF*ANO«IJ«EQ»3)GO TO 6 
IF<N.GT*UAND, IJ.LT.3)G0 TC 7 
IF(N»EQaleAND«IJoEQ«3)GO TC 6 
IF<N.EC.1.ANC,IJ.E0.1)GO TC 1 
NN=N+ I J-2 
Gl= 1 o 
GC TO 2 

1 NN=2 
Cl=-1 o 

2 DO 5 L=LLO,LF 

IF(M,EC®1«ANC.L*GT »LCORN* AND® L ®L T • LHI )GO TC 5 
IF ( NF®LE®2) GC TO 40 
IF (L« EQ® 1 ) GO TO 3 
Y= ( AM-1, )*DY 

ALAMB=1*+CURV?(L)*<Y*DELTA<KK,L,NN>-YPEF (L,NN) ) 

L M 1=L- 1 

!F(M*EC.1.AND®L.EG®LC0RN) AL AMB=l .-CURV2I LN1)*YREF(L ,NN) 
PAD=RREF(L)+( Y*DELTA ( KK ,L , NN) -YREF (L , NN) >*S IN (THETA (L ) ) 

C2=AL AMB/RAC 
C-0 TO 4 

3 G2=0® 

4 C 1A ( L , I J ) =G1*G2*W ( L » M » NN) 

C2A(L,IJ) 3 =D1A(L,IJ)*U(L,M,NN)/R(L,M,NN) 

C3A(L,IJ)=D1A<L, IJ)*V(L,M,NN)/R( L,M,NN> 


C4A(L,I J|=G2*(P(L,M,NNH-W(L,M,NN)**2/R<L,M,NN) ) 
F5ACL, I J)=D1A(L, I J)*H(L,M,NN) 

GC TO 5 

40 C1A(L, I J)-0« 

C2A(L ,1 J)=0. 

C3A(L, I J )~0« 

C4AI L , I J )-0« 

DEA(L,IJ) =0* 

5 CONTINUE 
GC TO 10 

6 IJJ=1 

, Gl=-lo 
GO TO 8 

7 IJJ=IJ+1 
Gl*l. 

8 CO 9 L~LLO»LF 

C1A(L, I J)=G1*D1A(L, IJJ) 

C2A ( L , I J ) =G 1*C 2A ( L , I J J ) 

C3A(L, IJ)=G1*C3A(L, IJJ ) 

C4A(L,I J)=D4A(L, IJJ) 

C5A(L,TJ)=G1*D5A(L» IJJ) 

9 CONTINUE 
10 CONTINUE 

RETURN 

END 



SUBROUTINE DELTA! 

COMMON IGAS, ISND, J,K,KK ,KURVE 1 ,KURVE2 ,KUPVE3 ,KURVE4 ,LF , M , MF , MM, N, N 
lF,LLC,DENSO,VELO,KP»KPP»KF,LHI,LCORN 
COMMON ALPHA,BS,CSALP,CURVE1»CURVE2»CURVE3,CURVE4,CURVE5,CX0,0PHI , 
ICT, OX ,CY,EINF»ENERS»EP$I, GAMM, GAMMA, PINF , PM l N , RMC , R S , SNA LP » SND , XO , 
2YREFO » YREF1 , YREF2 , YREF3 » YREF4 

COMMON Al(17, 6,9) ,A2< 17,6,91 ,A3( 17,6,9) ,A4( 17,6,9) ,A5<17,6,9),H( 17 
1,6,9) ,P( 17,6,9),R(V7,6,9) ,U(17,6,9),V(17,6,9) f W(17,6,9),ClA(17,3,9 
2) ,C2A( 17,3,9) ,C3A( 17,3,9) ,C4A( 17,3,9) ,C5 A ( 1 7, 3 ,9 ) ,C 1 AO (4, 3, 9 ) , C2 AO 
3(4,3,91 ,C3AO(4,3,9),C4AO( 4,3,9),C5AO( 4,3,9) 

COMMON C1A( 17 ,3 ) »D2 A ( 17,3 ) *D3A(17 ,3) , D4A ( 17,3) ,D5A( 17,3) ,CURV2(17) 
1 , THETA ( 17I,RRFF(17) ,XREF( 17) , YREF ( 1 7 , 9 ) , PYBPX { 17 ,9) , PYBP PM 17 , 9 ) , 0 
2ELTA(2, 17,9), PD PT( 2,17,9) , PDPX (2 , 1 7 , 9 ) , PDPPHC 2 , 17, 9 ) 

COMMON AAA (5, 5) ,AAB<5) , B1 A ( 3 ) , B2 A ( 3 ) , B3A (3 ) , B4A ( 3 ) , B5A ( 3 ) , RRXF ( 5 ) , 
1XRXF(5) ,CT,TIME,IYREFO, I YREF 1, IYREF2 , I YREF 3 , I YREF4 
CCMMCN/ERRCR/EP(12) , 8LK(l 2 ) 

900 FORMAT { * SUBROUTINE DELTAT * , 1A 6, A I 5/ ( BE 16.8 ) ) 

I SND= 1 
C T = 100« 

DC 7 M= 1 , KF 
A M=M 

Y= ( AM- 1 • ) *DY 
L LC=1 

TO 6 N=1,NF 
CO 5 L=LLO,LF 

FO=Y*DELTA( 1 , L ,N ) -YREF ( L , M 

IFtM.EQ.l.AND.L.GT.LCORN. AND.L.LT.LHI )G0 TO 5 
LP = l 

IFtM.EQ. 1 # AND,L« EG, LCORNI LP=L-1 

ALAMB=i.+CURV2(LP)*F0 

Fl=(Y*POPX(i,L,M+PYePXU,M l/ALAMB 

CS1=ALAMB*DX 

TF( L* EC«1 )G0 TO 3 

RAD=RREF(L»+FC*SIN(THETA(U ) 

F2= ( Y *P DP PH ( 1 ,L ,N ) + PY 8PPH ( L ,N ) )/RAD 
CS2-R AD*DPHI 
CO TO 4 

3 F 2=0 • 

CS2=100. 

4 CC=SQRT(Fl**2*F2**2 ) 

F 3=SQRT (4.+CC**2) 

CT1=1« /SORT ( 1 5*CC*(CC+F3I ) 


CS=DEL7A( l,L,N)*DY 
IF(DS1.LT.DS)DS=DS1 
IF{DS2.LToDS ) DS=DS2 

FNERS-H(L,M,N)-( P( L,M,N)+„5*(U(L,M,N)**2+V ( L , M , N ) ** 2 + W ( L ,M,N>**2 ) f 
1P(L,M,N) )/R(L,M,N) 

PHC=R(L,M,N) 

CALL GAS 

DEN0M=SCRT(U(L,M,N)**2-»-V(L,M f N)**2+W(L,M,N 1**2) /R(L ,M,NI + SND 
IF(DENOM,GTo l.E-12) GO TO 40 

PRINT 900 , BLK (l),K,L,M f N,FC,AL AMB, DELTA! 1 ,L,N) ,RAD , 0 X ,D Y , DPHI ,DS1 , 
1CS2,DS,CC,ENERS,RHC,SND 
K = K- 1 

CALL RESULT 
STOP C5C1 
43 CT2=1 • / DENOM 

CT1=CT*CT1*CT2*DS 

TF(DT1.LT.DT)DT*DT1 

5 CONTINUE 
LLC*= 2 

6 CONTINUE 

7 CONTINUE 

T IME=TIME + DT 

RETURN 

END 
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SUBROUTINE GAS 

COMMON IGAS,ISND, J , K , KK ,K LRVE1 ,KUR VE2 * KUP VE3 , KURVE4 ,LF , M , MF, MM,N, N 
IF ,LLO,DENSO,VELO,KP,KPP,KF,LHl ,LCCFN 
CCKMCN AlPHA,BS,CSALP,CURVEl,CURVE2,CURVE3,CURVE4,CURVE5,CXO,DPHI , 
1CT,DX,DY, EINF,£NERS,EPSf ,GAMM,GAPMA,PINF,FMN,RHC,PS,SNALP,SND,XO, 
2VREFC,YPEF1,YREF2,YREF3,YREF4 

CCMMCN Al<17,6,9) ,A2<17,6,9) ,A3< 17,6,9),A4< 17,6,9), A5(17,6,9) ,H<17 
1, 6, 9), PI 17, 6, 9), R (17,6,9) ,U( 17,6,9) ,V( 17,6,9) ,W( 17,6,9) ,C1A( 17,3,9 
2) ,C2A( 17,3*9) , C3A(17,3,9) ,C4A< 17 , 3 ,9 ) ,C5A ( 17, 3 ,9) ,C 1 AO { 4 , 3, 9 ) , C2AO 
3 (4, 3, 9) ,C3AG( 4,3,9) ,C4A0( 4,3,9 ) ,C5 AO (4, 3,9 ) 

CCMMCN E 1 A ( 17,3),C2A(17,3),03AC17,3 ),D4A (17,3) ,D5A< 17,3) ,CURV2C17) 
1 ,THET A( 17) ,RREFC17) ,XREF{ 17) , YREF < 17, 9), PY8PXI 17,9) ,PYBPPH( 17,9) ,0 
2FLTA(2?17,9), PCPT (2,17,9) , PDPX 1 2*17, 9), PTP PH 12, 17,9) 

CCMMON AAA(5»5)»AAE(5)»B1A(3)-*B2A(3I*B3A(3) ,B4A(3) , B5AI3) ,RRXF(5) , 
1XRXFC5) ,CT,TIME,IYREFQ, l YREF 1,IYREF2,IYRFF3,IYREF4 
CCMMCN/ERR0F/ER(12) »BLK(12) 

IF I IGAS.GT.CIGC TO 2 
IF( ISN0«EQ.0)G0 TO 101 
GAME*0. 

GAMR=0. . 

GO TO ICO 

2 I F C IGAS.GT.l ) GC TO 7 
PHO= ABS (RHO ) 

Y2=AL0G 10 ( RHC*OENSO/ 1 .292 ) 

Z3=ABS ( ENERS ) 

Z2=AL0G10(Z3*VELG**2/78400. ) 

IF(Z2.GT..8C1)G0 TO 3 
G AMM= 1 • 405 

IF ( ISNO.EC. 0 ) GG TO 101 
GAME=0. 

GAMR=0. 

GO TO ICO 

3 1F(Z2.GT.2.3)GC TO 5 
TESTY2=(3.255-2.278*Z2)/< l.-.822*Z2) 

IF(Y2.GT.TESTY2) GO TO 4 

GAS1 = 1 * 637~.0404*Y2 
GA52=.2175~.0332*Y2 
GAS 3= .136 6- .036 6* Y2 
CAS4=.0833-.0248*Y2 
GAS5=.0404-.0332*Z2 
CAS6=.0366-.0248*Z2 
CAS7=EXP(-18.3*Z2+.586*Y2*35.57> 


GAS8=-18.3 
G AS9=. 5C3 
GO TC 6 

4 GAMM=l.5055~.1255*Z2 
IF 1 1 SNO.EQ.O) GG TO 101 
G AME=-. 2892 

GAMR=0* 

GO TC 100 

5 GAS1= 1. 5004-. CC38* Y2 
CAS2=.1342~.CC 84*Y2 
GAS3-.3274+.0C91*Y2 
GAS4=. l342+.001fc*Y2 
GAS5=.0C38-.0C84*Z2 
GAS6=-( .0091-.CC16*Z2) 

GAS7=EXP(-15.t 5*Z 2+»507*Y2+42.91) 

G AS8=-l 5.85 

G AS9=. 5C7 

6 GAS1C=1.Z(1.*GAS7) 
0AMM=GAS1-GA$2*Z2-(GAS3-GAS4*Z2)*GAS10 
IFUSND.EQ.OGO TO 101 

G AS1 1 = (GAS3-C-AS4*Z2)*G AST *GAS 10**2 
0 AME= 2. 304*<-GA$2+GAS4*GASlC+GASll*GAS8) 
GAMR=2.304*(-GAS5 + GAS6*GAS10+GASU*GAS9) 
GO TO 100. 

7 GO TO 101 

100 SNDSQ=ENERS*< (GAMM-l. )*(GAMM+GAME) -»GAMR ) 
SND=SQRT( ABS ( SNDSQ) ) 

101 RETURN 
END 
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SUBROUTINE PDE 

CCMMON IGAS, ISNO, J , K,KK,KURVE 1 ,KURVE2 ,KUPVE3 >KURVE4 ,LF ,M,MF ,MM,N,N 
IF , LLOiDENSO* VELC »KP f KPP»KF *LHI ,LCOFN 

COMMON ALPHA, 9S»CS ALP, CURVE 1 ,CURVE2,CURVE3 »CURVEA ,CURVE5 ,CXO,DPHI , 
lOTrOX ,DY, El NF , EKERS, EPS I * GAMM, CAMM A , P INF, PM IN , RHC,R S , SNALP , SND , XO , 
2>REFO,YREFl,YREF2,YREF3,YREFA 

CCMMCN A1U7,6,9),A2(17,6,9),A3U7,6,9),A4{ 17,6,91 ,A5<17,6,9) ,H(17 
1,6,9) ,P(17,6,9) ,R( 17,6,9) ,U( 17, 6,9) ,V( 17, 6,9) ,W< 17,6,9) ,C1A( 17,3,9 
2 J,C2A< 17,3,9) ,C3A( 17,3,9) ,CAA( 17,3,9) ,C5A ( 17, 3 ,9) ,C 1AO { A ,3,9 ) ,C2AO 
3 ( A ,3 ,9 ) ,C3AO( A,3,9),CAAO{ A , 3 , 9 ) , C 5 AO ( A , 3 , 9 ) 

Ct«MON DIAt 17 ,3) ,D2Al 17,3 ) ,D3A (17 ,3) ,CAA 117,3) ,D5At 17,3) ,CURV2(17) 
1 ,THET A( 17 I , RR EF( 1 7 )‘, XREF( 17) ,YREF( 17,9) ,PY0PX( 17,9) ,PYBPPH( 17,9) ,0 
2ELTA(2,17,9) , PDPTI 2 , 17,9) , POP X (2 , 17 , 9 ) , PCP FH { 2 , 1 7, 9 ) 

COMMON AAA (5, 5 ) , AAB ( 5 ) , B1 A { 3) , B2A ( 3) , B3A ( 3 ) ,BAA ( 3 ) , B5A < 3 ) ,RRXF(5) , 
1 XRXF ( 5 ) ,CT*TIME,IYREFG, IYREF1, IYREF2, I YR EF 3 , I YREFA 
CCMM0N/EPR0R/ER(12) ,0LK(12J 
AM=M 
PH I = Q» 

I SND=0 
XKK=3— KK 
LLC-1 
L 1-2 
M=NF 
Fl*-1. 

CO 33 N«1 ,NF 
CALL 0 
IKURV=0 

DC 32 L=LLO,LF 

IF (M. EQ.l • AND#L.GT.LCORN. AND«L.LT. LHI )G0 TC 32 
F17 = l . 

IKURVL*IKURV 

IFIL.NE.KURVE1 .AND.L*NE.KURVE2.AND»L.NE. KUR VE3. AND* l* NE* KURVEAJGO 
1T0 2 
TKUR V- 1 

IF ( M. NE. 1 )G0 TO 102 
I SND= 1 

FHC=R (L ,M,N) 

ENERS=H(L,M,N ) -( P ( L,M, N )+ . 5*( U (L , M , N I **2+V ( L, M, N ) **2+W C L , M, N ) **2 ) / 
1P(L,M,N) )/R(L,M,N) 

CALL GAS 
ISND=0 

CHECK9-U(L,M,N)/(R(L»M,N) *SND) 


IF (L. EQ.LCORN )G0 TO 99 
CHECK8=CURV2(L-1)-CURV2(L ) 

IF(CHECK9.GT.l..AND.CHECK8.GT.O. )GC TC 100 
IF(CHECK9.LT.-1..AND.CHECK8.LT.0.)G0 TO 1C1 
IF(IYREFC.NE.O)GO TC 102 
IKURV=~l 
GC TO 3 

99 IF(CHECK9.LT« 1*)F17=0« 

100 IKURV=3 
CC TC 3 

101 IKURV-4 
C-0 TO 3 

102 CHECK=(CURV2(l~l)-CURV2(L) )*( ( AM-1 • ) *CY*CELTA ( KK ,L , N ) -YREF(L ,N ) I 
IF(CFECK.LT.C.)GO TO 3 

IKUR V=2 
GO TO 3 

2 IKURV=0 

3 DO 12 1-1*3 

IFIM.EQ.l.AND.L.EQ.LHI >G0 TO AC 

IFIL.NE.LLC.AND. I.EC.2. AND. IKURVL.NE. 1 )G0 TC 10 
IF( I . EG.3. ANO.L.EQ*LF)GO TC 11 

IFIL. EO.l.ANC. I.EQ.DGO TC A 

TF(L,NE,L10.AND.I.EQ.1.AND.IKURV.NE.2)G0 TC 10 
AO M«N 
F 1 = 1. 

L 1 = L+ I - 2 

IFIM. EC.l.ANC.L.EQ.LCCRN. AND . I .EO. 3 ) GC TC 12 

A IF((IKURV.E0.1.AND. I.E0.3 ) .OR. ( IKURVo EQ. 2. AND. I . EO* 1 ) )G0 TO 5 
eiAt l )-Fl*UCL i,M,Nl ) 

B2A( I )= s P(Ll,M,M)+U(Ll»M»M)**2/R(Ll,M»Nl ) 

P3A( I ) = E1A( I )*V(L1,M,N1)/R(L1,M,N1) 
EAAU)*ElAU)*Viai,H,Nl)/R(Ll,M,Nl) 

B5A( I ) =B1 At I ) *H(L1 ,H,Nl ) 

GO TO 12 

5 F201 = { AM-1. )*CY*DELTA(KK,1,N1)-YREF(L,N1 ) 

ALAMP=1.+CURV2(L)*F201 

L C-L- 1 

ALAMM=l.+CURV2(LQ)*F20l 
L G=L + 2- I 

IFU.EC.l I GC TC 6 
FA=ALA M M/ALAMP 
GC TO 7 

6 fa=alamf/alamm 
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7 IF(F4»GT.l<j ) F 4 = 1 • 

F4=2. *F4/< 1.+F4) 

F5=l.-F4 

TH£TAA=2.*THETA(L)-THETA( LD-THETA (LG) 
THETAB=2**(TFETA(L)“THETA(LG) ) 

I F ( I a EQ » 1 ) GO TO 8 
RNEW=R( L 1 »M,N1 ) 

UNEW=U(Ll ? MtNl)*COS(7HETAA)-V(Ll,M,Nl )*SIN(7HETAA) 
VNEW=U(L1 tM,Nl)*SIN(THETAA)+V(Ll T M,N1 ) *COS (THETAA J 
V*NEW = ML1 *M*N1 1 
FNEW=P(L1»M»N1) 

HNEW=H( L1,M,N1) 

UCLO=UOLD l*COS ( THET A8 ) -VOLDl*SIN(THETA8) 

VOLO-UOLDl*S I N (THET A B )+VOLDl*CGS (TFETAB) 

GO TO 9 

8 FNFW-ROLD 

UNEW=UaLDl*CCS(7HETAA)-VOLDi*SIN(THETAA> 

VNEW=UOLDI*S I N( THETAA )+VOLDI*COS( THETAA) 

V»NEfc = WOLO 
FNEW=POLD 
HNEW=HOLD 
RCLO=R(LG,M,Nl) 

UOLD=U(LG,M f Nl)*COS(THETAB)-V(LG,M,Nl >*$IN( THE TAB) 
VOLD=U(LGtH,Nl )*S IN (THETA B) +V( LG »M ,N 1 ) *CCS ( THETA 8 ) 
V»OLD=W(LG,M t M) 

FCLD=P(LG,M,N1I 
HOLD=H(LG.M,Nl ) 

9 eiA( I)=F4*UNEW+F5*U0LD 

B2 A ( I )=FA*( PNEW+UNEW**2 /RNEW)+F5*( POLD+UCLC** 2/ROLO ) 
R3A ( I ) = FA*UNEfc*VNEH/RNE W+F 5*UOLD*VCLD/RCLD 
P4A ( I >=F4*UNEW*WNEW/RNEW+F5*U0LD*WCL0/R0LD 
R 5A ( I }=F4*UNEW*HNEW+F5*U0LD*H0LD 
CC TO 12 

10 IG=l+l 

R 1A ( I )=B1A(IG) 

P2A ( I )=E2A ( IG ) 

P3A( I)=B3A( IG) 

PA A ( I )= BAA( IG) 

P5A( I)=65A( IG) 

GO TO 12 

11 E1A(3)=2.*B1A(2)-B1A(1) 

B2A(3)=2.*B2A(2)-B2A(1) 

B3A(3)=2.*B3A(2)-B3A(1) 


P4A(3)=2.*B4A(2>-B4A(1) 

F5A(3)=2.*B5A(2)-fl5A( 1) 

12 CONTINUE 

IF(Ll.NE.KURVEl.AND,Ll.NE.KURVE2.AN0.Ll.NE.KURVE3aAN0.Ll.NE.KURVEA 
1 )GC TO 13 
PCLD=R( L,M*N) 

CCLD1=U(L,M,N) 

\»0LD1=V(L»M,N) 

WOLD=Vr ( L * M* N ) 

FCLD= P ( L » H » N ) 

FCLD=H( L,M,N) 

13 F6=.5/0X 

IF(IKURVaGE.3)F6=2.*F6 
IF(IKURV.EQ.-1)G0 TO 2A0 
F7=F6 

IF ( L. EQ* 1 • AND, IYREF1.NE.0 ) F6=l ./OX 
I H~3 

TF( IKURV.EC.3 ) IH=2 
IL= 1 

IF ( IKURV.EQ.A ) IL=2 
16 CIFB1 = F6* (B1A ( IH) -B1A ( IL) ) 

CIFB2 = F7*(B2A( IH)-B2A< IL 1 ) 

DIFB3-.F 6*(B3A ( IH)-B3A(IL) ) 

C IFB4=F7* ( BAA ( IH)-BAA(IL) ) 

CIFB5-F6*(B5A(IH)-B5A(IL) ) 

IF(L.NE.1.0R.NF.LE.2)G0 TC 17 
N11=NF-1 

0IFBl=DIFBl+F7*(W(2,Mt?)-W(2,M f Nll ) )/OPHI 

CIFB2 = DIFB2+F7*(U( 2,M »1)**2/R( 2* P« l) -U( 2 « P »NF ) **2/R (2« M v NF)+( U(2 9 m 
1,2)*W(2*M»2>/R(2,M,2>+U(2,F,N11)*W(2,M,N11)/R(2,M,N11)) /DPHI ) 
OIFB3=DIFB3+F7*(V(2»Ht2)*Vi(2,M,2)/R(2fM,2)-V(2,M,Nll)*W(2,M f Nil) /R 
l(2»M,Nll) )/DPHI 

CIFB5=OIFB5+F7*(W(2,.H,2)*H(2,H,2)-h( 2.H, Ml )*H(2,M V Nil) ) /DPHI 
IF(F17.GT..1.0R.NF.LE.2)GC TO 17 
L AX=L— 1 

ALAX=1.-CURV2 ( LAX ) *YREF < L » N ) 

RLAX=RR EF (L )-YREF(L,N)*SIN( THETA (L ) ) 

V»LAX = W(L,H,N) /R(L»M,N) 

L'LAX-U(L*H»N)/R(L*M V N) 

CIFW=WLAX*( (WLAX*(01A(L,3)-01A(L, 1) )-CAA(L,3) +OAA ( L ♦ 1 ) ) / ( 2. *DPHI * A 
1LAX*ULAX)-WLAX*C0S(TFETA( t ) l/RLAX) 

ClFBl=DIFBl-*DIFVi/ULAX 

CIFB?=DIFB2+2.*DIFW 
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CIFB5=DIFB5 + 0IFVi*H(LffM f NI /ULAX 
17 F8=«5/(DY*DELTA(KK f LtM ) 

IF(MoEQ<,l )F8=2.*F8 
I H=3 

IF(M.EG.1)IH=2 
IL= 1 

20 IFdKURV.EQ.l.aR* I KL‘RV. EQ. 3 ) GO TO 21 
C!FC1 = F8*(C1A(UIH, NJ-C1A (L, IL,N) ) 

CIFC2=F8*<C2A<L»IH,N)-C2A(L» ILiNU 
CIFC3-F8*(C3A(L»IH*N)-C3A<L»IL*N)) 

CIFC4=F8*(C4A<L f IH*N )-C4A(L, IL,N )) 

DIFC5=F8*(C5A(L,IH,N»-C5A <L, IL*N) ) 

F9=« 5/DPH I 

F10=F9 
GC TO 22 

21 IF(L«EQ«KURVE1 ) J J J= 1 
IF(L*EC«KURVE2)JJJ«2 
IF(L.EQ.KURVE3)JJJ=3 
IF(L*EC*KURVE4IJJ,<=4 

DIFCl=Fe*(ClAC{JJJ,IH,N)-ClAC( JJJ, IL,M ) 

CIFC2=F8*(C2AC(JJJ » IH,N )-C2AQ( JJ J * IL «N) ) 

CIFC3=F8*(C3AC< JJJ, IF,N)-C3AC< JJJ, IL,N) ) 

CIFC4=F8*(C4A.G( JJ J* IH,N)-C4A0< JJJ,IL*M ) 

CIFC5=F8*(C5AC(JJJ, IH, N )-C 5AO ( JJ J , IL,NM 
IF(W»EQ*1*AND*L#EQ» LCORN) GC TO 22 
IF(NF„LE.2)GC TO 22 
l M 1=L- 1 
fvM 1=N- 1 

f^Pl=N4i 

IF(N*EQ« 1)NM1*2 
IF(N.EC.NF)NPl=NF-l 
Y= ( AM— 1 • ) *DY 

F9=.5*( 1.4CURV2<L U)*(Y*OELTA(KK,L,NPl>-YREF<L ,NP1 1 I I / ( (1.+CURV2U 
1)*(Y*DELTA(KK,L,NP1I -YREF ( L , NP1 ) ) MDPHI ) 

F 10** 5* ( 1.+CURV2I LM1 1*( Y*DELTA CKK ,L*NR1 )“YREF ( L * NMl )) )/( (l.+CURV2( 
1L)*(Y*DELTA(KK,L,NM1) -YRE F (L ,NM1 )) )*OPHl ) 

22 CIFD1=F9*01A(L,3)-F10*D1A(L,1> 

C!FD2=F9*D2A(L,3)-F10*D2A(L, 1) 

CIFD3=F9*C3A(L,3)-F10*D3A(L, l) 

CIFD4=F9*C4A(L,3)-F10*04A(L f 1) 

CIFD5=F9*C5A(L,3)-F1O*05A(L,1) 

IF(IK«JRV*EQ.-1)G0 TC 241 

Y= ( AM— 1 # ) *0Y 


LC=L 

IF C IKURV.EQ.1.CR.IKURV.EQ.2)LQ=L-1 
F14=CUR V2 CLQ ) 

ALAMB=1<»4F14*{ Y*DELT A ( K K, L , N ) -YREF (L f M ) 

PAD=RREF(U + (Y*DELTA(KK,L ,N)-YREF(L f NJ )*S INC THE TAR) ) 
IFfL.EC.l.OR.lYPEFl.EQ.OIGC TO 23 
F11=ALAMB*PDPPH(KK,L, N) /RAD 
F12=ALAMB*CCS ( THETA < L >) /R AD 
F 13= AL A MB+SI N C THETA ( l ) I /R AC 
GO TO 24 

23 f 1 1 = Oc 
F 12=0 • 

F 13=F 14 

IF<IYREF1.EQ.0)F13=0. 

24 El=(PCPT(KK,L,m*ALAMB*R(L,M,N)4PDPX<KK,L,M*UCL,M,imFl 1*W(L,M,N) 
1) /DELTA CKKtL,N)+F12*UILtM,N)+F13*V(L,H,N) 

E2=( (El+F14*V(L,M,N) ) *U( L , M, N J -F 12*W { L , M » N >**2 ) /R< L ,M,N ) +PDPX { KK , L 
1,N)*P<L ,M,N)/ DELTA <KK,L,N> 

E3=(El*V<L,M,N>-F13*V<(L,M,M**2-F14*U<L,M,N)**2)/R(L,M,N>-F14*P(L , 

1M, N ) 

F4=<El4F12*U( L,M,N)4F13*V(L,M,N) ) *M L » M, N ) /R ( L * M , N ) 4F 1 1 *P ( L , M , N ) /D 
1ELTA ( KK »L » N ) 

E5=E1*H (L ,M,N)-PDPT(KK,L, N ) *ALAMB*P <L , M, |v ) /DELTA (KK ,L,N> 

GO TO 242 

240 ALAMR= 1 » 

L J^L- 1 

ALAMM=l ,-CURV2(LJ) *YREF<L*K) 

ALAHP=1 »-CURV2(L)*YREF( L» K ) 

D IFB1=F6*< (R1A<3)-B1A(2) I /ALAMP+ ( B l A ( 2 1-e 1 A { 1 ) ) / AL AMM I 
CIFB2 = F6*UB2A(3)-B2A<2 » ) /ALAMP4 ( B2A( 2 I-B2A ( 1 ) ) /ALA MM) 

CIFB4=F6*( (B4A(3)-R4A<2) ) / ALAMP+ ( B4A ( 2 )-84 A (1) ) / AL A MM ) 

C1FB5 = F6*((B5A(3)-B5A(2I I / ALAMP4 ( R 5A ( 2 ) -E 5A ( 1 ) ) /ALA MM) 

F8-C3 A ( L,2,N)/{DY*0ELTA(KK,L,N) ) 

C!FCl=R(L f M,M*F8 
CIFC2=U(L,M,N »*F8 
riFC4=W(L»M,N)*F8 
CIFC5=DIFC1*H(L,M,N) 

F9=.5/(CPHI*ALAMP) 

F 1 0-F 9 
C-0 TO 2 2 

241 FAQ=RREf m-YREFU , N )*S IN ( THETA ( L) ) 

F 30= COS ( THET A ( L ) ) /RAD 

I F ( IYREFl.EC.9)F3O=0. 


38 



f l*F39*U(L f M,Nl 

E2=F30*(U(L*M f N»**2-W(l_,M,N|**2l/R <LtR,N I 
E4 = 2.*F30*U(L,M,NI*W{L,M,M/R<L,N t M 
£5=E1 *H (L,M,N> 

?A2 IF ( KK.EGs 2 IGG TC 25 
A10=ALAMB*R (L ,M,N) 
t?S = ALA RB*Ut L *M f N I 
A 30- ALA MB* V ( L » M t N1 
A4C=ALAMB*W( L ,M,N ) 

A5C=A1C*H(L »M ,N)-ALAMB*P( L»M,N) 

GO TO 26 

25 AIC=ALAMB*A1(L,M,N1 
A2P=ALAMB*A2( L »M,N) 

A 3C= AL AMB* A3 ( L » M » N I 
A40=AL&MB*AM L,M,N) 

A50~ AL AMB*A5(L f MtN) 

26 ALAMB*l.+CURV2(LQI*(Y*0ELTA(KKK,Lf N)-YREF(L»N) ) 
IF(IKURV«EQ*-l)ALAMB~l. 

P(L,M,N>=< A10-CT*(0IfBl+F17*(DIFCl+DIF0l4El) 1 ) /ALAM8 
l (L»M,N)=<A2C-CT*(DIFB2+F17*(CIFC2+DIF024E2I I I /ALAMB 
WtL v M v N)*( AAQ-OT*(OIFB^+DIFCA+OIFDA+EA) ) /ALAMB 
IF { M« EC» 1 J GC TC 27 

V<L»M t NI-{A30-CT* (QIF83+0 IFC3+0I FC3+E3 II /ALAMB 
GO TO 30 

27 F15=PYBPX(LiNI /ALAMB 
IFIL.EQ.1IG0 TC 28 
FI6=PYBPPH<L»N)/RA0 
GO TO 29 

28 FI 6=0* 

29 VIL,M,Nl*FI5*Ua,M, N I +F16*W( L »R » N) 

30 ENERT=( A50-DT*(DIFB5 + FI7*(0IFC5+DIF05 + E5 II l/(R(L*M t M* ALAMB} 
ENERS=ENERT-.5*(U(L,M»N)**24-V(L,M f M**24-W(L,M,N|**2)/R(L ,M,NI**2 
IF{ IGAS.EQ.O.AND. < M # NE. 1. CR.L. NE.LCORM IGC TO 31 

I SND= 1 

RHO=RtL,M,N) 

CALL GAS 
I SN0=0 

IFCM*NE.1*OR.L.NE.LCORN)GC TO 31 
GAMOG=G AMM 

BETASQ=lUIL,M,M/(R(L*M f N)*SNO) I **2-1. 

31 P(L*M,NI=( G A M M- 1 * I *ENERS*R (L ,M ,N I 
Ff L fM,N)*ENERT+P(LtM» M/RCL»M»N) 

IF(M*NE»1.0R.L.NE.LHI IGC TC 32 


LHf M1=LHI-1 

AL-AMB 1= 1*-CUR V2( LHI Ml )*YREF ( LHIM1 » N) 

F17=PYBPX(LHIM1,NJ/ALAMB1 
F18=1./SQRT( 1. 4F 1 7** 2 I 
THETA8= AT AN ( F 17 1 
Fl9=SQRT((GAM0C+l.»/(GAMG0-ni I 
IF(BETASQ.LT«1«E-09)G0 TO 210 
F 20= SORT ( BET ASCI 
F 2 1=F 20 /F 19 

ANU=F19*ATAN( F2 1 1 -AT AN ( F2 0 I 
GC TO 311 

310 ANU=0 » 

311 AMLHI = 2*/( IGA MOO-1 • ) * (THETA( LCGRNI -THETA (LHIM1 1 +THE TA8+( F 19- 1. 1*1. 
157C796A-ANU) I 

EXPF=1./(GAH00t1. ) 

F22=2.*(GAMOC-l.)*H(LCORN t M,N)/(2.4(GAMOC-l. )*AMLHI**2I 
F23=(F22*R(LCaRN f M,N)/(GAMC0*P(LCCPN,M,N) I )**EXPF 
RtLHIMl ,M,N)=R(LCORKM f N) *F23 
FtLHIMl » M , N I = P( LCOR N t M , Nl *F23**GAMC0 
KLHIM1 ,M f N) = H(LCORN, M,N) 

L (LHI Ml *M,NI=F18*SQRT(F22 1 *AMLHI*R ILHIM1 , M , M 
V(LHlMl,M t NI = F17*UaHlMl* M,N> 

W(LHIMl',M f NI = W (LCOR A » M f N) *F23 

32 CONTINUE 
LL0=2 

IFfN.FG.llGC TO 33 
PHI=PHI +DPHI 
F(l,M f N»-R(l,M f l) 

LU f M,N|=U(l • M, 1)*CCS(PHI> 

V( 1,M,N| = V( l*M,n 
M 1 ? M » N ) =— U { 1 ? M 1 1 I+SIMPHI) 

P(ltM*N) =P ( 1 » M 1 1 ) 

MltMfNI*H(l f H # U 

33 CONTINUE 
RETURN 
END 



SUBROUTINE RESULT 

CCMMCN IGAS, I SND, J , K , KK , K URVE 1 ,KURVEZ,KURVE3, KURVE4 r LF , M ,MF ,KM,N, N 
IF f ILC, DENSO f V ELC, KP f KPP ,KF ,LH I , LCCFN 

COMMON ALPHA, 3 S, C SALP, CUR VE1, CURVE 2, CURVE3, CUR VE4, CURVES, CXO,DPH I, 
1CT.OX ,DY,EINF, ENERS, EPS I, GAMM, GAMMA, P INF, PM IN,RHC»RS, SNALP , SND , XO , 

* 2VREFn,YREFi ,YREF 2, YREF3, YREF4 

CCMMCN A1!17»6»9)»A?( 17*6»5),A3! 17,6,9), A4( 17,6,9) ,A5(17,6,9) , H ( 1 7 
1,6,9) ,P11T, 6, 9) »R(17 f 6, 9) ,U( 17,6,9) ,VC 17,6,9) ,W! 17, 6,9) ,C1A( 17,3, 9 
2 ),C2 A (17, 3,9) ,C3A( 17,3,9) ,C4A( 17,3,9) ,C5A { 1 7 , 3 ,9 ) ,C 1 AO (4 , 3,9 ) , C2 AO 
2 ( 4 ,3 , 9 ) ,C3AC!4,3,9),C4A0(4,3,9) f C5A0!4,3,S) 

COMMON C1A{ 17,3) ,D2A(17,3) , C3A ! 1 7 ,3) ,D4A ( 1 7 , 3 ) ,D5A( 17,3), CURV2 ! 17 ) 
1, THETA ( 17),RREF(17) , XREF! 17) ,YREF( 17,9),PY8PX( 17,9) ,PYBPPH( 17,9) ,D 
2FLTA(2, 17,9) , PDPT (2 , 17 f 9 ) , PDPX <2 , 1 7 ,9 1 , PDPPH ( 2, 17, 9 ) 

COMMON AAA (5,5) » AAB! 5 ) , B1 A ( 3 > , B2 A l 3) ,B3AC3) , 84A ( 3 ) , B5A I 3 ) ,RRXF (5) , 
1XRXF (5) ,CT,TIME, IYREFO, IYREF1, IY RE F 2 , I YR EF 3 , l YREF4 
CCMMCN/ ERROP/ER( 12) ,8LK(12) 

CC MMON/ BLOCK l /K PO 

900 FORMAT!* GRIC*//4H LF= t 3, 2X3HM F= 1 3 ,2X3HNF= 13/ / 

1 AH DX=E16.8,2X3HDY=E16. 8, 2X5HDPHI*E16. 8/// 

2* REFFRENCE-SURFACE PROFILE*// 

3FH KURVE1=I3,2X7HKURVE2=I3,2X7FKUPVE3=I3,2X7HKURVF4=I3// 

A eH CURVEl=Elfc.8,2X7HCURVE2=E16.8,2>7HCURVE2=E16. 8, 2 X7HCURVE4=E16. 8 
5,2X7HCURVE5=Fl6.8/) 

931 FORMAT! 3H l* I 3 ,2X 9HTHE TA ( L ) *E 16. 8 , 2 X9HCUP V2 ( L I = E16. 8 , 2X8HXREF ( L ) = E 
116*8, 2X8HRREF (L)=E16.8) 

992 FORMAT!//* BCDY SURFACE*//7H YREFC=E16 .8 , 2X6HYR EF1= 616 . 8 ♦ 2X6HYREF2 
1-E16* ft, 2X6HYRGF3=E16.8» 2X6NYREF4=E16.8/) 

903 FORM AT (3H L= l 3 ,2X2 HN= 1 3 , 2 X 10HYREF !L » N ) =E 16 . 8 , 2X1 1HP YBP X ( L ,N ) = E 16. 8 
1,2X12HPYBPPH(L,N)=E16.8I 

904 FCPMAT!//* INITIAL SHOCK*// 

14H RS=E16.8,2X3HBS=E16.8, 2X3HXC=E 16. 8 , 2X4HCX0=E 1 6. 8// / 

2* GAS MODEL PARAMETERS*// 6H IGAS=IE// 

37H DENSC=E16.8,2X5HVEL0=E16.8/// 

4* NCNDt MENS ION AL FLIGHT CCNDI T IONS*/ / 

57H GAMMA = E16* 8 »2X5HG AMM=E 16. 8 , 2X 5HP INF=E 16 . 8 , 2X6HAL PHA*E 16. 8/ // 

6* PRINTING AND CUTOFF PARAMETERS*// 

9 5H KPO- I5,2X4HKPP=I5,2X3HKF*I5/// 

8* COMPUTATIONAL PARAMETERS*//* I YREFQ=*, 1 3// 

S AH CT=E16.8,2X5HEPSI=E16. 8,2X5HPMIN=E16.8////) 

905 FCPMAT (41 5, 6E 16. 8/ (20X6E16.8) ) 

906 FCPMAT { 20X6F16.8) 

907 fCRMATl/) 


90 8 FORMAT! 4XIHK', AX1HL , 4X1HM, 4X1HN ,7X2HUN, 14X2HVN, 14 X2HWN, 13 X3HRH0, 1 1 X 
18HP(L,M,NI,8X8HH(L,M,N)/26X3HPHI *13X3 HR AC, 13X3HXA0, 12X4HGAMM, 1 1X6H 
2ENTP0P,11X5HAMACH/22X12HDELTA( 1,L,N) *4X11HPCPT! 1*L* N ) ,5X1 1HPDP X! 1 , 

3L,N), 5X12HPDPPH! 1,L,N )////) 

909 F CRM AT ( * TIME**,E16.8/// I 

913 FORMAT ( 8H I YREF2 = I 3 , 2 X7HI YR6F 3 = 1 3 , 2X 7HI YRE F4* I 3/ / ) 

911 FORMAT!//* TYPE OF FLOW*//* I YREF 1=*, I 3// ) 

912 FORMAT!//) 

913 FORMAT! 1H1//) 

914 FORMAT (3H L = I 5 ,2X2HN= 1 5 , 2X 13FDELT A ( l ,L ,N ) = E 16 . 8 , 2X4HRAD=E 16. 8, 2X4H 
1 XAO=E 16.8) 

PRINT 913 

ALPHAD= ALPHA* 57. 2 5577951 

TF,(K.GT.O)GO TC 3 0000C568 

PRINT 9 11, IYREF1 

PRINT 900,LF,MF t NF,DX,DY,DPHI , KUR VE1 ,KURVE 2 ,KURVE3 , KURVE4 ,CURVE1 , C 
IURVE2, CURVE3, CURVE4,CURVE5 
LHIMl=LHI-l 
CO 1 L* 1, LF 

THETAD-THETA(L)*57.25577951 

1 PRINT 9C1,L,THETAC,CURV2!L) ,XREF(L) ,RREF(L) 

PRINT 902,YREFO,YREFl,YREF2,YREF3,YREF4 
PRINT 910,IYREF2, IYREF3,I YREF4 

CO 29 N»1,NF 
CC 2 L=1,LF 

2 PRINT 903,L ,N ,YPEF(L,N) ,PY8PX(L,N) ,PYBPPH(L,N) 

PRINT 9C7 

?C CONTINUE 

PRINT 9G4,RS,BS,X0,CX0, IG A S , CENSC , V6L0 ,G AMMA , GAMM , P INF , ALPHAD , KPO 
IKPP,Kf , IYREFC,CT,EPSI,PMIN 
PRINT 913 
IF (K* EQ.O IGG TC 3 
CO 39 N=i»NF 
CO 30 L=l ,LF 

F1=0ELTA( 1,L»NJ-YREF(L,N) 

P AC =RREF(l)+Fl*SIN( THETA! L) I 
XAD=XREF!L)+F1*C0S( THETA! L) I 
39 PRINT 914,L,N»0ELTA( 1 ,L,N) ,RAC,XAO 
RETURN 

3 PRINT 909 , T I M E 
PRINT 9C8 

CC 6 N= l, NF 
AN=N 
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PHI=57, 2957795 1*< AN-1. )*OPHI 

CO 5 F*1 f MF 

AM=M 

DO 4 L= 1 * LF 

IF<M.EQ.1.AND.L.GT.LCGRN.AND.L.LT.LH!P11GC TG 4 
F 1= C AM-1» ) *DY *DELTA ( 1,L,N l-YREFC L *N) 

PAD=RREF(L»*F1*SINCTHETAUU 
XAD=XREF(L)+F l*COS( THETA ( L ) ) 

GN*U(LtR f NI/RCL,M,KI 
VN=V(L,P»N) /R ( L»M* N ) 

ViN-W(L»F*M/F(L*M*N) 

1 SND= 1 

PHO=R(L,M,N) 

F2=UN**2+VN**2+WN**2 

ENERS= c H(L»M*N)-P(LtPfN)/RCLfK»N)-«5*F2 
CALL GAS 

AMACH=SQRT(F2 l/SND 

EMRQP=P(L»M, NI/A8S(R<L»M,N) )**GAMH 

PRINT 905 *K*L » M»N*UN»VN*WN*RHO»P (L»M T N)*H(L»M f NI f PH I , RAD , XAD, GAMM 
1FNTR0P t AMACH 

IF ( EQ*MF) PR INT 906 , DELT A ( 1 , L t N ) , PDPT ( I # L » N ) ,PDPX( i,L,N> ,PDPPH( 1 
1L,N) 

4 CCMINUE 

IF(P.NEoMF) PR INT 907 

5 CONTINUE 
PRINT 912 

6 CCNTINUE 
RETURN 
FND 



SUBROUTINE SHOCK 

COMMON IGAS, l SNO, J , K ,KK ,KURVE 1 ,KURVE 2, KUR VE 3, KUR VE4 ,LF * M , MF ,MM ,N , N 
IF , LLO, DENSO* VELO,KP, KPP*KF,LHI,LCCFN 
COMMON ALPHA, B$,CSALP, CURVE 1* CURVE 2, CURVE 3, CURVE4, CURVE5 ,CXO,DPHI , 

. lCT,DX,DY,EINF,ENERS,EPSl,GAMM,CAMMA,PINF,PKIN,PHO,R$,SNALP,SND,XO, 
2YREFOtYREFl,YREF2,YREF3,YREF4 

COMMCN Al(17, 6,9),A2( 17,6,9),A3( 17, 6,91, A4( 17, 6,91, A5(l 7,6*9) , HI 17 
It 6, 9) ,P( 17, 6,9), R< 17, 6,9) ,U( 17, 6,9 I, V( 17,6, 91, W< 17, 6,91, C1M 17,3,9 
2),C2A( 17,3,9) ,C3A { 17,3,9) ,C4A( 17,3,9) ,C5A( 1 7, 3 ,9) ,C 1 AO( 4 , 3,9 ) ,C2AO 
3 (4,3,9) ,C3AO{4,3,9 ) ,C4A0< 4, 3, 9 ) ,C5AO(4 ,3 ,9 ) 

COMMON D1A(17,3),D2A(17,3) ,03A ( 17 , 3 ) , 04 A ( 1 7 ,3 ) ,D5A( 17,3 ) ,CURV2(17 ) 
1,THETA( 17),RPEF(17) ,XREF( 17) ,YREF( 17 , 9 ) , P YB PX ( 1 7 , 9) , P YBPPH( 17 ,9) , D 
2ELTA I 2,17,9) , PCPT ( 2 , 1 7, 9 ) , PDPX (2 , 1 7 , 9 ) , PCP PH( 2 , 17, 9 ) 

COMMON AAA <5, 5 ) , AAB I 5 ) , B1 A ( 3) ,B2A( 3) ,B3A(3) ,B4A{ 3) , B5A (3 ) ,RRXF (5 ) , 
1XRXF15) ,C7,T IME,IYREFC,1YREF1, IYREF2, I YR EF 3 , I YREF4 
CCMMON/ERROR/ER(12) , BLK { 1 2 ) 

900 FORMAT ( * SUBROUTINE SHOCK *, 1A6, 41 5/ ( BE 16. 8 ) » 

KKK= 3-KK 
I SND= 1 
M-MF 
MM 1=M*-1 
MM 2=M-2 
LLO= 1 

CO 66 N- 1 , NF 
DC 65 L-LLO, L F 

IF(L.NE.KURVEl.AND.L.NE.KURVE2.AND.L.AE.KUPVE3.AND,L.NE,KURVE4)GC 
1TC 41 
IKllRV-1 

CHECK-(CURV2(L-1)-CURV2(L ) )*( DELTA IKK ,L, M-YREF ( L ,N ) ) 
IF(CHECK*LT*0* )GO TC 4? 

IKURV=2 
GO TO 42 

41 IKURV=0 

42 ALAMB = l*-*-CURV2(U*(OELTA(KK,L,M-YRFF 
SNTHE=SIN(T HE TAIL) ) 

CSTHE=CCS(THETAIL) ) 

F AC=RREF(L) + ( DELTA! KK,L ,N)-YREF(L,K) ) * SN THE 
IFIL.EC.DGC TO 43 

F1=(PDPPHCKK,L,M+PYBPPHCL,N) ) /RAD 
CO TO 44 

43 F 1-C • 

44 F2=(PCPX(KK,L,N)+PYBPX(L, N) )/ALAMB 


F 3= SORT I l.+Fl**2+F2**2) 

CSRX=-F2/F3 

CSBY-1. /F 3 

CSBP~-F 1/F3 

IF ( I KUR V* KE* 1 ) GO TO 45 

ALAMB=1.+CURV2(L-1 ) * (OELTA IKK ,L , N ) -YREF ( L , A ) ) 

F2 = (PCPX<KK,L ,M+PYBPXtL,NU/ALAKB 
F3=SQRT ( l.+Fl**2+F2**2) 

45 FHO=R(L,M,N) 

ENEPS = H(L,M,M-(P(L,M,N)+(U<l,M,N)**2+V(l,M,N)**2+W<L,M,N)**2> /RIL 
1 , M ,N ) ) /R ( L , M , N ) 

CALL GAS 

VNl=tCSBX*U(L,M,N>+CSBY*V<L*M,K)+CSBP*W<l,M,N) )/R(L,M f N) 

VN.2=!CSBX*U(L,MM1,NDCSBY*V(L,MMI,N)+CSBP*W(L,MM1,N) ) /R (L ,MM1 , N) 
VN3=(CSBX*U(L,MM2,N)+CSBY*V(L,MM2,NI+CS8P*W(L,MM2,N) )/R (L,MM2,N) 
LG=L 

IFIIKURV.EQ.l >LG=L-1 
IFIKK.EO. 1) A4(L,M,N)=P(L,M,N) 

IF IKK* EC* 1) A5(L,M,N)=VN1 
6HL,M,M=R(L,M,N)*SNO 

A2IL , M , N ) =A4 ( L ,M , N ) *A 1 ( L , M , N ) * A5 I L , M , N ) -CT / IDELTA I KK,L , N ) *OY ) * (F 3* 
1 1 VM+SNOI-POPT IKK, L,N) ) *. 5*1 3. *P<L ,M,N )-4. *PI L ,MM1 , N ) +P( L ,MM2, N) +A 
Z1(L,M,N)*(3**VN1~4.*VN2+VN3> ) 

A3(L,M,N) =-0T * SND*CURV 2(LG)/ALAMB*( SND* VI L , M , N ) -U( L ,M,N)/R(L,M,N)* 
1(CSBY*UIL,M,N)-CSBX*V(L,M,NI ) ) 

IF(L*EQ.LF)GC TO 46 
IF I IKUR V* EQ* 1 ) GC TO 49 
LPl=L+l 

VN5=(CS8X*U(LP1,M,N) ♦CSBY*V( LP1,M,N)+CSBP*V((LP1,M,N) )/R( LP1,M,N) 
LTX5=(UCLP1 ,M,N)-CS8X*(CSBX*UI LP1 , M, N )+CS0Y*V( LP1 ,MI,N )+CSBP*W( LP1 , 
1M,N) ) )/RILPl,M,N) 

46 IFIIKURV*E0.2)GG TO 49 
IFIL.EC.DGC TO 47 
LM1=L-1 

M=N 

F4-1* 

GO TO 48 

47 LM1=2 
M=NF 
F 4=- 1 • 

48 VN4=f F4*CSBX*U(LM1,M ,N1 ) + CSBY*\iI LM1,M,N1 ) + CSBP*W I LM 1 , M, N 1 ) )/R(LMl , 
IM,N1) 

L7X4=(F4*U( LM1 ,M,N1)-CSBX*(F4*CSBX*U(LM1,M,M )+CSBY*V(LMI,M,Nl)+CS 
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lepsVULNWtSNU ) 1/RCLPlfWtMI 
60 TO 54 

49 l l=L+3-2* IKUR V 
LQ=L+1-IKURV 
CHECK=ABS (CURV2C LG I l 
T F (CHECK. LT . 1 . E-09 ) CC TO 50 

PRl=RREF(L)+< DELTA ( KK ,L , N l-YREF ( L , N) ) *SI N( THETA { U ) 
XXl=XPEF(L)+(DELTA(KK f L»N l-YREFCL, N) )*COS(THETA( L) > 

RR3=RREF( L1) + (0ELTA(KK»L1 ,N) -YREF < LI ,N ) ) *S IM THETA { L 1 ) ) 
XX3=XREF(Ll)+t DELTA (KK,L1*M-YREF(L1*N))*C0S( THETA( LI) ) 

RADA=SQRT { ( PR1-PR3 ) **2+ (X XI— XX 3 1**2) 

IF( IKURV.EQ. 1)F5=1. 

IF ( IKURV.EQ. 2 >F5=-1. 

R ADB= 1. /CURV2(LG)+DELTA(KKfL»M)-YREF(LfN) 

RR2 s =RRl-RA0B*SIN(THETA{L) ) 

XX2=XX1-RAD8*CCS(TFETA(L) ) 

RADC=SQRT((RR2-PR3)**2+(XX2-XX3)**2) 

IF(CUPV2(LG) .GT.0.1F6=1. 

IF(CURV2(LG) . LT.0.IF6=-1. 

CTHETA=F5*F6*AC0S(F6*(RADB**2+RADC**2-RACA**2> /( 2.*RADB*RADCl ) 

THE TA A= THET A ( L ) +DTHET A -THET A ( L 1 1 
IF(IKURV.EQ.1)L2=L-1 
IF(IKURV*EQ.2 )L2=L + 1 
F7=(THETA(L)-THETA(L2 ) J/OTHETA 
GO TO 51 

53 TFETAA = THETA(L )— THET A (LI I 

P ADB= 1, /CURV2(LQ)+BELTA(KK»L1»N)-YREF(L1»N) 
RACC=RACB*ABS(SIN(THETAA) ) 

F7=0X/RA0C 

51 IF(F7.GT. 1. )F7=1. 

UNEW=U(L1 ,M,N)*COS ( THETAA )-V ( L 1 ♦ M f N ) *S IN < TH ETA A ) 

VNE Vi = U ( LI » Mt N ) *SI N < THETAA ) + V (LI » M » M *COS (THETA A ) 

VNN=(CSEX*UNEW + CSBY*VNEW+CSBP*M Ll,M f M ) /R (L1«N«N) 

L'TXX= { UNE V*-CS BX* ( CSBX+UNE W+CSBY*VNEW+CSBP*W(Li»M f N) ))/R(Ll.M f N) 

IF l IKURV.EQ* 2) GO TO 52 

VN5-VNN 

LTX5=UTXX 

GO T 0 46 

52 VN4=VNN 
LTX4=UTXX 

54 IP1=L+1 
l M 1 = L- 1 

IFIL.LT.LFIGC TO 55 

CIF2=P(L»M f N)-P(LMl f M,N) + Al(L,M ( N)*(VM-VN4) 

LTX1=U(L,M,N)/R(L ,M,N)-C5BX*VN1 

CIF3=LTX1-UTX4 

GO TO 59 

55 F 8=« 5 

IF ( IKURV.EQ. 1. OR. IKUR V. EQ.2 ) F8=F7/ 1 1.+F7 ) 

M=N 

IF ( L. NE.l ) GO TO 56 
N1=NF 
L V 1 = 2 

56 0IF2=F8*(P(LPltM,N)-P(LMl*N,Nl)+Al(LfF f N)*(VN5-VN4) ) 

C IF3=F8*(UTX5~UTX4 ) 

59 A2(l»M,N)=A2<L,MfN)-DT*DlF2/CALAMB*DX)*(U(L,M»N)/R(L,M t NKCSBX*SND 
,1) 

A3(L#MtN)»A3(L»«tN)-DT*DlF3/C ALAf»B*DX )*R (L . N, M )*SND**2 
1FUYREF1.EQ.OIGO TO 65 
IF (N.NE.l.ANO.N.NE.NF)GO TC 63 
CAP-0. 

IF (L. EQ.l )GC TC 62 
IF( N# EC.NFIGO TC 60 
N 1=2 
F8= 1 • 

GO TO 61 
63 N 1=NF- 1 

1F(NF.EC.1)N1=N 

F8=-l. 

61 CIF3=F8*(1.-CSBP**2)*W{L*H,N1)/RCL,M,M) 

DA3=-DT*SND**2/PAD*(R<L»M, M *0 IF3/0PHI+V U * M # N ) *STN (THETA(L) )+U(L, 

IR *N) *COS(THFTA(L ) ) ) 

GC TC 64 

62 M=NF~1 

IF( NF. EC. 1 )N1=N 

DA3=-DT*SND**2/ALAMB*<R(LfPf M *( U (2 , M, 1 ) /P ( 2, M* 1 )+U ( 2,M, NF) /R( 2,M, 
INF ) ) * • 5/DX+V ( L » Kf N I *CURV2 ( LG ) ) 

IF(NF.LE.2IGG TC 64 

CA3=OA3-OT*SNO**2/ALAMB*( l.-CSBP** 2) *( W( 2,N,2l /R(2, M , 2 >-W (2 , M, Nl )/ 
1R(2*M»N1 H/DPFI *.5/DX*R(LtM,N) 

GO TO 64 

63 KP1=N+1 
NM1=N-1 

VN7=(CSBX*U(L ,K,NP1)+CSBY*V(L,M, NPi)+CSBP*W(Lf M,NPl > )/R (L.M« NP1 ) 
VN6=(CSeX*U(L»M f NMl ) +CSBY *V( L , M , NM 1 | +CSBP*W ( L , M , NM1 ) > /R < L ,M» NM1 ) 
lTP7=(W<LtM,NPl)-CSBP*(CSBX*U(L,M^Pl )+C$BY*V< L*M,NP1 ) +CSBP*W ( L* M , 



1NP1I ) l/R(L t M f NPl)*RIL,M,N) 

UTP6=(W(L»M»NMl)-CSBP*(CSex*U(LtM,NMl >4CSBY*V( L ,M,NM1 ) +CSBP+W ( L , M , 
IN Ml) ) )/R(LjM»NMl)*R(L»M T N) 

CA2=-.5*OT/RAC*(W(L, M t NI/R(L,M t N)4CSBP*SNC)/DPFI*(P(L,M f NP1)-P(L»M 
ltNMIl +A1(UM,N)*( VN7-VN6) ) 

CA3=-0T*SNC/RAD*(SND*{ UTP7— UTP6 1 * • 5/0PHI+ { SND*V(L»M*N) — W { L» M * N ) /R( 
1L,M,N)*(CSBY*W(L»M, N)-CSBP*V(L,M, NH )*SIN(THETA(L))4(SND*U(L,M,N)- 
2W(L,M f N)/R(L*M f N)*(CSBX*W(L»M,N)-CSBP*U(L,M,N) ) ) *C0S( THE TA( L > ) I 

64 A2(L»M»N)=A2(L»M,N)4CA2 
A3(L,M t N)=A3(L,M,N)4CA3 

65 CCNTINUE 
LLC=2 

66 CCNTINUE 
LLC= 1 

CO 4 N=1,NF 
AN=N 

PHI=( AN-1, 1*0 PHI 
5NPHI = S IN ( PH I I 
CSPHI -CCS ( PHI ) 

CO 3 L=LLOfLF 
SNTHE«SIN(THETA(L> I 
C$7HE*C0S(THETA(U ) 

C SAX-SNTI-E*CS ALP+GSTHE*CSPHI * SNALP 
CSAYs— CSTHE*C SALP + SN THE*CS PH I *SN AL P 
CSAP=-SNPHI*SNALP 

ALAMB=l.4CURV2(L )*( DELTA! KK,L,N) -YREFIL,N 1 ) 

P AD*RREF ( L ) ♦ ( DELT A(KK*L*N)-YREF(L»N) )*SNTFE 
IF(L.EQ.l)GO TC 1 

F1-CPCPPH{KK,L,N)+PYBPPH(L ,N) ) /RAD 
GO TO 2 

1 F1=0. 

2 F2 = (P0PX(KK»L»M + PYBPX(l»N))/ALAHB 
F3=SQRT( l.+Fl**2+F2**2) 

CSBX=-F2/F3 

C SB Y= 1, / F3 
CSBP-— F 1/F3 

CSSHK=CSAX*CSBX+CSAY*CSBY4CSAP*CSPP 

ES=H(L, YF,N)-(P(L,MF, N)+. 5*(U(L,MF ,N»**24V(L,MF,N)**2+WIL,MF,N)**2 
l)/R(L,MF,Nn/R(LtMF,NI 
I TERU=4 

IF ( IGAS.EC.O) ITERUM 
no 32 I TER=1 , ITERU 
PHO=R (L *MF»N) 


ENER S=F S 
CALL GAS 

IF { ITER #GT, 1 1 GC TO 30 
F4=A1 (L,M,N) 

F5 = A2(L»M»N) + A3(L»M#M 

GC TO 31 

33 IF C ITER.EQ.ITERUIGO TO 32 

31 F6= (GAMMA-GAP M )/( GAMMA-l. )*PINF/ (F5- ( CAMP- 1 . ) / ( G AMM A-l . ) *P I NF-F4*C 
1SSHK ) 

F9=( l.-F6 + F6 3 *<2l*F4**2+l. l.-F6)*( (GAMM+1. ) * ( F 5-F 4*CSSHK) 4( GAMM- 

11. )*PINF) 

F 10= ( 2 • * ( 1.4F6**3)-3.*F6*( l.+F« ) ) *F4**2*4. 5 *( 1. -F6**2 ) * ( G AMM41 . ) * ( 
lF5-F4*CSSHK)-4.5*( ( 5. *GAMP4l . ) -6 • *G APP*F64 ( GAPM-1 . ) *F6**2 )*PINF 
IF(F9.GE.1.E— Cfc) GO TO 311 

PRINT 900,BLK( 1) , K , L , M, N, F4 , F5 ,CSSHK , CSBX , CSBY ,CSBP , F 9,F 10 
CALL RESULT 
STOP 1101 

311 Fll=, 5* F4*F1 C/F9**l • 5 

IF(ABS( Fll) .LE.l. )G0 TO 310 

PRINT 900 * BLK (2 ) ,K,L*M,N» F4,F5 ,CSSHK , CSBX, CSBY, CSBP ,F9 ,F 10,F11 
CALL RESULT 
STOP 1 1C2 

310 FSA=ACCS(Fin/3. 

F7= ( SQR T( F9 ) * ( COS ( PS A ) 4 1. 732C 50 8*S IN ( PS A ) )-F4*(l.+F6) )/3. 

F12-P INF/F7**2 

F8=F7*(GAPM*(F12+1. ) -SORT ( (GAMM*F12-1. I **2+2.* (GAMM+l. )* (GAMMA-GAM 
IP) *F 1 2/ (GAMMA-1.) ) ) /(GAMM+1. I 
P(L,MF,N)=F7/F8 

ES=(F8/F7*PINF+F8*(F7-F8) )/(GAMM-l. ) 

32 CCNTINUE 
PHO=R ( L *MF »N ) 

F(L»MF t N) = (GAMM-1. ) *RHO*ES 
U(L,MF,N)=(CSAX+(F7-F8) *CSBX)*RHO 
V(L,MF,N)=(CSAY4(F7-F8)*CSBY)*RHC 
W(L,MF,N)=(CSAP4(F7-F8)*CSBP)*RH0 
F9=1.4(F7-F8) *(2.*CSSHK4F7-F8) 

ML,MF,N) = ES4.5*F94p(L,MF,M/RH0 
kS=(CSSHK4F7)*CSBY 

DELTA (KKK,L, N )=OELTA ( 1*L,N)4.5*DT*(WS4PDFT (1,L,N) ) 

PDPT(KKK*L,N|=kS 

IF( KK* EQ. 1 )GC TC 3CC 

F10=SQRT(F9) 

F11=CSA X*C$BY— CSAY*CSBX 
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F12=CSAP-CSSHK*CSBP 
F 13=1 * 

F 14= 1 a 

IFtFll.LT.O. J F13=-l • 
IF(F12.LT.O. )F14=-1. 
PDPX(KKK,L,N)=F13*F10/SND 
FDPPH(KKK f L,N) =F14*F 10/SND 
GO TO 3 

300 PDPX(KKK # L,M=POPX<KK,L»N) 

POPPHtKKK,L,M=POPPH(KK,L ,NI 

3 CONTINUE 
L LG= 2 

IF(N.EQ.11G0 TC 4 
R(l»MF,Nt*R(lt*F,l) 

Ul.HFf NJ«U( 1 *MF* l)*CSPHI 
V ( 1 » MF t N )=V C 1 » MFf 1 1 
V«< 1, MF ,N)=-UC 1 ,MF » 1 ) *SNPH I 
P{l f MF»N)=PC ltFFfl) 

K 1 ? MF * N ) =H ( 1 * MF » 1 ) 

CELT A ( KKK , 1 » N ) =DELTA { KKK, l , 1 ) 
PDPT(KKK,1,N)=PDPT(KKK, It 1) 
IF( KK«EC« 1 )GO TC 40C 
P0PX(KKK, 1»N|=PDPX(KKK, 1, 1) 
FDPPH<KKK,1 ,N)=PDPPMKKK, 1,1) 
GO TO 4 

400 PDPXtKKK, 1,N)=PDPX(KK,1,N) 

FDPPH(KKK,1 ,M=PDPPH<KK,1,M 

4 CONTINUE 

IF ( KK • EC«2 ) C ALL SLOPE 

RETURN 

END 



SUBROUTINE SLOPE 

COMMON IG AS » I SND , J , K »KK , KURVE1 ,K URVE 2 » KURVE 3 ,KURVE4 , LF , M , MF , MM ,N , N 
IF, LL3, DENSO, VELC,KP, KPP ,K F ,LH I »LCCF N 
COMMON ALPHA, BS,CSALP , CUR VEl,CURVE2»CURVF3,CURVF4,CURVE5,CX0tDPHI , 
IDT, OX ,DY,EINF,ENER$,EPSl,GAMM, GAMMA, P INF, PV IN, RHC,RS,SN ALP, SND,XO, 
2YREFO,YBEFl,YRFF2,YREF3,YREF4 

CCMMCN A1(L7,6,9) , A2 ( 17,6,9) ,A3( 17,6,9) » A4 ( 17,6,91 ,A5( 17,6,9) ,H( 17 
1,6,91 ,P(17,6,9) ,R(I7,6,9) ,U(17,6 ,9 > , V ( 17,6,9 ) ,W ( 17, 6, 9 » , C 1A( 17, 3, 9 
2 ),C2A( 17,3,9) ,C3A( 17, 3,9) ,C4A( 17,3,9) ,C5A( 17,3,9) ,C 1 AO ( 4 , 3 , 9 ) , C2 AO 
3(4,3,91 ,C3AC<4,3,9) , C4 AO( 4, 3 ♦ 9 ) ,C5 AQ ( 4 , 3 , 9 ) 

COMMON C1A( 17,3) ,D2A( 1 7,3 ) , D3 A ( 1 7 , 3 ) ,C4A(17,3) ,D5A( 17,3) ,CURV2(17) 

1 , THET A ( 17),RREF(17) ,XREF( 17 ) , YREF ( 17 , 9 ) , PYB PX ( 17,9) , PYBP PH( 17,9) , D 
2ELTA(2,17,9), PDPT (2,17,9) ,PDPX (2 , 17 , 9 ) , P DPPH ( 2 , 1 7,9 ) 

COMMON AAA (5, 5 ) , A AB ( 5 ) , B1 A ( 3 ) , B2 A ( 3 ) , B3A { 3 ) ,B4A(3) ,B5A(3) , RRXF (5) , 
1XRXF (5) ,CT,T (ME, IYRFFC, IYREF1, IYREF2, IYREF3, l YREF4 
CCMMCN/ERRQR/ER(12) ,BLK(12) 

IF(K,EQ.O)GC TO 1 
KKK=3-K K 
GO TO 2 

1 K KK= 1 

2 DC 13 N=1*NF 
CO 12 L=l,LF 
CC 6 ILG=1 , 5 
LLQ=L-3+I LQ 

IF(L.EG.l«OF« ILQ.EC.5IGO TO 3 

I LQP 1= I LQ+ 1 

PRXF ( ILQ)=RRXF(ILQP1) 

XRXF ( ILQ)=XRXF(ILCP1 ) 

CO TO 6 

3 IF(LLQ.GT.LF) GC TO 5 
IF ( LLQ* LT* 1 1 GO TO 4 
M=N 

F 100= 1* 

GO TO 40 

4 M = NF-N+1 
LLQ=2-LL0 
F 100=- 1 • 

40 FRXF (ILQ) =RRE F (LLQ )+ (DELTA (KKK,LLQ,N1 ) -YP F F ( LL C , N1 ) ) *S IN ( THETA ( L L 0 - 
1) » 

FRXF ( ILC)=F100*RRXF ( ILQ) 

XRXF ( I LQ) =XREF ( LLQ ) + 1 DELT A (KKK ,LLQ , N1 ) -YREF ( LLQ ,N1 ) )*COS (THET A(LLQ 
1 ) ) 


GO TO 6 

5 lLQM2=ILQ-2 

RRXF( ILQ)=2.*RRXF(ILQI-RRXF(ILCM2) 

XRXF( ILQ) = 2.*XRXF ( ILQ)-XRXF( ILCM2) 

6 CONTINUE 
IF(K*EQ.0)GC TC 8 

IF(PCPX(KKK,L,N)*GE.-U.ANQ.POPX(KKK,L,N).LE, l.)GO TO B 

!F(PDPX(KKK,L ,M*LT.-1. )GC TC 7 

TGTME=- (XRXF( 3 J-XRXF ( 2) ) / ( RRXF ( 3 ) -PR XF ( 2 ) ) 

CO TO 9 

7 TGTME=- ( XRXF( 4)-XRXF(3) ) / ( RRXF (4 ) -PRXF (3 ) ) 

GO TO 9 

9 TGTME=t XRXF (4) -XRXF (2) )/( RRXF (4) -FRXF (2) )- ( XRXF( 4 )-XRXF( 3))/(RRXF( 
14 )— RR XF ( 3 ) )-( XRXF ( 2 ) -XRXF ( 3) I / ( RR XF ( 2 )-RP XF ( 3 ) ) 

9 TGE=( SIN(THFTA(L) )-CCS(THETA(L)J*TGTME)/(CCS(THETA(L))*SIN( THE TAIL 
1 ) ) *TGTME ) 

PCPX(KKK,L,N)=TGE*( l.+CURV2(L )*( DELTA (KKK, L,N)-YREF (L,N) ) )-PYBPX(L 
1 ,N) 

IF(K.EQ.O)GC TC 11 
IF(N#EC*l*OR*N*EC*NF)GO TO U 

IF(PDPPMKKK,L,N) *GE*-1*. AN'O. PDPPH (KKK ,L , N ) . LE* 1 . ) GO TO 11 

JF(PDPPMKKK,L,N) .LT.-l. )GC TO 10 

NM1=N-1 

IF(NM1.LT.I)NM1=2 

PDPPH(KKK,L ,K)= (DELTA (KKK,L,N) -DELTA (KKK,L,NM1 ) ) /DP HI 

GO TO 12 
13 NP1=N+1 

IF(NPl.GT.NF) NP1=NF-1 

PDPPH ( KKK, L , N »=( DELTA (KKK ,L , NP 1 ) -DELTA ( KKK ,L,N) ) /OPHI 
GC TO 12 

11 NP 1=N+ 1 

IF(NP1.GT.NF)NP1=NF-1 
1F(NF.LE.2)NP1=N 
NM 1=N— 1 

!F(NM1,LT.1)NM1=2 

IF(NF*LE*2)NM1=N 

PDPPH ( KKK,L, N ) = .5*( DEL TA(KKK,L,NP1 )-DELTA( KKK, L,NM1) ) /DPHI 

12 CONTINUE 

13 CONTINUE 
RETURN 
END 
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SUBROUTINE START 

COMMON IGAS, I SND , J , K , KK , KURVE1 ,KUR VE2 ,KUPV E 3 , KUR VE4 , LF , M , MF , MM ,N , N 
IF, LLG, DENSO, VELO,KP,KPP,KF,LHl ,LCORN 

COMMON ALPHA»BS»CSALP»CURVE1,CURVE2 , CURVE 3 »CURVE4,CURVE5»CXO»DPHI , 
1DT,DX,DY,EINF , ENER S » EPS I , G AMM » GAMM A , PI NF , PM IN, RHO, R S , SN ALP , SND , XO , 
2YREFO,YREFl, YREF2,YREF3,YREF4 

COMMON A1U7, 6,91, A2 (17,6, 9 1, A3t 17, 6,91, AM 17,6,9) ,A5( 17 ,6,9) ,H( 17 
1,6,9) ,P<17,6, 9),R(17,6,9) ,U< 17 ,6 ,9 ) , VI 17 ,6 ,9 ) ,W( 17, 6,9) ,C1A( 17,3, 9 
2),C2A(1?,3,9) , C3A( 17,3,9) ,C4A< 17, 3 ,9 ) »C5 A ( 17, 3 , 9 ) ,C 1ACI 4 , 3, 9) , C2 AO 
3 (4,3,9) ,C3AC( A, 3, 9) ,CAAO( A, 3, 91 ,C5A0(A,3,S ) 

COMMON C1A( 17,3) ,02A( 17,3 ) ,D3A( 17,3 ) ,D4A ( 17,3) ,D5A( 17,3) ,CURV2(17) 
1 , THET A ( 17 ),RREF(17) ,XREF{ 17 ) , YREF < 1 7., 9 ) , PYBPX ( 17,9) ,PYBPPH< 17,9) ,D 
2ELTAI 2,17,9) , PCPT < 2 , 17, 9 ) , PDPX ( 2 , 17 ,9 ) ,POPPH< 2, 17,9 ) 

CCMMON AAAI 5, 5),AAB( 5),BIA(3),B2A( 3) , B3A (3 ) , BAA ( 3) , B5A ( 3 ) ,RRXF(5) , 
1XRXFC5) ,CT,T I ME, I YREFC, IYREF1, IYREF2, I YREF 3 , IYREFA 
COMHC N /ERR OR/E R ( 12 ) ,BLK(1 2 ) 

900 FCPMAT ( * SUBROUTINE START * , 1A6 , A 1 5/ ( 8E 16, 8 ) ) 

AMF=MF 

ANF=NF 

IF(NF*ECcl)ANF-2« 

CY=1* /( AMF-1* ) 

CPHI=3« 1A 1592654/ ( ANF-1. ) 

FINF*PINF/(GAMMA-1* ) 

LCCRN^O 

LHI=0 

7IME=C. 

DO 1 L= 1 » LF 

IF ( L* LT « KURVE 1 )CURV2 ( L ) -CUR V El 

IF(L*GE*KURVEl.AND*L.LT.KURVE2)CURV2(l )=CURVE2 
IF ( L.GE.KURVE2*AND*L*LT.KURVE3 )CURV2(L )= CURVE 3 
IF(L.GE.KURVE3.AND»L*LT.KURVE4)CURV2IL)=CUFVE4 
IFIL.GE.KURVEA )CURV2 CL )=CURVE5 
1 CONTINUE 
THETA! 1 )=0.0 
FREF C 11=0. 

XREF ( !)=YREFO 
TC 3 L= 2 , L F 
L M1=L — 1 

THETA{L)sTHETA(LMl)4CURV2(LMl)*DX 
IFCABS(CURV2(LMl)).GT.i.E-C9)GC TO 2 
PREF (L)=PREF(LM1)+C0S( THE TAIL) )*DX 
XREF(L)=XREF (LMl) -SIN (THE TAIL) )*DX 


GO TO 31 

2 RRFF(L)=RREF(LM1) + !SIMTHETA(L) ) -SIN (THETA (LMl) ) )/CURV2(LMl) 
XRFFI L) = XREF( LMD + (COS (THETA (L ) ) -COS ( THE T A ( LMl ) ) )/CURV2f LMl) 

31 F99« ABS ( 1 ®-YR EFC*CURV 2 ( LM 1 ) ) 

TF(LCORN.NE.O)GC TO 30 
IF (F99.LT. l.E-C6)LCCPMLMl 
GO TO 3 

30 I F ( LH I • NE • 0 ) G 0 TO 3 
IF(F99*GE.i*E-C6ILHI=L 
TF(L*EQ*LF)LHI-99 

3 CONTINUE 

CC 5 N=1,NF 
• AN=N 

PHI={ AN-1* )*DPFI 
DO A L= 1 » LF 
YRFF(L,N)= YREFC 
A C.CNTTNUE 

5 CONTINUE 

CO 12 N= 1 ,NF 
DO 11 L=1 ,LF 

IF { N *NE • 1 *AND *N,NE *NF ) GO TC 6 
F YBF PH ( L , N ) =0 • 

GO TO 7 

6 NP1-N+1 
NM1=N-1 

PYBPPH(L»N )=— *5*(YRFF(L»NP1)-YREF(L ,NM1) ) /CPHI 

7 IF(L*NE*1 )G0 TO 8 
YREFPl-YREFI 2 ,N) 

NXY-NF-N+ 1 

LP1=2 
LM1-2 
GC TC 10 

8 LM 1=L- 1 
L P 1 = L + 1 
NXY=N 

IF ( LP 1 • L E *LF ) GC TO 9 
YRrFPl=2**YREF(L,N)-YREF( LM1,N ) 

GO TO 1C 

9 YRFFP1=YREF(LP1,N) 

10 FYBPX(L,N)=-*5*(YREFP1 -YREF ( LMl , NXY ) )/DX 
IF(LP1.E0*LHI )PYBPX«L,N)=2.*PYBPX(L,N) 

11 CONTINUE 

12 CCNTINUE 



SNALP=S IN(AIPHA) 

CSALP=COS! ALPHA) 

CC 18 N* 1 ,NF 
AN=N 

PHI*( AN- 1 . )*CPFI 
5 NPHI =S I N ( PH I ) 

CSPHI =-COS ( PH I ) 

CO 17 L* 1 ,LF 
SNTHE=$IN(THETA!L) ) 

CSTHE=CCS(THETA(L) ) 

R 1 =RRFF !L)-YPEF!L,N)*SNTHE 

Xt=XREF<L)-YREFa,N)*C$THE-YREFG + YREF<l,l) 

fi=cstfe*snalp-snthf*csphi*csalp 
F 2 =CSTFE*CSALP+SNTHE*CSPH I*SNALP 
F 3 -R 1 *C STFE-I Xl+XO)*SNTHE 
F 4 =(RS-BS*CXC)*CSTHE 

F 5 -SNPH l** 2 +C SPHI** 2 * ( C $ALP** 2 +B S* SNAL P **2 ) 

F 6 =F 1 ** 2 +(SNTFE*SNPHI )** 2 4 BS*F 2**2 

F 7 =F 2 *F 4 - < (F 1 *CSALP-F 2 * 8 S*SNALP) *CSPHI— SfsTFE*SNPHI **2 ) *F 3 
F 8 =F 3 ** 2 *F 5 + 2 • *F 3 *F 4 *C$PH l*SNALP- ( RS*C$TFE*F 4 ) *CXO*CSTHE 
IF!ABS(F 6 ).LT„ 1 .E- 09 )GO TC 13 
F 9 =SQRT(ABS( F 7 ** 2 -F 6 *F 8 ) ) 

IF(CSTHE«LT« 0 #)F 9 --F 9 
X 2 = (F 9 -F 7 ) /F 6 -XO 
GO TO 130 

13 X 2 =-* 5 *F 8 /F 7 -X 0 

133 !F(ABS(CSTFE).LT.. 01 )GO TC 14 
F 2 =R 1 +$NTHE*( X 2 -X 1 )/CSTHE 
GO TO 16 

14 F 10 -SNPHI** 2 +(CSPHI*CSALP)** 2 *BS*(CSPHI*SMLP )**2 

FI 1 = ( BS-U)*(X 2 +X 0 )*CSPHI*SNALP*CSALP + (RS-B$*CX 0 )*C$PHI*SNALP 
F 12 -<X 2 +XC)** 2 *(SNALP** 2 +eS*CSALP** 2 ) 42 .*(RS-BS*CX 0 )*(X 2 +X 0 )*CSALP 
1 - 2 * *RS*CX 04 BS *CX 0**2 
IF(ABS(F 10 ).LT. 1 .E-C 9 )G 0 TC 15 
F 13 =SQRT ( ABS ( F 1 1 ** 2 -F 10 *F 12 ) ) 

P 2 =(F 13 -Flt)/F 10 
GO TO 16 

15 P 2 =-« 5 *F 12 /F 1 1 

16 DELTA( 1 ,L,N)=SCRT( (X 2 -X 1 ) ** 24 (R 2 -R 1 )** 2 ) 

CELTA( 2 f L»N ) = DELTA ( l , L , N) 

FDPT ( 1 , L , N) =0 • 

17 CONTINUE 

18 CONTINUE 


K = 0 

CALL SLOPE 
C$RDMX=0. 

LLO= 1 

DO 22 N= 1 » NF 
AN=N 

PHI = ( AN-1 • )*CPFI 
SNPHT =S IN ( PHI ) 

CSPHI = COS < PHI ) 

CO 21 L = LLO# LF 
$NTHE=SIN(THETA(L) ) 

CSTHE=CCS (TFETA(L) ) 
C$AX=SNTHE*CSALP4CSTHE*CSPHI*SNALP 
CSAY=-CSTHE*CSALP*SNTHE*CSPHI*SNALP 
CSAP=-SNPHt*SNALP 

ALAMB=1.4CURV2IL)*( DELTA! 1 ,L , N ) -YREF ( L ,N ) ) 
RAD=RREF(L)+( DELTA! 1 , L , N ) -YREF (L ,N )) *SNTFE 
IF(L«EQ«1)G0 TC 19 

F 14- ( PDPPH! 1,L,N)+PY8PPH! L,N) > /RAD 
GO TO 20 

19 F 14=0 • 

20 F15=! PDPX!1,L,N)*PYBPX(L,N) )/ALAHB 
FI 6- SORT { 1« + F14**2 + F1 5**2 ) 

CS8X=-F 15/F16 

CSBY=1./F16 
CSBP=— F 14/F 16 

CSSHK=CSAX*CSBX+CSAY*CSBY+CSAP*CSBP 
IF!C,SSHK.LE.-1.E-06)GC TO 201 
PRINT 900, ER( 1) ,K , L , KF , N , CSSHK 
K = - 1 

CALL RESULT 
STOP 1301 

201 IFU.GE.LCORMGO TO 200 

F98=-PYBPX(L,N)/l l.-YREF<L,N)*CURV2!L> ) 
F97=l»/SQRT(l#-fF98**2)' 

CSBD= { C SAX*F9 8+C SA Y ) *F97 
IF(CS8C,GT.CSBDHX)G0 TO 200 
CSBDM X=CSBD 

200 F 17= 1 « + P I NF/C S SHK**2 

F18=l. + 2.*GANMA/( GANN A- 1. ) *PINF/CSSHK**2 
F19=SQRT(F 17**2-! l.-l. /GA N K**2 ) *F1 8 ) 

P(L, MF t M = <GAMP4l. )/! G AMM * { F 17-F 19 )) 
IF(R!L,PF,N).GE.1.)GC TO 202 
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PRINT 900,ERI2).K,L,FF,N» P(LtMF,M 
K =- 1 

CALL RESULT 
STOP 1302 

232 l<L»MF,N)=R<LtNF r N)*(CSAX-CSSHK*CSBX)+CSSHK*CS8X 
Va*MF t N)=R(L T MF,N>*(CSAY-CSSHK*CSeYJ+CS$HK*CS8Y 
MLt MF,M=R(L f FF t N)*{CSAP-CSSHK*CSEP)+CSSHK*CSfiP 
F(L?HF»N)=PINF+CSSHK**2*( 1 •- 1. /R ( L * MF * N ) ) 
ML»MF,N)=.5 + GANMA/(GAMMA-I. J+PINF 

21 CONTINUE 
LLC=2 

IF ( N« EQ . 1 ) GC TO 22 
R(1«MF«N!*R< 1 V MF, II 
UltMF,N)=U(l*f'F t l)*CSPHI 
V(1»MF»N)*V< 1,PF,1) 

M1»MF,M*-U( 1 1 KF» 1 ) *SNPH I 
P(1*MF»M«P(1 9 FF V 1) 

M1,MF,N)=H(1,MF,1) 

22 CONTINUE 

F 17*1 »+P INF 

F18=1.+2.*GAMMA/(GAMNA-1. )*PINF 
F19=SQRT(F17**2-( l.-l. /GANM**2 )*F18) 

FNS=(GAMM41* ) /(GAMM*(F17-F19) ) 

FNS— PINF+1® — 1 * /PNS 
AMSQ=1. /( GAPH*RNS*FNS ) 

FXPG-GAPM/(GAMM-1« ) 

F20=l#+« 5*1 GANN-1. )*AMSQ 
PN5T=PNS*F20**EXP0 
EXPC=EXPO/GAPM 
RNST«RNS*F20**EXP0 
FNST=GANM/(GAHM-1. )*PN$T/RNST 
EXP0*1 • /GAMM 
CO 27 N=1 »NF 
Afs = N 

PHI = ( AN— 1* )*DPHI 
SNPHI =S IN ( PHI I 
CSPHI=COS (PHI ) 

CO 26 L*ltLF 
SNTHE=SIN(THETA(L) I 
CSTHE-COS ( THE TA ( L ) ) 

CS AX= $NTFE*CS ALP+CSTF E*CSPHI*SNALP 

CSAY=:-CSTHE*CSALP4-SNTHE*CSPHI*SNALP 

CSAP=-SNPHI*SNALP 


L AL=L 

IF(L.GE.LCORN.AND.LCORN.NE.O)LAL=LCORN-i 

ALAMB=1.-CURV2(LALJ*YREF( L,NI 

PAO=RREF(U-YREF(L*N»*SNTHE 

IFU.EC.11G0 TO 23 

F21=PYBPPH(L t N)/RAD 

CC TO 24 

23 F2l=0. 

24 F22=PYBPX(L,N)/ALAMB 
F23=SQRT( l.+F21**2*F22**2 ) 

CSBX=-F22/F23 

CSBY=1./F23 

CSBP=-F 21/F23 

C$BD=C$AX*CS8X+CSAY*CSBY+CSAP*CS8P 
IFCL.GT.LCCRMGC TO 240 
P(Lfl»N)=PNST*(C$BD/CSBDMX)**2 
GC TO 241 

243 F(l, 1 tN)=PNST*CSBD**2+PINF*( l.-CSBD**?) 

IF(LCORN.EQ.O)GC TO 241 

IF(P(Lt ltM .GT.PUCCPN, 1, N M P ( L , 1 , N ) = P( L CORN , 1 f N ) 

241 L M 1=L - 1 

IF(CSB0.GT.0*0)P<L, 1,N)=P(LM1, 1,N) 

IFCPILf lfN).LT.PMIN)P(Ltl ,N)=PPIN 
P(Lt 1 ,N)=IP(L ,XtN)/PNST)**EXPO*RNST 

LT=SQRT (ABS(2.*(HNST-GAMM/(GAMN-l. )*P(L, 1 , N ) /R ( L » 1 t N Ml ) 

F 24= SORT! APS < l.-CSBD**2) ) 

IFIF24.LT. 1.E-C9 I F24=l. 

UT =UT /F 24 

L (1_ 1 1 » N )=R (Lt 1»M*UT*(CSAX-CSBD*CSBX) 

V { L» 1 iN )=R ( L 1 1 ♦ N) *UT * (CSAY— CSBC*CSBY ) 

MLt ltN) = R(L, 1,M*UT*(CSAP-CSBC*CSBP) 

IF ( L.NE.LCORN )GG TO 242 
WCHK=W(L,l,N)/R(L,i,M 
FCHK=HN$T-. 5*WCHK**2 
E XPCK=G ANP/ IGAPM—l.) 

PCHK=PNST*(2.*HCHK*(GAMM-1. )*RNST/(GAMM* (GANM+1. ) *PNST'I) * *EXPCK 
TF(P(L,1,N).LE.PCHK )G0 TO 242 
P ( L » 1 1 N ) =PCHK 

R (Lt 1*N )=RNST *( PCHK/PNST) **EXPC 
LT=SQRT (GAMM* PCHK/R ( L » 1 » N ) ) 

L ( L » 1 .N)=R(L, 1*N)*UT*C$BY 
V(Ltl»N)=-R(L ,ltN)*UT*CSB* 

ML,1 ,N»=ViCHK*R(L,l ,M 
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2 Ml HL,1 ,NI=FNST 
f*UP=MF-l 
CD 25 M=2,MUP 
AM=M 

Y=( AM-1. ) *DY 

R<L»M ? N)«Y*R <UHF f N) + < l.-Y)*R<Lt 1»M 
L(L,M,N)»Y*U(L,FF,NI+Cl.-YI*U(Ltl,N) 

V(L,M,N )=Y*V(L,MF f N)+< i.-Y)*V< L*1*N» 

W f L » M f N )=Y*W (LfMFfN) + fl.-*Y)*W(L»l»NJ 
P(L,M t N) = Y*Pa,MF,N)M l.-Y >*PC L* 1 , M 

ML*H* N)=(GAMM/(GAMM-1. I *P <L » M # N ) + . 5* < U( L» M, N ) ** 2+V (L ,M , Nl **2+W( L , 
1N»N)**2)/R(L»M»N) J/R(L*M»N) 

25 CONTINUE 

26 CONTINUE 

27 CONTINUE 
RETURN 
END 
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SUBROUTINE TRANSC 

COMMON I GAS, I $NC» J *K ,KK ,KURVE 1 , KUR VE2 , KUPVE 3, KURVE4»LF,M,MF,MM»N,N 
IF »LLG »OENSO » V ELC ,KP » KPP ,KF ,LHI ,LCORN 

COMMON ALPHA* BS,CSALP»CURVEl»CURVE2 , CURVE 3 »CURVE4,CURVE5 ,CX0,DPHI, 
1CT,DX,CY»EINF» ENERS, EPS I * GAMP * GAMMA* P INF, PH IN * RHO*R S, SNALP , SND , XO , 
2YREFO,YREFl,YREF2,YREF3,YREF4 

CCMMCN A1 (17, 6**5) , A2( 17,6,9) ,A3( 17, 6,9), A4( 17,6, 9) , A5( 17 ,6,9) *H< 17 
1,6,91 ,PU7,6,9),R<17,6,9) ,U( 17,6 ,9 1 ,V( 17,6, 9) ,W< 17, 6,9) ,CiA( 17,3,9 
21 ,C2A ( 11 , 3 , 9) ,C3A( 17,3,9) ,C4A( 17,3,9) ,C5A( 17,3,91 ,C1A0(4 ,3,9) ,C2AO 
3 <4,3,9 ) ,C3A0(4»3»9) , C4 A0( 4,3,9) ,C5AO(4,3»9 ) 

COMMON C1A(17,3)»D2A(17,3),03A(1?,31,C4A(17,3),05A( 17, 3 ) , CURV 2 ( 17 ) 
1,THETA( 17),RREF(17),XREF< 17) , YREF ( 17,9),PYBPX( 17,9) ,PYBPPH( 17,9) ,D 
2ELT A ( 2,17,9) , PORT (2, 17,91 ,PDPX (2 ♦ 17 , 9 1 , PCPFH ( 2 , 17, 9 ) 

COMMON AAA ( 5 » 5 ) , A A B ( 5),B1A(3) ,B2A( 3) ,B3AC3) ,B4A(3) ,B5A(3) ,RRXF{5) , 
1XRXF ( 5 ) ,CT,TIMF,IYREFC, IYREF1 , IYREF2, I YR EF3 , I YREF4 
CCMMCN/ERROR/ER(12) , BLK ( 1 2 1 
IJ=J-1 
LIC=1 

CO 3 N= 1 , NF 

CO 1 L= LLC , LF 

C1A(L,J,NI=C1A(L,IJ,N) 

C?A(L, J,N)=C2A(L, IJ,N) 

C3A(L,J,N)=C3ACL, IJ,N) 

C4A(L,J,N)=C4A(L, IJ,N) 

C5A(L»J»N)=C5A(L» IJ,N> 

1 CCNTINUE 
LLC= 2 

IFIKURVE1.GE. LF1GC TC 3 

IF (KURVE1.LT. LF)JJJH=1 

IF(KURVE2.LT.Lf)JJJH=2 

IF(KURVE3.LT.LF)JJJH-3 

IF(KURVE4.LT.LF)JJJH*4 

DO 2 JJJ=1,JJJH 

ClAOl JJJ,J,N)=C1A0< J JJ, IJ ,N) 

C2AO( JJJ,J,N)=C2AC(JJJ, IJ,N) 

C3AC< JJJ,J,N)=C3AOUJJ, IJ,M 
C A AO( JJ J , J , M = C4A0 ( J J J, IJ,M 
C5AC(JJJ,J,N) =C5A0 ( J J J , IJ,N) 

2 CCNTINUE 

3 CCNTINUE 
A RETURN 

END 



CJ1 

DO 


APPENDIX B 


SAMPLE INPUT 

The input for a sample case is presented in this appendix. The case involves the Apollo configuration traveling 
at 6935 m/s (22 750 ft/sec) and an angle of attack of 22° at an altitude of 45.87 km (150 500 ft) in the atmosphere of 
the earth. 


$NAME 



YREF3 

= 

0.1E+01, 

IYREF1 

= 

1. 

YREF4 

= 

0.1E+01, 

LF 

= 

10, 

RS 

= 

0. 1 253E+ 01 

MF 

* 

4, 

BS 

= 

0.1E+01, 

NF 

= 

7, 

XO 

= 

0. 1015E+01 

OX 

= 

0.14096201E+00, 

CXO 

= 

0. 107E+0 1 , 

KURVE1 

= 

6, 

IGAS 

= 

1. 

KURVE2 

= 

99, 

DENSO 

= 

0.1668E-02 i 

KURVE3 


99. 

VEL0 

= 

Q.6935E+04i 

KURVE4 

= 

99. 

GAMMA 

= 

0.1405E+01, 

CURVEl 

= 

0.5 7043 5 3 1 E *-00 , 

GAMM 

= 

0. 1142E+01, 

CURVE2 


0.16582451E+01, 

PINF 

= 

0.1583E-02 , 

CURVE3 

= 

o 

• 

o 

«• 

ALPHAD 

= 

0.22E+02, 

CURVE4 

= 

0. 1E«-01. 

KPO 

= 

100, 

CURVES 

= 

O.lEfOl. 

KPP 

= 

100, 

IYREF2 


99, 

KF 

= 

300, 

IYREF3 


99, 

I YREFO 

= 

1, 

IYREF4 

= 

99, 

CT 

= 

0.54E+00, 

YREFO 

= 

0 .553047 16E +00, 

EPSI 

= 

0. 3125E- 01 , 

YREF1 

= 

0. 1E*01 , 

PMIN 

= 

0.1E+00, 

YREF2 

s 

0.1E+01, 

SEND 





APPENDIX C 


SAMPLE OUTPUT 


In this appendix, the output for the calculation described in appendix B is given. This solution was calculated in 
10 minutes on the CDC 6600 digital computer. 


TYPE OF FLOW 
I YREF 1= 1 


GRID 

LF= 10 MF= 4 NF= 7 

0X = 1.40962010E-01 DY= 3.33333333E-01 DPHI= 5. 23598776E-01 


REFERENCE-SURFACE PROFILE 

KURVE1= 6 KURVE2= 99 KURVE3= 99 KURVE4= 99 

CURVE 1= 5.70435310E-01 CURVE2= 1 .65824510E+00 CURVE3= 0. CURVE4= 1 .00000000E+00 CURVE5= 1 .03300000E +00 

L= 1 THET A( L )= 0. CURV2(LI= 5.70435310E-01 XREF(LJ= 5.53047160E-0 1 RRE F ( L )= 0. 

L= 2 THET A ( L I = 4 . 6C713689E + 00 CURV2(L»= 5 . 70435310E-01 XREF!L)= 5. 47382856E-01 RREFI L) = 1 .40810156E-01 

L = 3 THET AI L I = 9 . 2 1 42 7379E + 00 CURV2(L)= 5. 70435310E-01 XREF(L»= 5 . 30426548E-0 1 RRE F l L> = ' 2 . 8071 0362E-01 

L= 4 THET A( L 1= 1 . 38214107E +0 1 CURV2(L)= 5. 70435310E-01 XREFI L 1= 5.02287811E-01 RREF(U= 4 . 18796552E-01 

L= 5 THET A( L I = 1 . 84285476E+01 CURV2IL)= 5 .70435310E-01 XREF(L>= 4 . 63 148486E-0 1 RREF ( L ) = 5 . 541 76379E-0 1 

L= 6 THET AI L ) = 2 . 303 56845E + 01 CURV21L>= 1 . 65824510E +00 XREF(LI= 4 . 1 32 61 499E-0 1 RREFI L ) = 6 . 85974986E-01 

L= 7 THET A( L 1= 3 . 64285478E + 01 CURV2(L)= l. 65824510E + 00 XREF(LI= 3 . 43 51 1075E-0 1 RRE F I L ) = 8 .08101384E-01 

L= 8 THET AI L ) = 4. 982141 12E+0 1 CURV2!L)= 1 . 65824510E + 00 XREF(LI= 2 . 47369793E-0 1 RREFI L) = 9 . 10750459E-01 

L= 9 THET AI L >= 6.32 142746E + 01 CURV2(L)= 1 . 65824510E + 00 XREF(LI= l . 30066829E-0 1 RREFI L I = 9 .88339074E-01 

L= 10 THE T A( L I = 7. 66C71380E + 01 CURV21L)= 1 .65824510E+00 XREFIL)= -2.01764858E-03 RRE FI L ) = 1 .0366471 4E+00 


BODY SURFACE 

YREF0= 5. 5 304 71 6 OE —01 YREF1= 1 . OOOOOOOOE+OO YRFF2= 1 . OOOOOOOOE +03 YREF 3= 1 . 00 000000E+00 YREF4= 1 . 000000 30 E+00 
I YREF2= 99 IYREF3= 99 I YREF4= 99 


cn 

w 



ai 


L = 

1 

N = 

1 

L = 

2 

N= 

1 


3 

N= 

1 

L = 

4 

N= 

1 

L = 

5 

N = 

1 

L = 

6 

N= 

1 

L = 

7 

N= 

1 

L = 

8 

N= 

1 

L = 

9 

N= 

1 

L = 

10 

N= 

t 


L = 

1 

N= 

2 

L = 

2 

N= 

2 

L = 

3 

N = 

2 

L = 

4 

N= 

2 

L = 

5 

N= 

2 

L= 

6 

N= 

2 

L = 

7 

N= 

2 

L = 

8 

N= 

2 

L = 

9 

N= 

2 

L = 

10 

N= 

2 


L = 

1 

N= 

3 

L- 

2 

N= 

3 

L = 

3 

N= 

3 

L- 

4 

N= 

3 

1 = 

5 

N= 

3 

L = 

6 

N= 

3 

L = 

7 

N— 

3 

L = 

8 

N= 

3 

L = 

9 

N= 

3 

L = 

10 

N= 

3 


YREF ( L » N ) = 5. 

YREFIL,N)= 5. 
YREFfL,N)= 5. 
YREF(L.N)= 5. 
YREF( L • N 1= 5. 

YREF t L t N ) = 5. 

YREF(L,N)= 5, 
YREF( L » N )= 5. 

YREF t L . N ) = 5. 

YREF(L,N)= 5. 


YREF1 L • N )= 5. 

YREF ( L • N ) = 5. 

YREF(L,N»= 5. 
YREFIL,N)= 5. 
YREFtL,N)= 5. 
YREF(L.N)= 5. 
YREF f L » N ) = 5. 

YREF ( L » N ) = 5. 

YREF(L,N»= 5. 
YREFtL,N)= 5. 


YREF t L . N ) = 5. 

YREF( l_» N)= 5. 

YREF( L »N )= 5. 

YREF(L.N)= 5. 
YREF ( L » N ) = 5. 

YREF t L • N > = 5. 

YREFIL,N>= 5. 
YREF(L,N)= 5. 
YREF ( L « N ) = 5. 

YREF < L t N ) = 5. 


530471 60E— 01 
53047160E-01 
53047160E-01 
53 0471 6 0 E-01 
530471 60E-01 
530471 60 E-01 
53047160E-01 
530471 60E— 01 
530471 6 OE-01 
53047160E-01 


530471 6 0E-01 
53047160E— 01 
530471 60E-01 
53047160E-01 
53047160E-01 
530471 60E—0 1 
530471 60E— 01 
53047160E— 01 
53 0471 60E— 01 
53047160E-01 


53047160 E— 01 
53047160E-01 
530471 60E-01 
53 047 1 60 E— 0 1 
530471 60E—01 
53047160E-01 
530471 60E— 01 
530471 60E-0 1 
530471 60 E-01 
53047160E— 01 


PYBPX t L t N )= 0. 

PYBPX(L,N)= 0. 
PYBPX(L,N)= 0- 
PYBPX ( L t N )= 0. 

PYBPX(L.N)= 0. 
PYBPX t L « N 1 = 0. 

PYBPX ( L»N)= 0. 

PY8PXtL,N)= 0. 
PYBPX(L,N)= 0. 
PYBPX(L,N)= 0. 


PYBPX(L,N)= 0. 
PYBPX C L « N J= 0. 
PYBPXtL»NI= 0. 
PYBPXtl,NI= 0. 
PYBPX(L»N)= 0. 
PYBPX t L » N ) = 0. 

PYBPX(LtN)= 0. 
PYBPX(L,M(= 0. 
PYBPXIL,N)= 0. 
PYBPXf L »NI = 0. 


PYBPX t L »N ) = 0. 

PYBPX ( L * N ) = 0. 

PYBPXtUN) = 0. 

PYBPXIL.N)= 0. 
PYBPX ( L » N )* 0. 

PYBPX ( L * N ) = 0. 

PYBPX ( L « N ) = 0. 

PY8PX(L?M)= 0. 
PYBPX t L » N ) = 0. 

PYBPX < L ♦ N I — 0. 


PYBPPHt L »N ) = 
PYBPPHt L # N > = 
PYBPPHt LtN) = 
PYBPPH(L,N)= 
PYBPPH(L.N) = 
PYBPPHtL.N) = 
PYBPPHtL,N)= 
PYBPPHtL,NJ= 
PYBPPHIL,N>= 
PYBPPHt L *N) = 


0. 

0. 

0. 

0. 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 


PYBPPH(L,N)= 0. 
PYBPPHt L ,N ) = 0. 

PYBPPH<L,NI= 0. 
PYBPPH(L,N)= 0. 
PY3PPHtL,N»= 0. 
PYBPPHIL,N)= 0. 
PYBPPH«L,N)= 0. 
PYBPPH(L,N)= 0. 
PYBPPHt L tN ) = 0. 

PYBPPH t L »N ) = 0. 


PYBPPHt L »N) = 
PYBPPH(L,N)= 
PYBPPHt L t N J = 
PY8PPHtl»N>= 
PYBPPH t L »N J = 
PYBPPHt L *N> = 
PYBPPHt L *N ) = 
PYBPPHt L «N) = 
PYBPPHt L I - 
PYBPPHt L *N ) = 


0. 

0. 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 



L = 

1 

N= 

4 

YREFIL,N)= 

5. 530471 6 0E- 01 

PYBPXIL, N»= 

0. 

PYBPPHI L »NI = 

L = 

2 

N= 

4 

YREFI L » N ) = 

5.53047160E— 01 

PYBPXIL, N>= 

0. 

PY8PPHIL,N)= 

L = 

3 

N= 

4 

YREF|L,N)= 

5. 53047160E— 01 

PYBPXIL.NI= 

0. ■ 

PYBPPHI L »N ) = 

L= 

4 

N= 

4 

YREFI L , N ) = 

5. 53047160 E- 01 

PYBPXI L , N ) = 

0. 

PYBPPHI L »N) = 

L = 

5 

N= 

4 

YREFI L,N) = 

5. 530471 6 OE— 01 

PY8PX ( L » N ) = 

0. 

PYBPPHI L .N ) = 

L= 

6 

N= 

4 

YREFIL.N)= 

5. 530 47 160 E- 01 

PYBPX IL » N ) = 

0. 

PYBPPHI L *N) = 

L = 

7 

N= 

4 

YREFI L,N) = 

5.530471 60 E-01 

PYBPX I L » N ! = 

0. 

PYBPPHI L »N 1 = 

L = 

8 

N= 

4 

YREF(L.N)= 

5. 530471 6 OE— 01 

PYBPXIL, N>= 

0. 

PYBPPHI L «N ) = 

L= 

9 

N= 

4 

YREF 1 L « N ) = 

5. 530 47 l 60 E- 01 

PYBPXIL, N)= 

0. 

P Y8PPH I L »NI = 

L = 

10 

N= 

4 

YREFI L » N ) = 

5.530471 6 OE— 0 1 

PYBPX(L.N»= 

0. 

PYBPPHI L »N ) = 


L = 

1 

N= 

5 

YREFI L » N ) = 

5. 530471 60 E-01 

PYBPX(L.N}= 

0. 

PYBPPHI L iN ) = 

L= 

2 

N= 

5 

YREFI L « N 1 = 

5 . 530471 60E-01 

PYBPXIL.N»= 

0. 

PYBPPHI L.N) = 

L = 

3 

N= 

5 

YREFI L • N ( = 

5. 530471 6 OE— 01 

PYBPXIL.N»= 

0. 

PYBPPH(L»N)= 

L = 

4 

N= 

5 

YREF 1 L » N ) = 

5. 530471 60 E-01 

PYBPX IL . N 1 — 

0. 

PYBPPHI L »N) = 

L= 

5 

N = 

5 

YREF 1 L . N 1 = 

5 . 530471 60 E-01 

PY8PXI L.N ) = 

0. 

PYBPPHI L .N ) = 

L = 

6 

N= 

5 

YREFI L » N ) = 

5. 530471 60E-0 1 

PYBPX ( L »N 1 = 

0. 

PYBPPH ( L »N ) = 

L = 

7 

N= 

5 

YREFI L »N ) = 

5.53047160E-01 

PYBPX l L » N 1 = 

0. 

PYBPPHI L .N) = 

L = 

8 

N= 

5 

YREFI L . N )~ 

5. 53 0471 60 E-01 

PYBPXI LpN 1= 

0. 

PY8PPHIL,N)= 

L = 

9 

N= 

5 

YREFI L • N )= 

5. 53047160 E-01 

PYBPXIL »HI = 

0. 

PYBPPHI L.N) = 

L= 

10 

N- 

5 

YREFI L.NI = 

5. 530471 60 E— 01 

PYBPXIL.N)= 

0. 

PYBPPHI L .N ) = 


L = 

l 

N= 

6 

YREFI L « N) = 

5.53 0471 6 OE— 01 

PYBPXIL. N)= 

O'. 

PYBPPHI L .N ) = 

L = 

2 

N= 

6 

YREFIL,N)= 

5. 53047160 E— 01 

PYBPXIL. N)= 

0. 

PYBPPHIL,N>= 

L= 

3 

N = 

6 

YREF|L.N)= 

5.530471 60E— 01 

PYBPXIL, N)= 

0. 

PYBPPHI L.N) = 

L= 

4 

N- 

6 

YREFI L • N 1 = 

5.53047160E-01 

PYBPXIL. N)= 

0. 

PYBPPH(L«N)= 

L= 

5 

N= 

6 

YREFI L.NI = 

5.53047160E-01 

PYBPXIL, N)= 

0. 

PYBPPHI L.N) = 

L = 

6 

N= 

6 

YREFI L . N ) = 

5. 53 04 7160 E-01 

PYBPXIL, N)= 

0. 

PYBPPHI L ,N ) = 

L = 

7 


6 

YREFI L . N ) = 

5.53047160 E-01 

PYBPXIL, N>= 

0. 

P YBP PHI L ,N ) - 

L = 

8 

N= 

6 

YREFI L . N 1 = 

5 . 530471 60 E— 01 

PYBPXIL. Nl* 

0. 

PYBPPHI L ,N ) = 

L= 

9 

N = 

6 

YREFI L , N ) = 

5.53047160 E-01 

PYBPXIL. N)= 

0 . 

PYBPPHI L »N) = 

L- 

10 

N= 

6 

YREFIL.N)= 

5. 53 04 71 6 OE— 01 

PYBPXIL. N)= 

0 . 

PYBPPHI L ,N ) = 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


CJ1 

CJJ 



L= 

1 

N- 

7 

YREF(L,N)= 

5.53047160 E-01 

PYBPXIL,N)= 

0. 

PYBPPHI L,N) — 

0. 

L = 

2 

N = 

7 

YREF(L,NI= 

5. 530471 60 E— 01 

PY6PXIL,N>= 

0. 

PYRP PHI L , N I = 

0. 

L = 

3 

N= 

7 

YREFI L i N ) = 

5.530471 6 OE— 01 

PYBPXIL,N)= 

0. 

PYBPPHI L , N 1— 

0. 

L= 

4 

N= 

7 

YREFI L,N) = 

5.53047160E— 01 

PYBPX(L,Nt= 

0. 

PYBPPHI L ,N) = 

0. 

L= 

5 

N= 

7 

YREFI L ,N)= 

5. 53047160 E- 01 

PYBPXIL,N1= 

0. 

PYBPPHI L ,N) = 

0. 

L = 

6 

N= 

7 

YREFIL ,NI= 

5. 530471 60 E- 01 

PYBPXIL,N)= 

0. 

PYBPPHI L ,N ) 51 

0. 

L = 

7 

N= 

7 

YREFIL.N)= 

5.53047160E-01 

PYBPX(L,NI= 

0. 

PYBPPHI L.N) = 

0. 

L = 

B 

N= 

7 

YREFIL, Nl= 

5.530471 50E-01 

PYBPX(L,N)= 

0. 

PYBPPHI L *N I = 

0. 

L= 

9 

N = 

7 

YREFIL. N)= 

5. 53 047160 E— 01 

PYBPXIL,N>= 

0. 

PYBPPHI L ,NJ = 

0 . 

L = 

10 

N = 

7 

YREFIL, N)= 

5.53047160 E— 01 

PYBPX I L , N t = 

0 . 

PYBPPHI L »N) = 

0 . 


INITIAL SHOCK 

RS= 1 . 25 300000E+00 BS= 1 . OOOOOOOOE+OO XO= 1 .0 1 500000E+00 CX0= 1 . 07000000E+00 

GAS MODEL PARAMETERS 
I GA S= 1 

DENSO= 1 .66800000E-03 VELO= 6.93 500000E+03 


NQND IMENSIONAL FLIGHT CONDITIONS 

GAMMA= 1.40500000E+00 GAMM= 1 . 1 42C0000E+00 PINF = 1 . 58300000E-03 ALPHA= 2.20000000E+01 

PRINTING AND CUTOFF PARAMETERS 
KPO= 100 KPP= 100 KF= 300 

COMPUTAT ICN AL PARAMETERS 
I Y REF C= 1 


CT= 5.40000000E-01 EPSI= 3. 12500000E-02 PMIN= 1 . OOOOOOOOE-O 1 



T I ME= 2.07263396E+00 


UN 

VN 

WN 

RHO 

PHI 

RAD 

X AO 

GAMM 

DELtA(l.UN) 

PDPT(1,L,N» 

PDPXI l.L.N) 

PDPPHI 1 . L * N) 


PU.M.N) 

ENTROP 


300 1 

300 2 

300 3 

300 4 

300 5 

300 6 

300 7 

300 8 

300 9 

300 10 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


8.98355684E-02 

0. 

1 . 28655837E-01 

0. 

1.63678001E-01 

0. 

1 . 95 00744 9 E- 01 

0. 

2. 22413795 E-01 

0 . 

2.51330998E-01 

0. 

3. 2303798 8E-01 
0. 

3.97762628E-01 

0. 

4.55612542E-01 

0. 

5.06550282E-01 

0. 


0. 

0. 

0* 

9.63877017E-02 

0. 

1.92152523E-01 

0. 

2.86675609E-01 

0. 

3.79346129E-01 

0. 

4. 69565224E-01 
0. 

4.79690999E-01 

0. 

4.88201866E-01 

0. 

4. 946 3491 3 E— 01 
0. 

4.98640244E-01 


3.95391614E-12 
3. 55271368E-1 5 
6.12473820E-13 
-3.87734284E-03 
-1.78137211E-13 
-1 .548431 51E-0 2 
-1. 5106843 8E -13 
— 3.47459096E-02 
-3 .5921773 4E -14 
—6. 15376536E— 02 
1.43607998E-14 
-9. 56864122E-02 
'6. 6600611 2E -15 
-1.014695 77E-0 1 
3.43889315E-15 
-1.09440868E-01 
2. 12030507E-15 
-1.19166721E-01 
1.33189849E-15 
- 1. 301 18144E-01 


1.29678747E+01 
1.14504882E+00 
1.2278C700E+01 
1 . 1453002 8E+00 
1 . 1551 871 7E+0 1 
1 . 1455161 5E+00 
1.052 07 478E+01 
1 . 1455 1904E+00 
8 . 80338647E+0 0 
1.14512952E+00 
6.42701840E+00 
1 . 1445 2572 E *00 
4. 62810039E+00 
1.14531C72E+00 
3.30694317E+00 
1.14725795E+00 
2.51926588E+00 
1. 1 493102 7E *00 
1.9807 1323E *00 
1.15 14471 5 E+00 


8.24937531E-01 
4.38663956E-02 
7. 743 468 8 7E— 01 
4. 3808 1 123E-02 
7.220 50969E-01 
4. 37 8067 13E— 02 
6.50646190E— 01 
4. 39 10768 4E -02 
5. 3551546 9E-01 
4. 4363488 2E-02 
3. 79568658E-01 
4.51338766F-02 
2. 61977733 F-01 
4.53076413E-02 
1.77735949E— 0 1 
4.50670189E-02 
1.29450343E-01 
4. 47625553E-02 
9.76166154E-02 
4.44374734E-02 


300 1 2 

300 2 2 

300 3 2 

300 4 2 

300 5 2 

300 6 2 

300 7 2 

300 8 2 

300 9 2 

300 10 2 

tn 

-7 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1.63291730E-01 

0. 

2.08745398E-01 

0. 

2.52080056E-01 

0. 

2. 934891 8 9E- 01 
0. 

3.28256928E-01 

0. 

3 . 53006400E-01 
0. 

3.66601997E-01 

0. 

3. 80698641 E-01 

0. 

4..1121201 6E-01 
0. 

4. 9 8841 89 1 E-01 

0. 


-1 • 68684C79E-02 

0. 

— 1 . 61573847E-02 
9. 81940048E-02 
-1.49552171E-02 
1. 959389 80E- 01 
-1 .483131 82E-02 
2.92733201E-01 
-8. 46 80 10 14E— 03 
3. 8 8 20 0491 E-01 
2 . 17906278E-02 
4.82671237E-01 
7.48911851E-02 
5. 0181 89 57E— 01 
1 . 23667917E-01 
5.22427035E-01 
1.55791440 E-01 
5 .469 92 33 IE— 01 
l . 671571 02E-01 
5. 80 97965 0E -01 


8.88681512E-12 
2 . 13209011E-02 
1. 47718268E-1 2 
1 . 85 379660 E-0 2 
— 5.08548331E-13 
7. 85 7131 60 E-0 3 
-3. 75893845E-13 
- 1. 0 1235375E-02 
-9.35505276E-1V 
—3. 496467 16 E-0 2 
3.96000456E-15 
-6.48639670E-02 
1.33562120E-14 
-7. 14872332E-02 
7.33668589E-15 
-8.05402740E-02 
2 . 83848387E-15 
-9.27354575E-02 
4.70573406E-1S 
-1.10512970E-01 


1.30359501E+01 
1.14639301E*00 
1.23236163E+01 
l. 1472390 8E+00 
1.16317757E+01 
1.1 48 27952 E+00 
1.08494739E+01 
1.14946946E+00 
9. 8415960 0E+00 
1.1 50 48873 E+00 
8. 6255406 IE* 00 
1 .1 5099842E*00 
7.48716929E+00 
1.1 5 102210E+00 
6. 05542377E+00 
1 . 1 5092769E+00 
4.37354151E+00 
1. 15115555E+00 
2.58232957E+00 
1.15502816E+00 


8. 19450798E-01 
4.31649360E-02 
7.6377936 IE-01 
4. 281 84036E-02 
7.09485669E-01 
4.23920213E-02 
6.49661407E-01 
4. 19291839F— 02 
5.777 50563E-01 
4. 1613062 2E-0 2 
4. 969 2419 5E -01 
4. 16103683E-02 
4. 247 51453 E-01 
4. 18577554E-02 
3.36142378E— 01 
4.22991322F-02 
2 .341 87371F-01 
4. 28415045E-02 
1.26819499E-01 
4. 2 39 3 6 640 E— 02 


H(L.M.N) 

AMACH 


5.06218276E-01 
3. 430 16968E-01 
5 . 05392870E-01 
4.93901867E-01 
5.05441028F-01 
6.31836570E-01 
5.05847797E-01 
7.57571933E-01 
5.0471 1700E-01 
8.72344400E-01 
4. 992769 98E-D1 
1 .00278 93 9E *00 
4.98333326F-01 
1. 322 00 51 9E +00 
4.978 34466F-0 1 
1.67895888E+00 
4. 99318999E-01 
1 • 9747791 8E + 00 
5.02997762E-01 
2 .249 70543 E+00 


5.05733052F-01 
6.31955587F-01 
5.04821625E-01 
8. 13260991E-01 
5.04234143F-01 
9.90748490E-01 
5.03671 277E-0 l 
1 . 16613352F+00 
5.02712566F-01 
1.31866998E+00 
5.01687351E-31 
1.43609897E+00 
5.02377668E-01 
1. 53286847E* 00 
5. 034 2193 5E -01 
1 . 66040691 E+00 
5.04476481F-01 
1.86 2287 93 E+00 
5. 042 87 2 59F -01 
2.34073966E+00 



300 

1 

3 

1 

2.42881190E-01 

0. 

-3.46192830E-02 

0. 

1. 40 78862 0E -11 
4.26418022E-02 

1.33845359E+01 
1 • 14891828 E<-00 

8. 244 88422 E-01 
4. 1 86 0982 4E -02 

5.053451 55E-01 
9.57318564E-01 

300 

2 

3 

1 

2.96377721E-01 

0. 

— 3.17647460E-02 
1.00000308E-01 

6.57362581E-13 
4. 09532749 E -02 

1.27807897E+01 
1 . 15072641 E<-00 

7.70451920F-01 
4. 105 82455E-02 

5.04649774E-01 
1.17887 93 8E»-00 

300 

3 

3 

1 

3 • 52635529E-0 1 
0. 

-2.800929C1E-02 

1.99725437E-01 

-9.44353936E-13 
3.1 1985783E-02 

1 . 2 1832780E+0 1 
1. 15319431E+00 

7. 144 J.8134F-0 1 
3.99811630E-02 

5. 039 8426 6E -01 
1.42277015E4-00 

300 

4 

3 

1 

4.04362561E-01 

0. 

-2. 29185177E-02 
2 . 9879 0793 E-01 

—3.91 116129E-13 
1.44988346E-02 

1 • 16105282E«-01 
1.15603755E+00 

6.60078245E-01 

3.87781014E-02 

5.032151 18E-01 
1.65914819F1-00 

300 

5 

3 

1 

4.41964565E-01 

0. 

~7 • 30919040E-03 
3. 9705485 3E-01 

-6.16739097E-14 
— 8.39168966E-03 

1 . 1 l 166Q64E+01 
1.1 5837732 E+00 

6. 15847132E-01 
3.78305292E-02 

5.02882148E-01 

1.83822820E4-00 

300 

6 

3 

1 

4.56907608E-01 

0. 

3.98904015E-02 

4.95777250E-01 

9.30672944E-1 5 
-3.40415218E-02 

1. 06552 150E+01 
1.15924949E<-00 

5. 83006266E-01 
3.75380968E-02 

5.03477506E-01 
1.92 09760 5E +00 

300 

7 

3 

1 

4. 60430446E-0 1 
0. 

1. 26002738E-01 
5.23946915E-01 

9 . 30325496E-1 5 
-4.15048891E-02 

9.9670081 8E+00 
1.16019252EI-00 

5. 369 19801 E-01 
3. 727 193 91 E- 02 

5. 040871 73 E-01 
2.01721574Ef00 

300 

8 

3 

1 

4. 56056670E-01 
0. 

2 . 18869655E-01 
5 . 5 66 5 2205 E-01 

3.57257339E-15 
-5. 16 396801 E-0 2 

8.94380287E+00 

1.16193982E+00 

4. 6 88 0588 3 E-01 
3. 67605078E-02 

5.04043635E-01 
2.1692 7766F4-00 

300 

9 

3 

1 

4.452 82 968E-01 

0. 

3. 14613541E-01 
5. 99349748 E-01 

8.99417431E-1 6 
-6.63041935E-02 

7 . 71674880E+00 
1.16489047E+00 

3. 8767612 4F-01 
3. 58678401E-02 

5.03544202E-01 

2.39338679Et00 

300 

10 

3 

1 

4.26718473E-01 

0. 

4.25865829E-01 

6.63319057E-01 

-2.70165735E— 18 
-9.09077963E-02 

6.4608 801 7E*-00 
1 . 1704427 2 E+00 

3.0257401 5E— 0 1 
3.40746669E-02 

5.03321922E-01 

2.74823650E«-00 

300 

1 

4 

1 

3.229156UE-01 

0. 

6.39627034E-02 

-4.98713933E-02 

0. 

-8.78227822E-04 

0. 

6. 39627034E-02 
4. 2 46 5751 4E -02 

1 . 394335 76E+0 1 
1 • 15268 699E <-00 
0. 

8.33603585E-01 
3.998 1787 3E-0 2 

5.04718488E-01 
l • 29940414E+00 

300 

2 

4 

1 

3. 84058445E-01 
0. 

6. 746391 05E-02 

-4.02636100E-02 
1.0180661 IE-01 
-1.21 073483 E-03 

0. 

6.33685837E-02 
5. 40 176403 E-0 2 

1 .3673303 1E+0 1 
1 . 1558423 1E*00 
0. 

7.9253540 3E-01 
3.85592063E-02 

5. 0445 2 101 E-01 
1. 56450 351E+00 

300 

3 

4 

1 

4. 5 720 092 4E- 01 
0. 

7.09397174E-02 

-3 .462336 7 IE- 02 
2.03511895E-01 
—1.42 196 186 E-03 

0. 

5. 45400249E-02 
5.39158290E-02 

1 .32788132E*-0 1 
1 . 1606 845 2 E+00 
0. 

7.338 58700E-01 
3. 64736206E-02 

5.04318307E-01 

1.90971239E+00 

300 

4 

4 

1 

5. 13863636E-0 1 

0. 

7.60697091E-02 

-2.186 16584E-02 
3.04848384E-01 
- 1. 628343 04E-03 

0. 

3.91212067E-02 
6. 58605292E-02 

1.29219890E+01 
1 • f6531 899E+00 
0. 

6.818 17127E-01 
3.45632597E-02 

5.0 419 638 6E- 01 
2.196 70400E4-00 

300 

5 

4 

1 

5 . 59035780E-01 
0. 

8.40280184E-02 

-2. 3 24 2 1368 E-03 
4.05909215E-01 
— 1.8772222 IE-03 

0. 

1.81812923E-02 
8. 633460 63E-02 

1 . 2601 1923E+01 
1.16970725EI-00 
0. 

6. 3 5 835 08 8E -01 
3. 28235099E-02 

5.04047237E-01 
2. 44 58 825 OEt- 00 

300 

6 

4 

1 

5. 57939571 E-01 
0. 

1.00479 43 3 E-01 

4.22860303E-02 

5.08883264E-01 

-1.92829294E-03 

0. 

-3.21907660E-33 

4.98406533E-02 

1.2 5966857E+01 
1 • 16977046E<-00 
0. 

6. 351 94793E-0 1 
3.27989259E-02 

5.03989816E-01 
2. 44 89131 7E <-00 

300 

7 

4 

1 

5.46606067E-01 

0. 

1.11 791 267E-01 

1.67392237E-0-1 

5.46074873E-01 

-1.99875467E-03 

0. 

-1. 15225450E-02 
1.18570333E-01 

1.250 26 563E<01 
1.171 09 953E+00 
0. 

6 . 21 87053 9E-0 1 
3.22850241F-02 

5.03840822E-01 
2 . 52021 647E<-00 

300 

8 

4 

l 

5 . 17 3 34470 E- 01 
•0. 

1.3 43 85486 E-01 

2. 88754003E-01 
5.90877374E-01 
-2. 24 590 178 E-03 

0. 

-2.27390862E-02 

2.20242822E-01 

1 . 23337573E+01 
1.17353643E+00 
0. 

5.981 05469E-0 1 
3. 135722 14E-02 

5.03442302E-01 
2. 646 8 1 242 E<- 0 0 

300 

9 

4 

1 

4.63788848E-01 

0. 

1.75952260E-01 

4. 050147 14E-01 
6. 51707 165 E-01 
-2. 19329304E-03 

0. 

-3. 987292 96E-D2 
4.08658497E-01 

1.2138251 8E+01 
1 . 1764390 5E+ 00 
0. 

5. 70855633E-0 1 
3. 02 7 48 96 9E- 02 

5. 031459 94 E-01 
2.7942 95 13EI-00 

300 

10 

4 

1 

3.98869770E-01 

0. 

2. 539 23740 E-01 

5. 16729123E-01 
7. 45658463E-0 1 
-7.97417623E-04 

0. 

-7.13026227E-02 
8. 10005915E-01 

1. 18332569E+01 
1 • 1 81 14942 E»00 

0 a 

5. 28823916E-01 
2.85637071E-02 

5.04443066E-01 

3.03742088E<00 



300 1 1 

300 2 1 

300 3 1 

300 4 l 

300 5 1 

300 6 l 

300 7 1 

300 8 l 

300 9 1 

300 10 1 

300 1 2 

300 2 2 

300 3 2 

300 4 2 

300 5 2 

300 6 2 

300 7 2 

300 8 2 

300 9 2 

300 10 2 


2 7.77998844E-02 

3.00000000E+01 
2 1.17679417E-01 

3. 000000 00 E +01 
2 1 . 53 16522 9E-01 

3. 000 000 00E +01 
2 1.84530709E-01 

3 . OOOCOOOOE+O 1 
2 2. 1 2 560 55 7E- 01 

3 . 000000 00 E+0 1 
2 2 .4572466 IE-01 

3.00000000E+01 
2 3.14637729E-01 

3.00000000E+01 
2 3.879U897E-01 

3 .000 000 00 E+0 1 

2 4.44308659E-01 

3. 000 000 OOE +01 
2 4.92084133E-01 

3.00000000E+01 


2 1. 4141478 6E-01 

3.00000000E+01 
2 1.889 13748E- 01 

3. 000 000 00 E+0 l 
2 2. 31834223E-01 

3.00000000E+01 
2 2.725 66 634E-01 

3.00000000E+01 
2 3.07772628F-01 

3 . 000 000 00 E+ 01 
2 3.34856800F-01 

3 . OOOOOOOOE+O 1 
2 3.52096303E-01 

3. OOOOOOOOE+O 1 
2 3. 6851 83 19 E- 01 

3 . OOOOOOOOE+O 1 
2 3.99268258E-01 

3.00000000E+01 
2 4.86059497E-01 

3.00000000E+01 


0. 

0. 

0. 

9. 6387701 7E-02 

0. 

1. 92152523E-01 

0. 

2 . 866 756 09 E- 01 
0. 

3. 793461 29E-01 
0. 

4.69565224E-01 

0. 

4. 79690999 E-01 
0. 

4.88201866E-01 

0. 

4.94634913E-01 

0. 

4. 9 8640244 E-0 1 


-1.68684C79E-02 

0. 

— 1.59 187 856 E- 02 
9. 8 2 03800 6E -02 
—1.47431 197 E-0 2 
1.95938714E-01 
— 1.48755 745 E-02 
2.926S3738E-01 
-9.85245599E-03 
3. 880549 69E-01 
1.68886951E-02 
4. 82254905 E-01 
6.64301950E-02 
5 .0093972 1 E-0 1 
1.1284 80 61E-01 
5.20664989E-01 
1 .43 464930 E-01 
5 .4341 4433E-01 
1 . 53063041E-0 1 
5.73562880E-01 


— 4.491 77842E— 02 
3.55271368E-15 
-4..15779783E-02 
-3.87734284E-03 
— 3.96747878E-02 
-1.54 8431 5 IE -02 
-4.10135007E-02 
-3.47459096E-02 
-4.46546353E-02 
-6. 153765 36 E-02 
-4.95465570E-02 
-9.56864122E-02 
-5.26956090E-02 
-1. 014695 77E-01 
—5. 5 2 4669 77 E-0 2 
- 1.09 4408 68E -01 
-5.87262835E-02 
-1.19166721E-01 
-6.3 562351 8E -02 
-1. 301181 44E -01 


-8. 16458648E-02 
2.13209011E-02 
-7.71887767E-02 
1.86595269E-02 
-7 .86 4 79 940 E-0 2 
7. 85549043E-03 
-8.24528730E-02 
-1 • 02 8 3 9447 E-0 2 
-8. 6598332 8E -02 
-3.54014012E-02 
-8. 9768245 6E -02 
-6.58430892E-02 
-8, 91 6 00945 E-0 2 
-7.26785577E-02 
-8.76974529E-02 
— 8.20281899E-02 
-8. 707460 49 E-02 
-9.45416647E-02 
-8. 8903644 8E-02 
— 1.12278917E-01 


1.29578747E+01 
1.14504882E+00 
1 . 2355 1229E+0 1 
1.14527414E+00 
1.17091759E+01 
1 . 14547 548E +00 
1.07568286E+01 
1 • 14547856E +00 
9. 09 l 9 5404 E+0 0 
l • 145 12419E+00 
6.78507511 E+00 
1.1445S986E+00 
4.87622807E+00 
l. 14524967E+00 
3.49537161E+0C 
1.14701 90 1E+00 
2.65573964E+00 
1 . 148921 80E+00 
2.09340263E+00 
1. 15079333E+00 


l . 30359501 E+0 1 
1. 1463 9 301 E+00 
1.24150862E+01 
1.1 47 14942 E+00 
1 . I 8067609E+0 1 
1.14810116E+00 
1.10851746E+01 
1 • 14919597E+00 
1 .009 7947 7E+01 
1 . 15016730E+00 
8. 84578801 E+00 
1 . 15072359E+00 
7 • 63 18.8068 E+00 
1.15080376E+00 
6.15215651E+00 
1. 15070464F+00 
4.45227560E+00 
1. 15082197E+00 
2.65191315E+00 
1 .15450128E+00 


8.24937531E-01 
4. 38663956E-02 
7.79954617E-01 
4. 38 132 142 E-0 2 
7.33286243E-01 
4.37827285E-02 
6. 6716512 IE -01 
4.38996586E-02 
5.55098 42 IE-01 
4.4318841 8E -02 
4.028 54419E-01 
4. 501 432 94E -02 
2.777 2001 3E-01 
4. 52458957E— 02 
1. 89328469E-01 
4.50628562E-02 
1.37637987E-01 
4.48107640E-02 
1.042 99 182 E-01 
4.45703129E-02 


8. 19450798E-01 
4.31649360E-02 
7.70723844E-01 
4. 28523844E— 02 
7 . 2252971 7E-01 
4.24562013E-02 
6. 6699183 4E-01 
4.20248914E-02 
5. 961 74304E-01 
4.17193145E— 02 
5. 12246910E-01 
4. 1691478 IE-02 
4.346 81 31 OE— 01 
4. 192 1 3480E-02 
3.42747266E-01 
4. 23678613E-02 
2.39525122E-01 
4.2948560 1E-0 2 
1.31167459E-01 
4. 2542827 7E -02 


5. 0621827 6E -01 
3.43016968E-01 
5. 05460 80 2E -01 
4.78 8462 99E- 01 
5. 056261 40E-01 
6.10050592E-01 
5.06223431E-01 
7. 33082704E-01 
5 . 05342236E-01 
8. 500404 94E-01 
5.01397638E-01 
9.97054856E-01 
4. 999508C5E-01 
1 . 300 8 51 33E + 0 0 
4. 99354387E-01 
1 .646 3473 OE+OO 
5.00268355E-01 
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1.15508411E+00 
-9. 23463 54 9E-04 
1 . 350 31 83 1 E+01 
1.15739554E+00 
-4.36761622E-03 
1 . 32763 177E+0 1 
1 . 16071607E+00 
-8 .2830890 IE— 03 
1. 3025990 3E+01 
1.16394256E+00 
— 1 . 6 76 2 39 76 E— 02 
1 . 2 9470 95 3 E+01 
1.1649846 IE +00 
-3. 33327 07 5 E— 02 
1.27410545E+01 
1.1 6776728E +00 
-4. 17595932E-02 
1.24533775E+01 
1.17180390E+00 
-5. 50295 4C 7 E-02 
1 . 22049520E+0 1 
1.17543872E+00 
-7 . 7 1075699E-02 
1.1 9648035E+0 1 
1.179090036*00 
-1. 16784685 E- 01 


8.24488422E-01 
4. 1 86 0982 4F -02 
7.954029146-01 
4.136967496-02 
7. 620349 82E-01 
4.06593029E-02 
7.227017806-01 
3.983344126-02 
6.797276116-01 
3.90653735E-02 
6.31589706E-01 
3.860726176-02 
5.67494509E-01 
3. 8152493 7E-02 
4.883699056-01 
3.747507106-02 
4.022920926-01 
3.65270376E-02 
3. 14404891E-01 
3.492981866-02 


8. 33603585E-01 
3.998 17873E -02 

8.021858116-01 
3. 88960903E-02 

7.67050178E-01 

3.766142896-02 

7.33492142E-01 

3.646038096-02 

6.970114206-01 

3.512887296-02 

6.854491496-01 

3.469856996-02 

6. 55784925E-01 
3.358449036-02 


6.149145426-01 

3.201496356-02 

5. 80 1 2 2492 E— 01 
3.06452470E-02 

5.468903366-01 

2.93057291E-02 


5.053451556-01 
9.573185646-01 
5.050067596-01 
1.092912096+00 
5.047763336-01 
l . 26273044E + 00 
5.044379996-01 
1.441959436+00 
5.042805786-01 
1.606020676+00 
5.043910036-01 
1.727832986+00 
5.044874846-01 
1.867673076+00 
5. 042 7691 6E -01 
2.056835006+00 
5.038575906-01 
2.294575066+00 
5.03277573E-01 
2.628766856+00 


5.047184886-01 

1.29940414E+00 

5.04562220E-01 

1.505221446+00 

5.04566819E-01 

1.719848806+00 

5.045 1 2410E-0 1 
1.91355 64 2E+00 

5.045472166-01 

2.117999496+00 

5.64625334E-01 
2. 18040346E+00 

5.046118606-01 
2 . 341 96267E + 00 

5. 0 460 800 7E -01 
2.56395811E+00 

5.04678744E-01 

2.754594306+30 

5.047742486-01 

2.938925436+00 



300 1 1 

300 2 1 

300 3 • 1 

300 4 1 

300 5 1 

300 6 1 

300 7 1 

300 8 1 

300 9 1 

300 10 1 

300 1 2 

300 2 2 

300 3 2 

300 4 2 

300 5 2 

300 6 2 

300 7 2 

300 8 2 

300 9 2 

300 10 2 


4 -1.84246256E-11 
9.00000000E+01 
4 4.37300957E-02 

9.00000000E+01 
4 8. 3072198 4E-02 

9. OOOOOOOOE+O 1 
4 1.18463001E-01 

9.00000000E+01 
4 1. 51473700E-0 1 

9. OOOOOOOOE+O 1 
4 1.90402504E-01 

9. 00000000E+0 1 
4 2 . 6016611 2E-01 

9. 000000 00 E +01 
4 3.29711382E-01 

9.00000000E+01 
4 3.82662140E-01 

9.00000000E+01 
4 4. 241 50 766 E- 01 

9. OOOOOOOOE+O 1 


4 -3.34899532E-11 
9.00000000E+01 
4 5 . 56283501 E-02 

9.00000000E+01 
4 1 . 0 1121 172E-01 

9. OOOOOOOOE+O 1 
4 1.45523363E-01 

9.00000000E+01 
4 1 . 86731390E-01 

9.00000000E+01 
4 2 . 2 5862857E-01 

9. 0000 00 00 E +01 
4 2.63644115E-01 

9. OOOOOOOOE+O 1 
4 2.97917736F-01 

9.00000000E+01 
4 3.34063327E-01 

9.00000000E+01 
4 4.00516776E-01 

9.00000000E+01 


0. 

0 . 

0 . 

9.63877017E-02 

0 . 

1 . 92152523E-01 

0. 

2.86675609E-01 

0 . 

3. 79 3461 2 9E- 01 
0. 

4.69565224E-01 

0. 

4. 7969C999E-01 
0. 

4. 88201 86 6E— 01 
0. 

4. 9463491 3E-01 
0. 

4. 98640244 E-01 


-1.68684079E-02 

0. 

-1 . 5756 1481 E-02 
9.82002102E-02 
-1.42 7471 64E-02 
1.95776444E-01 
-1.502912 75E-02 
2.92 177 7 82 E-01 
-1.62566185E-02 
3. 86615881 E-01 
— 8. 6746 514 4E -03 
4.78636503E-01 
1 . 96 8 803 09 E-02 
4. 93 9C 984 76-01. 
5.07505296E-02 
5 .0828221 5E-01 
7. 4165 48 16 E-02 
5.22 121 706E-01 
8.81503393E-02 
5. 36355892 E-01 


— 8.98355684E-02 
3.55271368E-15 
-8.47405574E-02 
-3.87734284E-03 
-7.88015434E-02 
-1.548 43 151E-02 
-7.75804327E-02 
-3. 474 590 96 E-02 
-8.04791942E— 02 
-6. I 53 7653 6E -02 
-8.54921075E-02 
-9. 56864122E-02 
-8. 852816 46 E-02 
-1.01469577E-01 
-9.0 145248 6E-02 
-1.09440868E-01 
-9.24392411E-02 
-1.19166721E-01 
-9.61704581E-02 
-1.301181 44E-01 


-1.63291730E-31 
2.13209011E-02 
-1.607 84 116E-01 
1.86149711E-02 
-1.57834413E-01 
6.85518481E-03 
-1.57601364E-01 
-1.238U574E-02 
-1.59643654E-01 
-3. 972 0272 9E -02 
-1.63192273E-01 
-7.43527698E-02 
-1. 621 40 58 7E -01 
-8. 2203705 8E -02 
-1 . 58937606E-01 
-9.24845164E-02 
-1.54789192E-01 
-1.05290738E-01 
-1.48290696E-01 
-1.21137974E-01 


1. 29678747 E+01 
1 • 14504882E+00 
1.2 90 57919 E+01 
1 . 145 0766 4E +00 
1.28285098 E+01 
1. 14513530 E+ 00 
1 . 25524591E+01 
1.14507 67 1E+00 
1.1 5479727E +01 
1. 14474406E+00 
9.45948848E+00 
1.1442C796E+00 
7.03144192E+00 
1 . 14440 453 E+00 
5. 1 5000641 E+ 00 
1.14532311E+00 
3. 92652613 E+00 
1. 14642 235E+00 
3.08665329E+00 
1. 1472915 9E +00 


1.30359501E+01 
1.14639301E+00 
1.30092793E+01 
1.14652869E+00 
1.29576826E+01 
1.14683768E+00 
1. 27376 03 2E+01 
1.14729033E+00 
1. 2 09 6295 4E +01 
1 . 14773877E+00 
1.07691177E+01 
1. 14809786E+00 
9.06751166E+00 
1.14827328E+00 
7.26225088E+00 
1.14827016E+00 
5. 43987639 E+00 
1 .14808955E+00 
3. 57248283 E+00 
1. 14902307E+00 


8. 2493753 IE -01 
4. 3 86 6395 6E— 02 
8. 2 03 2047 IE -01 
4.385 81 1 76E-02 
8.14397565E-01 
4.38353792E-02 
7.9531497 6E-0 1 
4.38945279E-02 
7 . 26940993 E-0 1 
4. 41 7758966-02 
5 .849 52946E-0 1 
4. 47 2 2430 7F -02 
4.19852572E-01 
4. 5 05 42 09 3E— 02 
2. 94859243E-01 
4.51195600E-02 
2. 16373133E-01 
4. 5 104380 6E -02 
1.646 6 62 5 2 E-01 
4. 51876035E— 02 


8. 19450798E-01 
4. 3 1 64936 OE -02 
8. 16557577E-01 
4, 309 8634 7E— 02 
8. 10634582E-01 
4.29473982E-02 
7.920144366-01 
4.274414576-02 
7. 44 7 1 26 7 8 E— 01 
4. 2597746 9E-02 
6.523302096-01 
4.260138296-02 
5.38081042E-01 
4. 2 79 466 79F— 02 
4. 2 0663 00 5E -01 
4.317097926-02 
3.05952788E-01 
4.376554446-02 
1. 907255316-01 
4.416043796-02 


5. 0621 82 76E -01 
3.430 1 6968E-01 
5. 06 2372 40 E-01 
3.64292998E-01 
5.074472296-01 
4.37 772 624E-01 
5 • 1011 5015E-01 
5. 420118066-01 
5. 12563515E-01 
6.58927497E-01 
5. 12427863E-01 
8.10082270E-01 
5.109688246-01 
1.08929474E+00 
5.09650315F-01 
1. 3895627 OE +00 
5.089392986-01 
1.63729711E+00 
5.101 16 104E-0 1 
1. 84371 235E+ 00 


5.0573 30 52E-01 
6.31955587E-01 
5.057264526-01 
6.584053726-01 
5.06230352F-01 
7. 25988393E-01 
5. 074540106-01 
8.33692597E-01 
5.085924796-01 
9.604169536-01 
5. 084490376-01 
1.098134006+00 
5.07652390E-01 
1.236521256+00 
5.068898076-01 
1.3806 807 IE +00 
5.O6559627E-01 
1.544724156+00 
5.067236896-01 
1.84893 93 6E+00 



300 1 3 

300 2 3 

300 3 3 

300 4 3 

300 5 3 

300 6 3 

300 7 3 

300 8 3 

300 9 3 

300 10 3 

300 1 4 

300 2 4 

300 3 4 

300 4 4 

300 5 4 

300 6 4 

300 7 4 

300 8 4 

300 9 4 

300 10 4 


4 -4.98131762E-11 
9.00000030E+01 
4 5.35315950E-02 

9 .OOOOOOOOE +0 1 
4 1.10725829E-01 

9 . OOOOOOOOE + 01 
4 1. 6582300 5E-01 

9.00000000E+01 
4 2.21611 49 8E-01 

9.00000000E+01 
4 2.72167337E-01 

9.00000000E+31 
4 3. 2C63044 1E-0 1 

9.00000000E+31 
4 3.64 1 91 640E-0 1 

9.00000000E+01 
4 3.99377211E-01 

9. OOOOOOOOE *-01 
4 4.27656652E-01 

9. OOOOOOOOE +01 


4 -6. 6227 657 2E- 11 
9.00000000E+01 
6. 39627034E-02 
4 7.21629490E-03 

9. 00000 000 E+0 1 
6.7 695674 7E-0 2 
4 1 . 1202272 8E-01 

9. OOOOOOOOE +01 
6.78945839E-02 
4 1 . 771 1 7353E-01 

9 .OOOOOOOOE 30 1 
6. 90948942 E-02 
4 2.60310525E-01 

9.00000000E301 
6.89900466E-02 
4 3.25224377E-01 

9.00000000E301 
6 . 95464708E-0 2 
4 3.94261330E-01 

9.00000000E301 
7.18341491E-02 
4 4.45367951E-01 

9. OOOOOOOOE 301 
7.38456994E-02 
4 4.62384847E-01 

9.000000006301 
9.23720766E-02 
4 4.41884518E-01 

9.00000000E+01 
1 . 16310023E-01 


-3.46192830E-02 

0 . 

—3 . 3 79490 C 3 E-02 
1 .000 127 19E-0 1 
-3.19808419E-02 
1.99400365E-01 
-3. 12791 32 9 E-02 
2. 97679954 E- 01 
-3. 009 1839 6E -02 
3. 93885633 E-01 
— 1.66 1159 32E-02 
4.87707782E-01 
2.88428703E-02 
5. 0812869 5E-01 
8. 9 14262 51 E-02 
5 . 2836 2563E-0 1 
1.591 26 101 E-01 
5.496C8499E-01 
2. 45705389 E-01 
5.74C71539E-01 


-4.98713933E-02 

0 . 

-8.78227922E-04 
-5.07057001E-02 
1 .0 1825227E-01 
-9. 37224323E-04 
-5 .096 5 9044 E-02 
2. 03024286 E-01 
-7. 8063 392 6E -04 
— 4.69865755E-02 
3.03182127E-01 
-7.50693349E-04 
-4.25064279E-02 
4.01155385E-01 
-3. 80022812 E- 04 
-3. 0264878 7E-02 
4.96779061E-01 
-L.68568368E-04 
3.551 47 845E-02 
5.22347543E-01 
-1.2877 735 5E -04 
1. 20782043 E- 01 
5.48442912E-01 
-4.39143028E-05 
2 . 22385074E-0 1 
5 .77095292E-0 1 
1. 33228100 E-04 
3.30888783E-01 
,6.1 17871 86E-01 
6.024686C8E-04 


-2. 428811 90E-01 
4.26418022E-02 
-2.45841 47 IE-01 
4.11072851E-02 
-2. 48 7645 05 E-01 
2 . 9 1 946847E-02 
-2.48370586E-01 
9.98359483E-03 
-2.46462231E-01 
-1. 790 2892 3E -02 
-2.43612519E-01 
-5.30191274F-02 
-2.43202772E-01 
-6.29378343E-02 
-2.43991167E-D1 
-7. 552 81649E -02 
-2 • 46 917 85 6E -01 
-9. 141 47 5 35E— 02 
-2. 53128058E-01 
-1.12157805E-01 


-3. 229156 11E-01 
6. 3962703 4E -02 
2.57108601E-03 
-3.26871272E-01 
6.35995991E-02 
5.56783130E-02 
-3.41992318E-01 
5. 15341846E-02 
3.05067394E-02 
-3.44290207E-01 
3.23483471E-02 
3.76538413E-02 
-3.37172624E-01 
3.91448839E-03 
2 .83745089E-02 
-3.16468081E-01 
-3.16854849E-02 
9.05656448E-03 
-3.12274645E-01 
-4.36719629E-02 
4.041 16068E -02 
-3.08414051E-01 
-5. 85718134E-02 
7.81709142E-02 
-3.02847674E-01 
-7.75387696E-02 
1 . 347221 51E-01 
-2 . 9 6 4090 5 6E -31 
-1.03177635E-01 
2.32568767E-01 


1.338 45 359E+01 
1. 14891828E+00 
1. 33708256 E+ 01 
1.1492C585E+00 
1. 3313068 8E+01 
1.14990011E+00 
1.31089 84 4E+01 
1.15084734E+00 
1.27235442E+01 
1. 152100826*00 
1. 19856653 E+01 
1. 1534023 9E+00 
1 . 0852 5778E *-0 1 
1 . 15503582E+00 
9.58103129E+00 
1. 157291C7E+00 
8. 232 83 28 8E +00 
1.16026515E+00 
6. 8201 1 188E+00 
1.16485904E+00 


1.39433576E+01 
1. 15268699E+00 

0 . 

1 . 39258056E+01 
1.15288829E+00 
-5. 562 0062 2 E-0 3 
1.37855960E+01 
1.15451507E+00 
-7. 65693410 E-03 
1.36593279E+01 
1. 1560090 4E+00 
-l. 07654 196E— 02 
X . 34643284E +01 
1. 15837182 E+ 00 
—1.7 927 2947 E-0 2 
1.33142761 E+01 
1. L6023743E+00 
-3.50829114E— 02 
1.300 55162 E+01 
1.16421 149E+00 
-4. 16496 839 E-02 
1.26498895E+01 
1. 16902701E+00 
-5. 0481 2980 E-02 
1 .2336 2602 E+01 
1. 17349985 E+00 
-6.46718203E— 02 
1.20859609 E+01 
1 . 1772306 5E+00 
-8 .7666 3 07 5 E-02 


8.24488422E-01 
4. 1 8609824F-02 
8.213 83584E-01 
4.172 13543E-02 
8. 1 22 6 89 8 2 E— 01 
4.13896995E— 02 
7.91581267F-01 
4.09585093F-02 
7.57109782E-01 
4.04146624E-02 
7.00263247E-01 
3.99144210E-02 
6. 1812616 8E -01 
3.93549650E-02 
5. 2 76 9303 7E- 01 
3.86012813E-02 
4. 345 0608 5E -01 
3. 764 57466E— 02 
3.38458834F— 01 
3.61623203E-02 


8.33603585E-01 

3.99817873E-02 

8.30909899E-01 
3. 988934 14E-02 

8.09516891E-01 
3.915091 85E-02 

7. 90431065E-01 
3.84855266E-02 

7. 6126865 8E -01 
3. 74560703E-02 

7. 39073 73 IE -01 
3. 666 16874E-02 

6. 9391 92 02 E-01 
3. 501 29550E-02 

6.4276404 8E-01 
3.308 89069E-0 2 

5.98456078E-01 
3. 13710157E-02 

5.6361121 8E-01 
2.99836876E-02 


5. 053451 5 5E -01 
9.57318564E-01 
5.05373885E-01 
9.92331410E-31 
5.05620034E-01 
1 .076349 52E +00 
5.05769400 E-01 
l« 1864565 6E +00 
5.061033 66 E-31 
1.326854 72 E+00 
5 .06 135580E-01 
1.47381800E+00 
5.05725491E-01 
1.650 7 6562 E+00 
5.05292469E-01 
1.86629399E+00 
5.04985004E-01 
2. 11951046E+00 
5. 04318485 E-01 
2.45031650E+00 


5. 0471848 8E -01 
1. 29940 414E+00 

5. 04665695 E-01 
1.31 73 5651 E+00 

5.04814816E-01 
1 . 46 1 17779E + 00 

5.04848496E-01 
1.58 164857E + 00 

5. 051723 90E-01 
1.759 73649E + 0 0 

5.0535 1933E-0 1 
1 . 88946 1 64E + 0 0 

5.05386333E-01 

2.14036483E+00 

5.05455576E-01 
2. 4193133 9E +00 

5.05605781E-01 

2.66116365E+00 

5.06060851E-01 

2.85561111E+00 



300 1 

300 2 

300 3 

300 4 

300 5 

300 6 

300 7 

300 8 

300 9 

300 10 


1 

1 

1 

l 

1 

1 

1 

1 

1 

1 


5 -4.49177842E-02 
1 . 20000000 Ei- 02 
5 -1.88978427E-03 
1.20000000E+02 
5 4.06434146E-02 

1.20000000E+02 
5 8.062 48 560E-02 

1 . 2 0000000 E* 02 
5 l. 17683821E-01 
1 . 2000003 0E + 02 
5 1 . 5 9949 99 2 E- 01 

1.20000000E+02 
5 2.31361437E-01 

1 .20000000E+02 
5 3.00783172E-01 

1 . 20000000E+02 
5 3.55512403E-01 

1.20000000E+02 
5 4. 033 66 171 E- 01 

1 . 20000000E+02 


0. 

0. 

0. 

9.63877017E-02 

0. 

1 .9215252 3E-0 1 

0. 

2 . 86675609E-01 

0. 

3 . 79346129E-01 

0. 

4.69565224E-01 

0. 

4. 79690999E-O1 
0. 

4. 88201866E-01 
0. 

4.94634913E-01 

0. 

4.98640244E-01 


-7. 77998844E-02 
3.55271 36 8E-1 5 
-7. 3330358 5E -02 
-3 . 87734284E-03 
-6. 6933577 5E -02 
-1.548431 51 E-02 
-6.45489690E-02 
-3.47459096E-02 
-6.50207332E-02 
-6. 1537653 6E -02 
-6.73901355E-02 
-9. 56 8 64 12 2E -02 
-6. 8725686 IE-02 
-1. 014695 77E-01 
-6. 9195314 8E -02 
-1 . 09 440 86 8 E— 01 
-6.97089898E-02 
— 1.19166721E-01 
-7.13973761E-02 
- 1 .301 18144E-01 


1.29678747E+01 
1.14504882E+00 
1. 32535697E+01 
1.14495658E*00 
1.35465388E+01 
1.14491127E+00 
1 . 37758094E+0 1 
1.14479243E+00 
1.33360211E+01 
1. 14443906E+00 
1 . 1562 1252E+01 
1 . 14387620E+00 
8 . 8 1317C65E+00 
1. 14391996E+00 
6.62076763E+00 
1 .1445764 1 E+00 
5 . 0953262 1 E+00 
1.1 4539375E+C0 
3 .91696663 E+00 
1 . 14620849E+00 


8. 2493753 1E-0 1 
4.38663956E-02 
8. 4596801 IE -01 
4.38865604E-02 
8.67229629E-01 
4.38824719E-02 
8.846 18714E— 01 
4. 39241599E-02 
8.56007372E-01 
4. 4152315 8E -02 
7. 33833258E-01 
4.46286657E-02 
5.42050491E-01 
4.4965992 0E-02 
3.921 18153E-01 
4.50634689E^02 
2.912 16932E— 01 
4.51051618E-02 
2.16178274E-01 
4. 5 20 306 31 E-02 
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3. 62 966 843 E-02 
6 4.77534164E-01 

1. 500000 DOE+32 
3.92413972E-02 
6 5.25978555E-01 

1 . 5 0000000 E+0 2 
4. 577933696-02 


-4. 9871393 3E -02 

0 . 

-8 . 78227822E-04 
-5. 71567339E-02 
1.01476785E-01 
-5 . 965001 77E-04 
-6. 1 2 87 5 132 E-02 
2. 0 2 3847 2 8 E— 01 
-2.561C5012E-04 
-6. 1 8 12049 8E— 02 
3.00100916E-01 
-1 . 5637 172 8E-04 
-7. 14469624E-02 
3.96638170E-01 
5. 75173320 E-04 
— 1.10764142E-01 
4.84842855E-01 
1 . 6305 304 IE- 02 
-6. 71540340E-02 
5.01212703E-01 
2.28561419E-03 
-1.47949067E-02 
5. 15933882E-01 
9.65331234E-04 
6. 332 73 791 E-02 
5. 2 966 56 3 5 E-01 
1 .2773436 8E-03 
1.63505982E-01 
5.43174600E-01 
1. 578130C6E-03 


-1.61457805E-01 
6.39627034E-02 
1.62443725E-03 
— 1.65654895F-01 
5.92756108E-02 
-2.51362711E-02 
— 1.76508697E-01 
4.75916582E-02 
-2 . 73049806E-02 
-1.61715600 E-01 
1.98241 105E -02 
-3.25665970E-02 
- 1.6769 86 39E-31 
-9.64220485E-03 
— 6.82475479E-02 
- 1 . 570 82549E-01 
-5.97568O21E-02 
-1.70005036E-D2 
- 1.68011 649E -01 
-7. 230 867 97E -02 
-2 . 6 2 434 8 9 5E -03 
-1.67429661E-01 
-8.60232558E-02 
1. 7187370 8F-0 2 
-1. 66660 833E-D1 
- 1.O1482384F-01 
3.91757360E-02 
- 1. 65428803E -01 
— 1 . 19514426E-01 
6. 8 8 826 942 E-02 


1.39433576E+01 
1. 1526869 9E+00 
0 . 

1.41238263E*01 
1.15064667E400 
-2 .4 593 697 IE— 03 
1.42760779E+01 
1.14896627E»00 
-2. 41059951 E-03 
1.44046463E+01 
1.14757556E+00 
- 8. 3 80 2 749 5E— 03 
1 .43 7 79 596 E*0 l 
1.14786222E+00 
-8 .5 3533475E-03 
1.41222451E+01 
1 . 1 50664 13E+00 
-1. 0524348 7E-02 
1. 38597849 E+ 01 
1.15364984E*00 
-1.1098 40 82 E-02 
1.34510809E+01 
1.15853486E+00 
-1. 23327 84 2E-02 
1 .2983921 1E*-01 
1. 16449648E+00 
-1 .41850246E— 02 
1.25157660E+01 
1.17091308E+00 
-1.71241297E— 02 


8.33603585E-01 
3. 998 17 87 3 E-02 

8. 615 17761E-01 
4.09331693E-02 

8 .8540652 8E-01 
4.173 8421 8E-02 

9. 05 8 56999 E-01 
4.24219251E-02 

9.015 886 89E— 01 
4. 2 2 79662 IE -02 

8.61272868E-01 

4.09249139E-02 

8. 20811605E-01 
3.95417877E-02 

7. 59300639E-01 
3.73860188E-02 

6.90785179E-01 
3. 489681 5 7E— 02 

6.237 24046E-0 1 
3. 23567746E-02 


5.04718488E-01 

1.29940414E+30 

5.04957952E-31 

1.09776689E«-00 

5. 052 592 76E -01 
9.01514960E-01 

5. 053479 92E-01 
6.961 15404F-31 

5.06020538E-31 

7.56764058E-01 

5. 20472507E-3 l 
1. 299 26 41 4F 4- 30 

5.07485030E-31 
1. 417770206 ♦■OO 

5.063061 5 OE- 01 
1.782043336+30 

5.06545838E-01 
2 . 1672331 3E+30 

5.06793595F-01 
2. 53284421 E+ 00 



300 1 1 

300 2 1 

300 3 1 

300 4 1 

300 5 1 

300 6 1 

300 7 1 

300 8 1 

300 9 1 

300 10 1 

300 1 2 

300 2 2 

300 3 2 

300 4 2 

300 5 2 

3 00 6 2 

300 7 2 

300 8 2 

300 9 2 

300 10 2 


7 -8.98355684E-02 
1.80000000E+02 
7 -4.70839953E-02 
1.80000000E+02 
7 -1.56504112E-03 
1.80000000E+02 
7 4. 257401 91 E- 02 

1 . 8000000 0E *02 
7 8. 43 59456 9E- 02 

1.80000000E+02 
7 1.31016685E-01 

1 • 80000000E+02 
7 2.04918638E-01 

1 . 80000000E + 02 
7 2.75212503E-O1 

1 . 80000000E+02 
7 3.34076879E-01 

1.80000000E+02 
7 3 • 93089420E-0 1 

1.800000 00 E+02 


7 -1.63 291730 E- 01 
l . 80000000E *-02 
7 -1.11919040E-01 
1.80000000E+02 
7 -5.21777061E-02 
1 . 80000000E+02 
7 1.01257782E-02 

1 . 80000000E+02 
7 6.93485890E-02 

1 . 80000D00E+02 
7 1.27267548E-01 

1 .80000030E + 02 
7 1.91961870E-01 

1 .80000000 E+02 
7 2.51977632E-01 

1.80000000E+02 
7 3.03793830E-01 

1 . 80000000E+02 
7 3.60625008E-01 

1 . 80000000 E+02 


0. 

0. 

0. 

9.63877017E-02 

0. 

1 .92 152523E-01 

0. 

2 . 86675609E— 01 
0. 

3. 79346129E-0L 
0. 

4.69565224E-01 

0. 

4. 79690999E-01 

0. 

4.88201866E-01 

0. 

4. 94634913 E-01 

0. 

4. 986402 44 E-01 


-1.68684C79E-02 

0. 

— 1. 58774050E-02 
9. 80495849E-02 
— 1.4099 392 2E-0 2 
1.95427712E-01 
-1.432155896-02 
2.91030079E-01 
-1 . 85279504E-02 
3.84684915E-01 
-2.382102 93E-02 
4.7393 90 02E-01 
-1.65083079E-02 
4.85714649E-01. 
-4. I 49 2 6 549 E- 03 
4. 95801299E-0 1 
1. 35943522E-02 
5.04101720E-01 
3 . 89381740E-02 
5. 10577584E-01 


3. 68479746E-1 l 
3. 55271368E-1 5 
9.33799115E-12 
— 3.87734284E-03 
-3.54756207E-12 
-1.5484315 IE— 02 
-5. 32660 76 8E -12 
-3. 474590 96E -02 
-4.78371127E-12 
-6.15376536E-02 
-4.33108709E-12 
-9.56864122E-02 
-2.69725581E-12 
-1.01469577E-01 
-2. 1 111528 6E -12 
- 1.09440 86 8 E-01 
-2. 21048554E -l 2 
-1.191 66721E-01 
-2.82432049E-1 2 
-1.30 118144E-01 


6.69775859E-11 
2. 1320901 l E-02 
1.54283592E-11 
1. 674578 76 E-02 
-8.74578362E-12 
4.70544177E-03 
-1.15550294E-11 
-1.70 462 405E-02 
-9.51693460E-12 
-4. 55153281E-02 
-7.93182332E-12 
-8.54002533E-02 
— 1.01 96 187 8E -11 
-9. 330781 47E -02 
-1.22428436F-11 
-1.03023716E-01 
-1.51611435E-11 
-1.143 87653 E-01 
—2. 15144568E-1 1 
-1.27275840E-01 


1.29678747E+01 
1 . 14504882E+00 
1.36423830E+01 
1. 1 4486 232 E *00 
l • 4369085 4E+ 01 
1. 14472289E+00 
1 . 5 1680872 E+01 
1. 14453026E+00 
1.54686724E+01 
1. 14416846E+00 
1 .43 128 61 3E+01 
1.14364146E+00 
1. 13149745E+01 
1 . 14365 218E+00 
8.81693237E+00 
1. 144 1788 4E i-OO 
6. 85699201E+00 
1.14482646E+00 
5 .0546 1250E+00 
1 . 1457132 1E+00 


1 .3035950 1 E+01 
1.14639301E+00 
1. 3 686 8 202 E+01 
1 • 14568845E *-00 
1.43038014E+01 
1. 1451459 7E+0C 
1.49926844E+01 
1. 14484935E+0C 
1.53838120E+01 
1. 14475275E+00 
1 . 46884446E+0 1 
1.14473095E+00 
1 . 24563883E+0 1 
1 . 14502914E *-00 
1.02829383E+01 
1.14549014E+00 
8.18740790E+00 
1 . 14576509E+00 
5.90531786E+00 
1. 14614066E+00 


8.24937531E-01 
4.38663956F-02 
8. 744 83 463 E-01 
4. 38993945E-02 
9. 278711 47E-01 
4. 39088492E— 02 
9.87511 53 2E-01 
4.39471912E-02 
1.0 12941 60E *-00 
4.41214470E-02 
9. 3 343 365 IE— 01 
4.449 86416E-02 
7. 17855566E-01 
4.47737898E-02 
5. 4 1 040 82 4E— 01 
4.48350592E-02 
4. 067 8 8 38 7E -01 
4.48890719E-02 
2 .878 5366 8E-0 1 
4. 49 7 2 5261 E-02 


8. 19450798E-01 
4.31649360F-02 
8. 70 8 74461 E-01 
4. 346 17120E-02 
9.19467729E-01 
4. 3689381 6E -02 
9.71786586E-01 
4.37891116E-02 
1.00107849EI-00 
4. 381001 42 E-02 
9. 5 08 5995 IE— 01 
4. 3 8 7763 71 E-02 
7. 888 80527E-01 
4. 39294046E-02 
6. 344 65 82 OE— 01 
4. 3957999 3E-02 
4. 91242 50 8 E-01 
4. 4161583 7E-02 
3.407 3893 6E-0 1 
4.451 11 156E-02 


5.06218276E-01 
3.43016968E-01 
5.07701456E-01 
1.78924889F-01 
5.10766877E-01 
5.92008810E-33 
5.16467038F-01 
1.60214803E-01 
5.232 56288E-01 
3.16223424E-01 
5.27821404E-01 
4.92093320E-01 
5.26081704E-01 
7.82764760E-01 
5.24844003E-01 
1.07308937E*-00 
5.24753792E-01 
1. 32965248E+00 
5.25035728F-01 
1. 603349 17E*- 30 


5.05733052E-01 
6. 3195558 7E -01 
5 .06 7625 94E -01 
4.31778807E-01 
5.08615863E-01 
2.05092874E-01 
5 . 12452562E-3 1 
6.62031052E-32 
5. 17199075E-01 
2. 7026186 9F -01 
5.2039 73 31F-01 
4.89036937E-01 
5. 18572873E-01 
7 . 37634555E-01 
5.17545208E-01 
9. 80358700E-D1 
5.17857046E-01 
1. 2029742 9E +00 
5.18311158E-01 
1.46850829E+00 


O) 

CO 



of 300 1 

300 2 

300 3 

300 4 

300 5 

3 00 6 

300 7 

300 8 

300 9 

300 10 


3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


7 -2.428811906-01 
1 .80000000EI-02 
7 -1.840932806-01 
1.80000000E+02 
7 -1.138866856-01 
1.8000000064-02 
7 -3. 296779046-02 
1 . 8000000064-02 
7 5.653636436-02 

i. 800000006+02 
7 1. 437530966-01 

1.800000006+02 
7 2.235746146-01 

1.800000006+02 
7 3.009054826-01 

1 -80000000E+02 
7 3.721262026-01 

1 . 80000000E+02 
7 4 . 440042 8 0E -0 1 

1.800000006+02 


— 3.46192830E-02 

0 . 

-3.548051666-02 
9.97114681E-02 
-3.5C712200E-02 
1.9870 29 01E-01 
-3.633262126-02 
2.953845486-01 
-4.33088130E-02 
3. 900237 02 E- 01 
-5.248935916-02 
4. 7831278 OE -01 
-4. 20816213 E- 02 
4.91738299E-01 
— 1.68936884E-02 
5.034007326-01 
2.410026476-02 
5.13568526E-01 
8.578663496-02 
5.225149246-01 


9.962290086-11 
4.264180226-02 
4.84493848E-11 
3.73689181E-02 
2.279576526-11 
2.48951986E-02 
1.926200306-11 
6.53428681E-04 
2.11 600490E-1 1 
-2. 9 49 3002 6E -02 
2.54956362E-11 
-7.51140943E-02 
2.953416206-11 
-8.51460522E-02 
3. 2 430 73 81E-1 1 
-9.660656316-02 
3.50 182942E- 1 1 
-1.096085856-01 
3. 856757 94E-11 
- 1 . 24433536E-0 1 


1.338 45 359E+01 
1. 14891828E+00 
1 • 3 90 2 404 9 E+ 01 
1 . 14739080E+00 
1.434362016+01 
1. 14611149E+00 
1. 4808806 2E+01 
1. 1453373 9E+00 
1 . 5 12 86279E+0 1 
1.14531204E+00 
1.489141666+01 
1.146C7923E+00 
1 . 33869 802E+0 1 
1 • 14735004E+00 
1.17929177E+01 
1.14928924E+00 
1.01412857E+01 
1.15181 069E+00 
8.24968624E+00 
1.15573546E+00 


8.24488422E-01 
4. 18609824E— 02 
8.72245523E-01 
4. 25664940E-02 
9. 13975248E-01 
4.31791930E-02 
9.54045830E-01 
4.35440544E-02 
9.77312885E-01 
4. 353063 54E-02 
9.53280172F-01 
4. 31462401E-02 
8.366294326-01 
4.264162926-02 
7. 138512856-01 
4. 18798347E— 02 
5.902698926-01 
4. 094698 5 2E— 02 
4. 535 50856E-0 1 
3.95791768E-02 


5. 0534515 5E -01 
9.573185646-01 
5. 059 90 52 5F- 31 
7.23O41451E-01 
5.069252876-31 
4.549311326-01 
5 .08901265E -3 1 
1.859436066-01 
5. 1 1698393E-31 
2.694738886-01 
5. 139493726-01 
5. 82606200 E- 01 
5. 125C5229E-01 
8. 7939 47 18E -01 
5.114162596-31 
1. 188323226+00 
5.111367756-31 
1.505641556+00 
5.102489476-01 
1. 8881742 4E +30 


300 1 

300 2 

300 3 

300 4 

300 5 

300 6 

300 7 

300 8 

300 9 

300 10 


4 


4 


4 


4 


4 


4 


4 


4 


4 


4 


7 -3.22915611E-01 
1 . 80000000E + 02 
6.39627034E-02 
7 -2 .56 51733 IE-01 
1.80000000E+02 
6. 2069946 9 E-02 
7 -1.86175303E-01 
1.800000006+02 
6.136 104876-02 
7 -9 . 2 72 6959 4E— 02 
1.80 0000 00E+02 
5.46823302E-02 
7 4.42075226E-02 

1.800000006+02 
5.06651554E-02 
7 2.49014456E-01 

1.800000006+02 
3.35322980E-02 
7 2.873 06 480E-01 

1.800000006+02 
3.043170246-02 
7 3 . 93857 31 5E-0 1 

1 .'80000000E+02 
2.98392536E-02 
7 4.811475786-01 

1.800000006+02 
3 . 181 41357E-02 
7 5.419787036-01 

1.80000000E+02 
3.68131635E-02 


-4.98713933E-02 

0. 

-8.782278226-04 
-5 . 67000820E-02 
1. 01 37 3351 E- 01 
-5 . 63973263E-04 
-6.15639481E-02 
2. 019780 90E-01 
-2. 78605136E-04 
-6. 36672157 £-02 
2.9973 9018E-01 
4.261176416-05 
-7.349302146-02 
3 .9536248 8E-01 
6.584166026-04 
-1.20333295E-01 
4.82686558E-01 
1.46817515E-02 
-7.846179276-02 
4.977619496-01 
1. 976781266-03 
-3.298682436-02 
5.11000165E-01 
1.082424556-03 
4.071501446-02 
5.230353326-01 
1.303066196-03 
1.377729596-01 
5 . 34452264E-0 1 
1 . 5 6 54 l 5 34 E-0 3 


1.324553146-10 
6.39627034E-02 
4.20764005E-03 
1 . 53658208E-1 0 
5.79920486E-02 
-3.39975109E-02 
1.53658208E-10 
4.50849554E-02 
-2. 7 6 5097 27 E-0 2 
1 . 53658208E-1 0 
1.83530978E-02 
-3. 8248301 IE-02 
1 • 5 3658208E-1 0 
-1.34706771 E-02 
-8.21538937E-02 
1 . 53658208E-10 
-6.48279353E-02 
-1 .80518371E-02 
1.53658208E-10 
-7. 69 842 8 96 E-0 2 
— 5.71401835E-03 
1.53658203E-10 
-9.01894106E-02 
1. 13340757E-02 
1 .53658208E-10 
-1.04829517E-01 
3.04169042E-02 
1 . 53658208E-10 
-1.21591232E-01 
5. 487747 19E-02 


1.39433576E+01 
l. 15268699E+00 
0. 

1 .4 1 799 769E +0 l 
1. 15002265E+00 

0. 

1 . 43685590E+0 l 
1. 14796334E+00 
0 . 

1.45259021E+01 
1 .146287116+00 
0. 

1.456393736+01 

1.145887606+00 

0. 

I • 43086927E+0 l 
1 . 14861C88E+00 
0. 

1.40868433E+01 
1. 15106002E+00 
0. 

1.36532587E+01 
1. 156081566+00 
0. 

1 .3158691 8E+0 1 
1.16221690E+00 
0. 

1.2636 3999E+01 
1.1 692 1495E+00 
0. 


8.33603585E-01 
3.998 17873E— 02 

8.70289396E-01 
4. 122 97 552 E-02 

9.0008902 1E-0 l 
4. 2229653 8E -02 

9. 2 54 l 1633 F-0 1 
4.30710332E-02 

9. 316 0404 7E-01 
4.3275751 3E-02 

8.90570259E-01 

4.19115289E-02 

8.55767128E-01 
4. 073 8225 1 E-0 2 

7. 89517682 E- 01 
3.R4535184E-02 

7. 162 79160E-01 
3.58355759E-02 

6. 40 84 269 IE- 01 
3.301 54281E-02 


5.0471 8488E-01 
1. 2994041 4E +00 

5.04984303E-01 

1.02787904E+00 

5 .0523672 IE-01 
7.572146296-01 

5.05531220E-01 

4.29520021E-01 

s.ofeicssiae^oi 

3.266044166-01 

5.19294515E-01 

1.072461506+00 

5.072540226-01 
1 • 1 72 78874E + 00 

5 . 06419994E -01 
1.60348372E+00 

5.06576843E-01 

2.02859165E+00 

5.06776971E-01 
2. 440062 53 E+ 00 
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